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XX. BAND DRITTES HEFT 1952 


Die Stabilitat erzwungener harmonischer Schwingungen gerader 
I-Trager im Verband eines Tragwerkes*. 


Von K. A. Reckling. 


1, Einleitung. Die bisher itber Stabilitatsprobleme auf dem Gebiete der Elastokinetik ver- 
offentlichten Arbeiten gehen in ihrer mathematischen Behandlung auf zwei verschiedenen 
Wegen vor: 

1. In der von E. Mettler entwickelten ,,Allgemeinen Theorie der Stabilitat der elastischen 
Bewegung‘! sowie in einer Reihe von Einzelarbeiten? wird fiir die ».Nachbarverschiebungen‘, 
welche fiir die Stabilitatsbeurteilung maBgebend sind, ein Ritz-Ansatz mit zeitabhangigen Koeffi- 
zienten eingefiihrt, der die jeweiligen Randbedingungen erfiillen mu8. Man gelangt dadurch zu 
einem System von unendlich vielen gewéhnlichen linearen Differentialgleichungen mit perio- 
dischen Koeffizienten, das in einigen konkreten Einzelfallen in mehrere getrennte Systeme oder 
auch Einzelgleichungen (Mathieu-Gleichungen) zerfallt. Die Integrale dieses Systems werden in 
trigonometrische Reihen und alle vorkommenden Koeffizienten mit der Methode der ,,Stérungs- 
rechnung“ nach Potenzen eines kleinen Parameters entwickelt, der die Amplitude der erregenden 
Kraft angibt. Auf diesem Wege ergeben sich zusammengehérige Werte von Erregerfrequenzen 
und Belastungsamplituden, welche bei beliebigen Anfangshedingungen zu unbegrenzt anwachsen- 
den Nachbarbewegungen fiihren kénnen. 

2. Dagegen geht W. Kucharski® in einer Arbeit itber die durch pulsierenden Schub belastete 
Rechteckplatte ganz anders vor, indem er die Methode der ,,Stérungsrechnung“ auf eine partielle 
Differentialgleichung iibertragt und bereits die Nachbarverschiebungen selbst nach wachsenden 
Potenzen eines kleinen Parameters entwickelt, wobei die Entwicklungs-Koeffizienten Funktionen 
der Zeit sind. 

Die vorliegende Arbeit will die durch pulsierende Belastungen mit vorgegebener Erreger- 
frequenz w erzwungenen stationaren Schwingungen eines I-Tragers im Verband eines Tragwerkes 
hinsichtlich ihrer Stabilitat untersuchen. Dadurch, daB Nachbarbauteile an den Enden des 
Tragers angeschlossen sind, muB man dort auferdem noch elastische, zeitabhangige Einspann- 
bedingungen beriicksichtigen. Dieses Problem wurde bisher nur von S. Woinowsky- Krieger? fiir 
den speziellen Fall der Belastung eines an den Enden frei aufliegenden Tragers durch gleichgroBe 
Endmomente nach der erstgenannten Methode geliést. Es soll im folgenden in gréBerer Allge- 
meinheit mit dem zweiten Verfahren untersucht werden, zu dessen Ausdehnung auf das hier 
vorliegende System zweier partieller Differentialgleichungen sich der Verfasser besonders deshalb 
entschloB, weil einige Vorversuche auf dem ersten Wege sehr umfangreiche — vor allem durch 
die Einspannung an den Tragerenden bedingte — Rechnungen in Aussicht stellten. 

Die Gestalt der fiir den allgemeinen Belastungsfall zu integrierenden Differentialgleichungen 
der Nachbarverschiebungen legte es nahe, diesen naiherungsweise auf den einfacheren Fall der 
Belastung durch gleich groBe pulsierende Endmomente zuriickzufiihren; fiir diesen wird dann 
die eigentliche Stabilitatsrechnung durchgefiihrt. Durch die Entwicklung der Nachbarverschie- 
bungen nach einem kleinen Parameter ¢, der die Belastungsamplitude miBt, erhalt man unendlich 
viele, sukzessiv aufzulisende partielle Differentialgleichungen. Diese lassen eine andere mathe- 
matische Behandlung zu als die entsprechenden Gleichungen der Arbeit von W. Kucharski?, 


* Gekiirzte Fassung einer von der Technischen Universitat Berlin-Charlottenburg genehmigten Disser- 
tation. Der Verfasser ist den beiden Berichtern, Herrn Prof. Dr.-Ing. I. Szabé und Herrn Prof. W. Kucharski 
fiir wertvolle Ratschlage sehr zu Dank verpflichtet. 

1 E. Mettler, Ing.-Arch. 16 (1947), S. 135; sowie 17 (1949), S. 418. 

2 E. Mettler, Mitt. Forsch.-Anst. G.H.H.-Konzern 8 (1940), S.1; E. Mettler, Forschungsh. a. d. Gebiete 
d. Stahlbaus, Heft 4; S. Woinowsky-Krieger, Ing.-Arch. 13 (1942), S.197; K. Kewjal, Dissertation 
Hannover 1949; F. Weidenhammer, Ing.-Arch. 19 (1951), 5S. 162. 

3 W. Kucharski, Ing.-Arch. 18 (1950), S. 385. 
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wodurch gewisse dort auftretende Schwierigkeiten umgangen werden kénnen und man auf eine 
Reihe von Randwertaufgaben fiir gewdhnliche inhomogene Differentialgleichungen mit kon- 
stanten Koeffizienten gefiithrt wird; nach deren Lésung kann man ,,Instabilitatsgebiete“ in einer 
(c, w)-Ebene angeben. Die kritischen Erregerfrequenzen, bei denen Instabilitat auftritt, kénnen 
aber auch ohne Lisung dieser Randwertaufgaben in einfacher Weise allgemein angegeben werden. 

An Hand von zwei Beispielen mit angefiigter Zahlenrechnung fiir konkrete Einzelfalle wird 
die praktische Brauchbarkeit des Verfahrens dargestellt. Am Beispiel eines axial pulsierend 
belasteten Stabes wird schlieBlich gezeigt, daB der hier eingeschlagene Weg auch in den bisher 
mit der erstgenannten Methode untersuchten Einzelfallen schnell zu denselben Ergebnissen fiihrt 
wie dort. 


2. Spezielle Voraussetzungen und Naheres zur Aufgabenstellung, Wir betrachten in Abb. 1 
einen geraden, zur y- und z-Achse symmetrischen I-Trager von konstanter Biegesteifigkeit, der 
aus einem durch harmonisch schwingende 

Belastung Krafte einer bestimmten vorgegebenen Fre- 

quenz w beanspruchten Tragwerk heraus- 
geschnitten sein mége. Wir setzen voraus, 
daB die Eigenschwingungen des Tragwerkes 
abgeklungen sind, so daB nur der stationare 
Anteil der erzwungenen Schwingungen mab- 
gebend ist. Der Trager mége auBerdem noch 
durch eine in der Steg-Ebene wirkende gleich- 
maBige, zeitlich verdnderliche Belastung p(t) 
von der gleichen Frequenz w beansprucht 
werden. Die in der (x,y)-Ebene wirkenden 
Schnittkrafte an den Stabenden sind dann 
harmonisch mit der Frequenz w veranderliche 
Krafte. Alle diese Schwingungen sollen die 
gleiche Phase haben. Dann wirken in der 


Nachbarbewegungszustand (x, y)-Ebene die folgenden vorgegebenen ein- 
gepragten Krafte bzw. Momente: 
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Abb. 1. Belastung des Triigers mit den zugehérigen Abb, 2. Langenelement der Stabachse 
Verschiebungen. im Grundbewegungszustand, 
P(t)= P+ Picoswt, p(t) = py t+ pcos wt, | 


Q:(t)=Qio+ Qi: cost, M,(t)= Miy+ My, coswt (i=1,2), | (1) 


wo der Zeiger 0 die statische Vorbelastung bezeichnet. 


Der durch diese Belastungen erzwungene stationare ,,Grundbewegungszustand“ in der (x, v)- 
Ebene setzt sich zusammen aus Langsschwingungen u(x,t) in Richtung der x-Achse und aus 
Biegeschwingungen v(x,t) in der (x,y)-Ebene. Diesem Grundbewegungszustand wird sodann 
ein durch eine beliebige kleine Stérung verursachter ».Nachbarbewegungszustand“ iiberlagert, 
der durch die Auslenkung w(x,t) in der (x,z)-Ebene und den Torsionswinkel O(x,t) gekenn- 
zeichnet wird. Da diese beiden letzten Verschiebungen miteinander gekoppelt sind, wird der 
Nachbarbewegungszustand als ,,Kippschwingung“ bezeichnet. Der raumliche Nachbarzustand 
soll sich vom ebenen Grundzustand nur um kleine Gréfen erster Ordnung unterscheiden. 

Da an den Tragerenden Nachbarbauteile angeschlossen sind, wirken dort noch zeitlich ver- 
ainderliche Einspannmomente My,;(t), deren Vektoren parallel zur y-Achse weisen. Der Trager 
sei dort auBerdem so gelagert, daB die Achse des kleinsten Tragheitsmomentes (d.h. der Trager- 


XX. Band 1952 Reckling: Die Stabilitat erzwungener harmonischer Schwingungen usw. 139 


steg) stets in der (x,y)-Ebene bleibt. Samtliche eingepriigten Krafte und Momentenvektoren 
sollen ihre Richtung bei der Verformung stets beibehalten. 

Der Einflu8 der Rotationstriagheit um die Achse des kleinsten und grifBten Tragheitsmomentes 
wird von vornherein vernachlassigt. Lassen wir schlieBlich auch die Dampfung auBer acht, so 
bleiben wir bei der Beurteilung der dynamischen Stabilitat bestimmt auf der sicheren Seite?. 


3. Die grundlegenden Bewegungsdifferentialgleichungen und ihre Randbedingungen; Einfithrung 
von Naherungen. Im Grundbewegungszustand wollen wir einem Stabquerschnitt die langs der 
ausgebogenen Stabachse gerechnete Koordinate x und die dazu senkrechten, mit den Haupt- 
tragheitsachsen zusammenfallenden Koordinaten y, z zuordnen; der ausgekippte Nachbarzustand 
werde durch die entsprechenden kérperfesten Koordinaten x, y, 2 gekennzeichnet (Abb. 1). 
Da beim Ubergang zum ausgekippten Bewegungszustand keine zusatzlichen Lingskriifte auf- 
treten, dandert sich die Linge der Stabachse nicht, und es ist A4x—Ax. Ferner soll naherungs- 
weise 1x = Ax gesetzt werden; denn nach 
Abb.2 ist das Lingenelement der aus- 
gebogenen Stabachse im Grundzustand 


Ax = | (2 “ ee a “Ax, 


wenn wir bei einer Entwicklung von Au 
und Av in eine Taylor-Reihe nur das je- 
weils erste Glied mitnehmen. Nun hat das 
Glied du/éx = ¢,=o,/E fiir Stahl mit 
E ~2-10°kg/cm? die GréBenordnung 10~°, 
Wenn wir die gréBte Tragerdurchsenkung 
Umax =f, wie es in der Praxis iiblich ist, 
durch die Forderung f < 1/500 begrenzen 
und eine Durchbiegungskurve von der 
Form v= f sin (zx/l) annehmen, so wird 


2 
andererseits 2 < 
Ox 


2 
ae und e) gegeniiber ,,Kius‘* vernach- 
lassigen kénnen. 

Die im Grundzustand an der Stelle x 
wirkenden Krafte und Momente sind von 
endlicher GréBe und rufen endlich grobe 
Verschiebungen in der (x, y)-Ebene hervor. Abb. 3. Krifte und Momente am ausgekippten Trigerelement. 
Die beim Ubergang in die ausgekippte 
Nachbarlage hinzutretenden Krafte und Momente Q., M,, M, sowie die durch sie bewirkten 
Verschiebungen w und # und deren Ableitungen sind dagegen nach Voraussetzung von erster 
Ordnung klein. 

Wenn alle am ausgekippten Tragerelement von der Linge Ax wirkenden Schnittkrafte und 
-momente einschlieBlich der negativen Massenbeschleunigungen ins Gleichgewicht gesetzt werden 
(Abb. 3), so ergeben sich nach Fortlassen von kleinen GréBen zweiter Ordnung und nach Division 
durch 4x beim Ubergang zur Grenze Ax -—> 0 die Gleichungen 


m2 


5 TI! ; 
3 10-4, so daB wir 


N'+ Mzg' = — oii, | 

— Ng’ + M+ p=ev, (2) 
Ny’ — My — 2 M79’ — M,9" — (M.q') + pO = ow, | 

M,y' + M,y' — M;— ped = — 01,8, 


worin die Ableitungen nach x mit Strichen und die nach t mit Punkten bezeichnet sind und die 
Konstanten 9 die schwingende Masse des Tragers je Langeneinheit, i, den Tragheitshalbmesser 
des Tragerquerschnitts, e die Entfernung des Lastangriffs von p von der Stabachse bedeuten. 


1 Uber den stabilitatsverbessernden EinfluB der Dampfung vgl. die in der FuBnote 2 von 8S. 137 an- 
gefiihrten Arbeiten von Mettler und Weidenhammer. 
IL Lie 
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Es gelten nun folgende Elastizitatsgesetze fiir die SchnittgréBen: Fiir die Langskrafte hat man 
N=— EF uw’. (3) 

Den Kriimmungen 
y = py’ bzw. y = w’! + v' (4) 


der neutralen (x,z)- bzw. (x,¥)-Ebene im ausgekippten Zustand! sind die entsprechenden Biege- 
momente M, bzw. M, proportional: 


M,~— B,v", M,~ B,(w’’ + v8) (5) 


mit B, = EI, als Biegesteifigkeit bezogen auf die z-Achse, B, = EI, als Biegesteifigkeit bezogen 
auf die y-Achse. Zu dem de Saint-Venantschen Anteil M,, = C(0’ —v w’) am Torsions- 
moment M, (mit C= GI,= Torsionssteifigkeit des Querschnitts)' kommt infolge einer zu- 


sitzlichen Ausbiegung w= = h® der Flanschachsen ein zusatzliches Torsionsmoment 
1 ttt 
M,. = — ve WB, 
hinzu?, Aus der Zusammensetzung dieser beiden Anteile folgt 
, TENE: h? B gy” (6) 
M, = C(#' — v ©) rare y : 


Einer Integration der sich aus (2) unter Beriicksichtigung von (3) bis (6) ergebenden gekoppel- 
ten Differentialgleichungen stellen sich nun sehr groBe Schwierigkeiten entgegen, eo) daB wir ge- 
zwungen sind, eine Naherungsannahme zuzulassen: Nach der oben erfolgten Abschatzung wollen 
wir die Kriimmung q’ = v’’ <7?/5001 als kleine GréBe erster Ordnung annehmen, wie das 
auch bei der Berechnung der statischen Kippstabilitat allgemein iiblich ist ?. Dann kénnen wir 
die Glieder Mz’, M,y’ und M,q’ in Gleichung (2) fortlassen und auBerdem die Glieder v’’ 3 
gegeniiber w’’ in (4) und v’’ w’ gegeniiber # in (6) als kleine Groen zweiter Ordnung vernach- 
lassigen. Das Glied Ng’ in der zweiten Gleichung von (2), durch das der Einflu8 der Langskrafte 
auf die Biegeschwingungen in der (x, y)-Ebene in die Rechnung eingeht, wollen wir vorerst noch 
nicht fortlassen; denn die Langskrafte kénnen fiir gewisse Erregerfrequenzen zur Instabilitat der 
Biegeschwingungen u(x,t) fiihren, selbst wenn ihr Einflu8 auf diese wegen der Kleinen Kriimmung 
sehr gering ist; wir werden darauf in Abschnitt 4b noch zuriickkommen. Die Gleichungen (2) 
gehen dann iiber in 


EFu’ =i, (a) 

EFu'v"’ — BY + p=ov, (b) 

EF uw’ — Byw!Y + B,(20'"# + v8") + pO8=0w, (c)f (7) 
2 a 

Bow" —™ BO + CO" + ped = itd. (a) 


Die beiden ersten gekoppelten Differentialgleichungen beschreiben die Grundbewegung [Langs- 
schwingungen u(x,t) bzw. Biegeschwingungen v(x,t) in der (x,y)-Ebene]. Die durch w(x,t) 
und (x,t) beschriebene Nachbarbewegung ist durch Einsetzen der Integrale u(x,t), v(x,t) der 
beiden ersten Differentialgleichungen in die beiden letzten zu ermitteln. 

Zu den Gleichungen (7) treten noch Einspannbedingungen in der (x, z)-Ebene fiir die Trager- 
enden hinzu, wo Nachbartrager angeschlossen sein migen, die mit zeitlich veranderlichen Ein- 
spannmomenten Mp ,(t)= B,w'’(0,t) und Mp (t) = B,w’’ (1,1) auf die Tragerenden zuriickwirken 
(vgl. Abb.1b). Diese Einspannmomente wachsen bei elastischer Einspannung proportional zur 
Neigung der Biegelinie an den Enden: 


Mri (t) = cy (0,t) = c,w’ (0,t), Mr 2(t) = — eny(l,t) = — c,w’ (1,t) , 


wobei ¢,, c, die Federkonstanten der elastischen Einspannung sind. Durch Gleichsetzen der 
Ausdriicke fiir die Einspannmomente erhalten wir 


c,w’ (0,t) — Bewen{0.t) == 0 sind cw’ (1,t) + Bow (La 0. 


Dazu tritt noch die Forderung, daB die Durchsenkungen w an den Tragerenden verschwinden 


* Vel. z. B. Pfliiger, Stabilitatsprobleme der Elastomechanik, S. 155. 
2 Vel. z. B. E. Chwalla, Forschungshefte aus dem Gebiete des Stahlbaus, Heft 2. 
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sollen, so daB wir als Randbedingungen 


w(0,t) = w(I,t) = w’ (0,t) — C,w'’(0,t) = w’ (I,t) + C,w’’(1,t) = 0 

mit | (8) 

Eee ee By | 

ae re: und Cia at | 

erhalten. [Fir freie Auflagerung an den Enden ist beispielsweise c,—=c,—0 und damit 

w'"(0) = w"’(I); bei vollkommener Einspannung miissen wegen 1/c, = 1/c,—=0 die Randbedingun- 
gen w’(0) = w’(l) = 0 erfiillt werden.] 

Bei der Aufstellung der Randbedingungen fiir den Torsionswinkel # beachten wir, daB sich die 


Ausbiegung der Flanschachsen um den Betrag + w(x,t) = ++ = h&(x,t) von der Ausbiegung w 


der Tragerachse unterscheidet (vgl. Abb.1), und daB dann die Tangentenneigung der Flansch- 
achsen an den Tragerenden (w-+- w)’ =w’ + S h® betragt, wodurch sich die Endquerschnitt- 
flachen des Tragers verwélben. Die elastische Einspannung durch die Nachbartrager verhindert 
diese Verwélbungen teilweise. Es erscheint naheliegend, fiir den Verlauf der zusatzlichen Aus- 
biegungen w der Flanschachsen dieselben Randbedingungen wie fiir die Ausbiegung w der Trager- 
achsen anzunehmen. Dann gelten aber wegen h = const fiir } dieselben Randbedingungen (8) 
wie fiir w. 


4, Der Grundbewegungszustand, a) Longitudinalschwingungen. Da nur die stationaren 
Schwingungen des Tragwerkes interessieren, kénnten wir fiir u(x,t) eine harmonische Bewegung 
mit der vorgegebenen Frequenz q@ ansetzen u(x,t) = uy(x) + u,(x) cos wt, wo up(x) die Ver- 
schiebung infolge der zeitlich unveranderlichen Lastanteile ist. Durch Einfiihren in (7a) wiirden 
wir zwei gewohnliche Differentialgleichungen zweiter Ordnung fiir die beiden Koordinaten- 
Funktionen uy(*) und u,(x) erhalten, zu denen infolge der an den Stabenden wirkenden Lings- 
krafte je zwei Randbedingungen mit zum Teil zeitabhangigen Koeffizienten hinzutreten. Die so 
gewonnene Lisung wiirde jedoch fiir die weitere Rechnung zu unhandlich werden, so daB wir 
u(x,t) auf einem anderen Wege bestimmen wollen: 

Wir gehen dazu vom Elastizitatsgesetz Gleichung (3) fiir die Langsverschiebungen u(x,t) 
aus und beachten, da die in einem Querschnitt x iibertragene Langskraft die GriBe P(t) zu- 


satzlich einer von der Tragheitswirkung herriithrenden Kraft N(x,t) hat: 

EFu' = — [P,-+ P, coswt+ N(x,t)]. (9) 
Wir wollen nun u(x,t) nach den Eigenschwingungsformen cos kz-x/l des um die Stabmitte 
symmetrisch schwingenden Stabes entwickeln: 


u(x,t) = >" u(t) cos *2*, 
= Dy Dieter 


Die an der Stelle x infolge der Tragheitswirkung des gesamten Stabes wirkende Kraft wird 


ae = eee en ol = 1 du, . kax 
N(x,t) = — 0 Be i Jess a és) = — =— oy i de 88 
0 


L 
k=1,3)5 000 | = p35 vse 


Nun entwickeln wir auch die am Ende angreifende Langskraft in eine nach sin kx x/l fort- 
schreitende Fourier-Reihe 
4 ee ee 22 
Po-+ P, cos wt me (P,-+ P, cos ee 2 =, sin - 


und erhalten mit diesen Entwicklungen aus (9) die Differentialgleichungen 
ae whem = Coen a helices 54 ct) 
mit 
ee i der Longitudinalschwingungen 
1, = = siieat k-ter Eigenfrequenz der Longitudinalsc gungen. 
Die Lésungen der homogenen Gleichungen sind nicht von Interesse, da die freien Schwingungen 
nach Voraussetzung abgeklungen sein sollen. Aus den partikuléren Lisungen der inhomogenen 
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Gleichungen erhalten wir 


/ oo cos bes P,cos at 

1 

u(x, t) hs ar ae Py rae 
PERS vane aoa | 


als Gesamtlésung in einer Form, die es gestattet, eine wesentliche Vereinfachung vorzunehmen : 
Wir wollen die Untersuchung namlich nur fiir solche Erregerfrequenzen durchfiihren, fitr die 
w <j, und damit erst recht @ <a, ist. Folgende Abschatzung zeigt die Berechtigung dieses 
Naherungsschrittes: Aus spiteren Rechnungen ergibt sich, da die kritischen Erregerfrequenzen oO) 
héchstens die GréBe der doppelten Biegeeigenfrequenz 2 w,, in der (x,z2)-Ebene oder die der 
doppelten Eigenfrequenz der Torsionsschwingungen 2 ,, haben. Fiir den in der (x,z)-Ebene 
an den Stabenden frei aufliegenden Trager sind diese Frequenzen nach Spaterem [vgl. Gl. (32) | 


na\" 1 /ETy ie hte.) GLa , mah? Ely 
2 iden = 2 a / oa bzwa ec aye ; (} tGhP ): 
Nehmen die Erregerfrequenzen @ diese Werte an, dann werden die (w/);)? mit G ~ 0,4 E und 
I, ~ F3/40 i, naherungsweise’ héchstens gleich 
DUC Up Ne ee 2 nzhipty\? 
nt ea oder’ 42) 10=4 74 ia L ( Fl £ ") | 
und kénnen fiir kleine n und Trager mit normalem I/h gegeniiber ,,eins“‘ vernachlassigt werden. 
Das gilt erst recht fiir die Quotienten (w/a,,)?, so daB wir als Naherungslésung der Gleichung (7a) 


Al ~ il kux 
u(x,t) = = (Po+ Pyeoswt) > 7, c08—7— 
| Ji BE aan 
erhalten und daher 
EF u' = — (P,-+ P, cos wt) = — P(t) (10) 


in (7b) und (7c) einfiihren kénnen. Wir hatten diese Naherung auch aus (7a) direkt gewinnen 
kénnen, wenn wir dort ii=0 gesetzt hatten; das wiirde bedeuten, daB EF u’ von x unabhangig 
ist und nach dem Elastizitatsgesetz (3) fiir x= 0, also auch fiir alle x den Wert — P(t) annimmt, 
wie es auch die Gleichung (10) aussagt. 


b) Transversalschwingungen in der Lastebene. Gleichung (7b) lautet mit (10) 
B,v'V + P(t)v’’ — p(t) = — ov. (11) 


Bevor wir diese Gleichung in dhnlicher Weise wie unter a) integrieren, verdient zundchst noch 
ein wichtiger Sonderfall kurz erwahnt zu werden: Falls p(t)=0 ist und die Momente an den 
Tragerenden verschwinden, haben wir die Differentialgleichung der Biegeschwingungen des durch 
die pulsierende Axiallast P(t) beanspruchten Stabes; dann beschreibt Gleichung (7a) die Grund- 
bewegung u, und (11) ist die Differentialgleichung der Nachbarbewegung, die bereits fiir ver- 
schiedene Einspannbedingungen an den Enden auf ihre Stabilitat untersucht wurde?. In diesem 
Fall ist iibrigens wegen B, <B, viel eher eine instabile Nachbarbewegung w(x,t) in der (x,2)- 
Ebene zu erwarten. Wir erhalten auch tatsadchlich mit v(x,t)=0 aus (7c) unter Beachtung 
von (10) eine Differentialgleichung fiir w(x,t), die derjenigen fiir die Nachbarbewegung in der 
(x, y)-Ebene genau entspricht: 


B,w'V + P(t)w’’ = — ow. (12) 
Wir werden auf diese Differentialgleichung spater noch zuriickkommen. 

Zur Integration von (11) gehen wir vom Elastizititsgesetz (5) fiir das Moment M,(x,t) aus, 
das sich nach Abb.4 aus folgenden Anteilen zusammensetzt: Aus einem Moment M,(x, 1) infolge 
der Belastung p und der durch p hervorgerufenen Querkraft Q,(t) an den Enden, und aus zwei 
durch die pulsierenden Endmomente M,(t) und M,(t) hervorgerufenen Anteilen, wobei im Falle 


M,(t) + M,(t) noch eine Querkraft [M,(t)— M,(t)] 1/1 an den Enden hinzukommt; ferner aus 
einem von der Langskraft herrithrenden Anteil P(t)-v und schlieBlich aus dem durch die Trig- 


* Diese Naherungsformel fiir Iy nach Féppl, Drang und Zwang II, S. 101, reicht fiir die Abschatzung 


aus, 


2 Vel. FuBnote 2 von S, 137. 
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heitswirkung hervorgerufenen Momentenanteil M(x,t). Es gilt also 


— Bo = M,(x,1) = M, (x1) + My) + [M,) — MQ) 2 + PQ) o(x,2) + Mix). (13) 
Wir entwickeln nun die Ausbiegung v(x,t) in eine nach den Kigenschwingungsformen sin kz x/I 
des an den Enden gestiitzten Tragers fortschreitende Fourier-Reihe: 


foo) 


v(x,t) = » v,(t) sin sae 


I9= NPA aA 


Zur Ermittlung der unbekannten Koeffizienten v(t) entwickeln wir alle auf der rechten Seite 
von (13) stehenden Momente ebenfalls in Fourier-Reihen: 


M, (x1) =} (x — 2%) p(t) = SEBO dgigtte 


me k=1, 3,5 # l 
4 td k 
M,(t) = — M,(t) > Sipe: ae ag 
(SSN Baaio.c 
2 Sab a le 
[M,()— M,()] = — = (M,(@) — MQ) >) SY" sin *2*, 
(p= Ty Pon 


Zur Entwicklung des Momentes M(x,t) beachten wir, daB wir die Tragheitskraft 


co 


—gi(x,)=— 9 >" a(t) sin 274 
b=, 2. «: 
als eine Summe aus kontinuierlich verteilten Einzelbelastungen q),(x,t) = — @ sin (kx/l) %,(t) 


darstellen kénnen. Das gesamte von der Tragheitswirkung herriihrende Moment wird dann 


F . | | [ea Eee ee 
M (x,t) = — em fl faces * dx|= —o EI ay ze ve (?) sin : 
k= ipoe LEB ine 

en) 


Wenn wir alle diese Entwicklungen in (13) einfiihren, ergeben sich k gewohnliche Differential- 
gleichungen zweiter Ordnung fiir die v;,(t) 


2p ; 5 5 4 (EY (y+ L,;, cos wt) (ho), 335 De ah 
v TU 
gat [at — 2 (BY Preos or eet 


el? 


(14) 


(—1)[M)—M,@)] (k= 2.4.6, ...,| 
wobei zur Abkiirzung 


2 
te 4 i 4 


eo (kt)? Piv kx 


2 : 
ee (eee Ni) EG 101) 


peg tee ee 
eingefiihrt wurde und durch @,,—= k ie ie (k?— ep) die k-te Eigenfrequenz der Biege- 


schwingung des durch die konstante Langskraft P,—= epPp, = epB, (x/l)? belasteten Tragers 
in der (x, y)-Ebene angegeben wird; ep ist dabei das : 
Verhaltnis der Last P, zur ersten kritischen Knick- 4) erat fede eae 1 edhe { y “lt 
last Pry. ' 7 a : fete. yy 
Wir haben damit k inhomogene Mathieusche By ae == SN 
Differentialgleichungen zur Bestimmung der Grund- p(t) ihe lant) 
bewegung v(x,t) in der (x,y)-Ebene. Mit der Inte- | ‘ Pee 
gration dieser Gleichungen wollen wir uns hier Lule) el Leto) | 
nicht beschaftigen, weil schon Mettler’ eine ahn- 
liche Differentialgleichung integrierte. Es ergab cae ; 
sich dabei folgendes: Die homogenen Gleichungen (14) beschrieben die Biegeschwingungen 
des durch pulsierende Axiallast beanspruchten, an den Enden freien Tragers und kénnen fiir 
die kritischen Erregerfrequenzen w= 2,,/r (mit r als ganzer Zahl) zur Instabilitat fiihren. 


Abb. 4. Belastung des Tragers. 


1 Vel. die an zweiter Stelle in FuBnote 2 von S. 137 aufgefiihrte Arbeit. 


144 Reckling: Die Stabilitat erzwungener harmonischer Schwingungen usw. Ingenieur-Archiv 


Partikulare Integrale der inhomogenen Gleichungen fithren dagegen zu einer Resonanz im 
iiblichen Sinne, wenn w in der Nahe von ,, liegt. Wir wollen solche Werte von , fiir die bereits 
die Grundbewegung instabil ist oder Resonanzverhalten zeigt, aus unseren weiteren Uber- 
legungen ausschlieBen. AuSerhalb dieser Instabilitats- bzw. Resonanzbereiche, in denen grofere 
Auslenkungen v auftreten, wollen wir fortan den Beitrag P, cos wt: v zum Biegemoment M,(x, t) 
als Nebeneinflu® vernachlassigen und erhalten als Naherungen die Gleichungen (14) ohne die 
Glieder mit P, auf der linken Seite. Die Berechtigung dieser Annahme geht aus den Bemer- 
kungen hervor, die im Abschnitt 3 zu der dort eingefiihrten Naherung (sehr kleine Kriimmung v’’) 
gemacht wurden. Die Lésungen der homogenen Gleichungen interessieren nicht, da sie nach 
Voraussetzung abgeklungen sein sollen; fiir die v,(t) kénnen wir dann eine harmonische Bewegung 
v,(t) = v9, + v1, cos wt mit der vorgegebenen Frequenz w ansetzen und erhalten aus (14) nach 
kurzer Zwischenrechnung 


212 | ~ 2sin 5" Lae M cos wt (ee M | +e 
Oa | > k (k2 — ep) (es Pot 10-4 ; Cay ae PH “yl 
eee rate, 
, (15) 
Oe Shean 
l Ju M,, 4 (M,, — My,,) cos w1| j 
ESRD an  ) | zs 4 ede) | 


by J 
Wir wollen noch vereinfachen und nehmen ¢p<1 an, da in der Praxis aus Sicherheitsgriinden 
meist nur Langskrafte P)<P,, zugelassen werden. Wir kénnen dann k? — ep = k* fiir alle 
k>3 setzen. Da uns ferner nur Erregerfrequenzen w < 2w,, oder w <2 @,, interessieren, nehmen 
die Quotienten (c/w, ,)? naherungsweise hichstens die GréBe 

\ by _3 PF? Gaara 
(2) = oder 4-10 i i 1-4 FI 
an und kénnen fiir kleine n in den meisten praktisch interessierenden Fallen fiir alle k >3 gegen- 
tiber 1 vernachlassigt werden. [Vgl. hierzu die entsprechenden Abschatzungen unter Ziff. 4a)]. 
Mit diesen Naherungsannahmen erhalten wir nach einigen Zwischenrechnungen aus (15) 


— B,v"’ = Myy+ (Moy — My) a i ae x) 


_ ep) 2 1\’ . I 9 
lane — |2 Po = + Myo + Moo} sin j 


males 


1 
4—ep 


(Moy — M,,) sin oe! aye 
+ cos on M+ (Mz, — My) 7 +7" — x) + 


I fh > ay amy 71 EP) + Mat Ma) sin 


(16) 
eT 


@ l 


(x) |a—er) 


W2¥ 


ie 


| (M,,— M;;) sin — ; 


5, Niherungsweise Zuriickfiihrung des allgemeinen Belastungsfalles auf die Belastung durch 
gleichgroBe pulsierende Endmomente, Es wire grundsitzlich miéglich, die Lésung (16) der 
Gleichungen des Grundbewegungszustandes im allgemeinen Belastungsfalle in die Gleichungen (7c) 
und (7d) einzufiihren und diese zu integrieren. Versucht man das, so iiberzeugt man sich jedoch 
bald davon, daB der damit verbundene Rechenaufwand in keinem Verhaltnis zum erzielbaren 
Erfolg stehen wird. 

In einem speziellen Falle vereinfacht sich jedoch alles erheblich: Wenn namlich nur zwei 
Endmomente M,-|-M,cos wt von gleicher GréBe wirken und au®erdem (w/a, )* gegeniiber 


eins vernachlassigt werden kann, ist einfach B, v’’ = — (M)+ M, cos wt) und die Differential- 
gleichungen des Nachbarbewegungszustandes lauten 
w ae wv |. ; (M,+ M, cos wt) #’ = 0, 
b+ ME or 2 yy eae = 
40 Us et, Que 
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Auf diesen Fall miiBte der allgemeine Belastungsfall reduziert werden, um eine méglichst grofe 
Vereinfachung zu erzielen. Das kann z.B. auf einem der folgenden beiden Wege geschehen: 


a) Die Durchbiegung der Grundbewegung kann in der Form 
v (x,t) = v,(x) + vy(x) coswt 
geschrieben werden, wo v,(x) die statische Durchsenkung und v,(x) die Amplitude der durch die 
periodisch veranderlichen Krafte hervorgerufenen (dynamischen) Durchsenkung bedeutet. Nun 
wird eine méglichst einfache Durchbiegungskurve 
v(x,t) = v,(x) cos wt 

angenommen mit der Forderung, daB die Maximalwerte ihres statischen und dynamischen 
Anteils gleich den entsprechenden Anteilen von v(x,t) sind. Wenn dann auch der Verlauf der 
Durchbiegung in Abhangigkeit von x in beiden Fallen verschieden ist, so wird man auf diese 
Weise doch im allgemeinen eine verniinftige Naherung erwarten diirfen. Fiir (x,t) nehmen wir 
die bei gleichen Endmomenten M, + M,coswt unter der Voraussetzung (w/a, »)? <1 ent- 
stehende Durchbiegungskurve 


(x,t) = = x) (M,-+ M, cos w1) 


mit 
= = sail! Bisa = Efe [2 ee 
Times = 9,(+, 1] = ge My, lmao = ta(-4 1) = 2 Mh. 
Die Bedingungen 2%, max==Vsmax UNG Udmax = Vamax liefern uns dann die konstanten Naherungs- 
momente M,und M,. 
b) Die durch zwei verschiedene Belastungen hervorgerufenen Grundzustande werden sich 


nicht wesentlich voneinander unterscheiden, wenn in beiden Fallen die im Trager aufgespeicherte 
Formanderungsarbeit dieselbe GréBe hat. Mit Gleichung (16), fiir die wir abgekiirzt 


B,v''(x,t) = M,(x) + M(x) coswt 


schreiben, hat die Formanderungsarbeit, welche durch eine allgemeine Belastung zu irgend einer 
Zeit tim Trager aufgespeichert wird, die GréBe 


1 I 
es | [v"(a,)P da = i [M,(x) + M,(x) cos wt} dx. 
x=0 x=0 


Andererseits ist die durch gleiche Endmomente M = M,-+ M, cos wt erzeugte Formanderungs- 
arbeit 
1 
1 a : 
oR | [M,-+ M, cos wt]? dx. 
=0 


— 


UH= 


Durch Gleichsetzen beider Arbeitsbetrage ergeben sich Bestimmungsgleichungen fiir die Nahe- 


rungsmomente. 

Der Nachbarbewegungszustand kann also naiherungsweise stets durch die Differentialglei- 
chungen (17) beschrieben werden, wenn dort die Momente M, bzw. M, durch die aus einer der 
beiden Naherungen folgenden Momente M, bzw. M, ersetzt werden. 

6. Der Nachbarbewegungszustand bei Belastung durch gleichgro8e pulsierende Endmomente, 
a) Lésungsansatz. Nach Abschnitt 5 geniigt fiir eine naherungsweise Behandlung aller Be- 
lastungsfiille die weitere Untersuchung der beiden gekoppelten Differentialgleichungen (17), d. h. 
der Gleichungen fiir Belastung durch gleich groBe pulsierende Endmomente. Wir transformieren 
sie in eine dimensionslose Form, indem wir die neuen Variablen 


Eyal ae a) Bee 18 
cath ened rhe >) 
und die dimensionslosen Tragerkonstanten 


= G Labs — Mir 1 
ii gw? Ho = By h’ [My = 4%} fas b2— Byh (19) 
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mit M,,—=erstes kritisches Moment des statischen Kipproblems einfiihren. Wir wahlen als 
Parameter 
M M. 
6) ee (20) 
Mz, Mkr 
und setzen voraus, daB beide klein gegeniiber ,,eins“ sein sollen; damit haben wir alle praktisch 
wichtigen Falle erfaBt, bei denen die Momentenamplituden weit unterhalb der kritischen Kipp- 
momente liegen. Die Stérfrequenz w geht bei der Transformation in die dimensionslose » otér- 
frequenz*‘‘ 


/ 
= | o 2x, 0 (21) 


iiber. Diese Umrechnung gilt fiir alle Frequenzen. Wir wollen M,,, das sich aus dem ersten 
Eigenwert der Gleichungen (17) fiir den statischen Sonderfall «= == 0 ergibt, fir 
unsere weiteren Rechnungen als bekannt voraussetzen. Die Gleichungen (17) gehen dann itiber 
in die Gleichungen 


n+ Hi NY + py(E + & cos at) = 0, | 
- 22 
B44 ow — 120" + (e+ e008 07)" =0, | ee 
zu denen nach (8) die folgenden Randbedingungen hinzutreten: 
(0,2) = (1,7) = 4'(0,t) — C,7'"(0,2) = (1,7) + Ca7"(,1) = 0, | (23) 


(0,7) = (1,1) = 90,7) — C, 80,2) = # (1,7) + C,0'(1,7) = 0. | 


Die Integration dieser Differentialgleichungen soll mit Hilfe eines ,,St6rungsansatzes‘* in der 
Weise durchgefiihrt werden, daB dem geschlossen lésbaren ,,ungestérten* Randwertproblem 
(c,=0) das ,,gestérte‘* Randwertproblem (¢,-+-0) mit denselben Randbedingungen superponiert 
wird, wobei fiir die Lésung des letzteren eine Entwicklung von 7 und @ nach Potenzen von ¢é, 
angesetzt wird: 


nt) =nolEt) + Solr), PEt) =HlEN+ SAlErer. 24) 


Wir setzen dabei die Entwickelbarkeit von 7 und # nach Potenzen von ¢, voraus und wollen 
annehmen, daf} fiir diese Entwicklung ein von Null verschiedener Konvergenzradius besteht ; 
wahrend die Konvergenz dieses Verfahrens fiir gewéhnliche Differentialgleichungen unter weit- 
reichenden Voraussetzungen gesichert ist, ist diese Frage allerdings fiir partielle Differential- 
gleichungen, soweit bekannt, noch ungeklart. Man wird aber sicherlich Konvergenz erwarten 
diirfen, da die Lésung fiir ¢,= 0 existiert, und es physikalisch plausibel ist, daB sie hinsichtlich ¢ 
stetig ist. Mit diesem Ansatz erhalten wir zuniachst die beiden Gleichungen: 


Ho + xin’ = — oe 8» 


9 1 giv 2 97 Za (29) 
Osim tee — #3 9) = — Me &oMo | 
fiir das ungestérte Problem; d.h. also die Differentialgleichungen der freien Transversal- und 
Torsionsschwingungen. Weiterhin ergeben sich die 2 n Differentialgleichungen 

Nn + 2x2 ntY = — w, cosat: ®_1, | 


n 


(alo ae (26) 


+. 1 fp a 
a, silts a any. paid x3 oe = — Ue COS AT °* Yn—1 p) | 


in deren rechte Seiten die Lésungen der jeweils vorhergehenden Gleichungen eingehen. Aus- 
gehend von den beiden Gleichungen fiir n=1, auf deren rechten Seiten die Lésungen von (25) 
einzusetzen sind, hat man die 2n Gleichungen sukzessiv zu lésen, wobei jede der Funktionen Nn 
und #,, die Randbedingungen zu erfiillen hat. Diese Gleichungen lassen sich mathematisch in 
anderer Weise behandeln als die ahnlich aufgebauten, jedoch von zwei Ortsvariablen abhangigen 
Bewegungsgleichungen der Arbeit von W. Kucharski!: In der II. Mitteilung jener Arbeit muBten 
die rechten Seiten der dortigen Gleichungen (5), die die Randbedingungen nicht erfiillen, nach 
den Kigenfunktionen der freien Schwingungen entwickelt werden. Das ist beim vorliegenden 


1 Vel. FuBnote 3 auf S. 137. 
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Problem nicht erforderlich; die dabei auftretenden Schwierigkeiten lassen sich hier durch Pro- 
duktansatze fiir die 7, und %, umgehen, welche die partiellen Differentialgleichungen (26) in 
Randwertaufgaben fiir gewohnliche inhomogene Differentialgleichungen mit konstanten Koeffi- 
zienten tiberfiithren. [Naheres dariiber unter c).] 


b) Freie Schwingungen. Wir wollen dabei zwei Falle unterscheiden: Fall 1: Es sei 
keine statische Vorbelastung durch ein Moment vorhanden (¢,—0). Man kommt dann 
mit den Ansatzen 


Mo(&>T) a Uor(E)pox(t) und (Et) = YS? Vou(E) pou(t) (27) 
= k=1 
zu einer Lésung: Aus (25) folgen die gewéhnlichen Differentialgleichungen 
Por + We Por = 0 und Wok + Wak Wok = 0 (28) 
fiir die Zeitfunktionen q(t) und yw ;(7) sowie 


2 
L 


w? 
Use—= Un=0 und Vo, — 428 Vi, — 40 Vor = 0 (29) 


fiir die Koordinaten-Funktionen U);,(&) und V,,(&). Die Gleichungen (28) haben die k Lésungen 
Pox(t) = Aj, cos WT + B, sin W447 ; Wor (T) = C, cos Wa, t + Dy sin Wg, T- (30) 

Die k Lésungen der Eigenwertdifferentialgleichungen (29) sind 
Uo (&) = a1, CofA, E+ ay, Sin dA, E+ ay, cos A, & + ay, sin A, , | 
Von(€) = b,;, Gof mF + by, Gin», & + b5, cos %,E + by, sin %E, | 


wobei die k-ten (dimensionslosen) ,.Kigenfrequenzen“ w,;, bzw. Wg, der freien Transversal- bzw. 
Torsionsschwingungen mit den Eigenwerten A,, »,, », durch die Gleichungen 


(31) 


vt a 
2 2 k k =u = 2 
One — 4, 4, OS eer H2V,, Ve, — 4x2 (32) 


verkniipft und mittels (21) auf die wirklichen Eigenfrequenzen w,;, bzw. w,, umzurechnen sind. 
Die Forderung, daB die Uj; bzw. V; die Randbedingungen erfiillen miissen, fithrt zu acht 
homogenen Gleichungen fiir die Koeffizienten a,, und 6b,, von (31) von der Form 


4 4 
Dra; 41, = 9; > pir 6 = 0 (== 2,3, 45 
r=l1 T= 


welche das Verschwinden der Koeffizienten-Determinanten A(a;,) und A(f;,) zur Folge haben. 
Die Ausrechnung ergibt die beiden transzendenten Frequenzgleichungen zur Berechnung der 
Eigenwerte A, und 1, bzw. 9, 


G (Ax) — 2 Ay (Cy + Ce) B (Aa) + A, G Co S (An) = 2 
2 v, cos V, Gof v% + (C, + Cy) (vg + v2) (cos v, Sin», — sin ¥, Cof v,) — | 
— = sin »;, Sin », [C, C, (vg + %)? + ve — vg] = 2 % | 
ke 

mit 1 . 

% (Ax) = Coj An sin An > Sin Ar cos Ax | 

© (A;) = )] Ak cos At > S(A,) =e Sin Ax sin At : | 
Fir spezielle Einspannungsfalle vereinfachen sich die Frequenzgleichungen: Wir erhalten z. B. 
fiir freie Auflagerung an den Trigerenden (1/C, = 1/C, = 0) erstens sind, = 0, d.h. A, ke 
und daraus nach (32) wj, = (ka)*x?, sowie zweitens sin CGinv,— 0, d.h. v,=—k2x, woraus 

oh = (ka) + 8 (kx)? folgt. 


Die Gleichungen (31) lassen sich nach Einsetzen in die Randbedingungen auch in der Form 


schreiben 


(33’) 


Uo, (£) = M, (CofA, & — cosA, E+ P, Sina, & + Q, sin A, €), 
Vou () = N, (Gof 4 — cos E+ R, Gin E+ S, sin 7,4) | 


1 Die Funktionen B(Jx) usw. sind bei Hohenemser- Prager, Dynamik der Stabwerke, tabuliert. 


(34) 
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mit den Konstanten 


M, = 4%, = beliebig, N, = 0, = beliebig, | 

Pee 1 a sin A, 

0, {  —sin Ak — Git Ax (Got Ae COE ata (ging : ’ (34’) 
R Les 1 jhe _— 2 2) sin Vi, / 

s — Vy, sin Vik a Vi, Sin Vy, (er) (Go} » Mk hone 4) ae fi (% ait t) (Sin | | 


Die Lésungen der Differentialgleichungen (25) des ungestérten Problems lauten also fiir e¢,=0 

a a! [Uon(&) (A, cos Wo, t + B, sin w,;,7)], | 

i (35) 

Do (E,7) = 3 [Vou (€) (Cy, cos wax t + D, sin wa), )]. | 
k=1 


Die Eigenfunktionen U,;,(&) und V,;,(&) bilden dabei fiir alle Randbedingungen zwei orthogonale 
Funktionen-Systeme, was direkt aus den Differentialgleichungen (29) geschlossen werden kann: 
Sind /,, und J, zwei verschiedene Eigenwerte und U,,,, Uy, die zugehorigen Eigenfunktionen der 
ersten Gleichung, so ergibt sich durch Multiplikation von 


USK (8) — Am Uom(§)= 0 mit Upn (E) 


und 


UIV(E) — A, Usn(&)= 0 mit Urn (&) 


und Bildung der Differenz sowie durch zweimalige Teilintegration 
1 1 

(AS — AS) f Uon() Uo (6) dé = f (UY Uy, — UNV Uy,) dé 
£=0 £=0 


— [Uom Usn — Uon US m — Uo m Ugn+ Uon Usnla 2 


Durch Einsetzen der Randbedingungen wird die rechte Seite dieser Gleichung Null; daraus folgt 
wegen /,, +A, die NEE NENG 


[onl sae) de= 0 (m + n). (36) 


Entsprechend kann auch gezeigt Haren daB die V,;,(&) orthogonale Eigenfunktionen sind. 
Die beliebigen Konstanten M;, und N, werden durch die Normierung 


LUOE=1, JVI e@ae= (7) 


festgelegt. 


Fall 2: Statische Vorbelastung durch Moment M, (é)=+-0). Wir kommen jetzt mit 
dem Ansatz 


no(Est) = 2 [Usi(é) ge], (67) = oy [Vou (2) Pox(z)] (38) 


aus. Durch Einsetzen in (25) ergeben sich die Differentialgleichungen @,,- wi Pox mit den 
k Lésungen 


Pox(t) = A; cos w, T+ B, sin w;, T (39) 
fiir die Zeitfunktionen und die gekoppelten Gleichungen 


Une + Se Vi,= 0, | 


Vox — 422 Vou — 4.0 Vou + 4 Me & U0, = 0 


fiir die Koordinaten-Funktionen, aus denen man durch Elimination von V,, die k Differential- 
gleichungen achter Ordnung 


IV = 
Uo. — 


(40) 


VIII yea he IV x5 a Ok 
Uok — 4 20, Uo — sa (WET 4 ta Ha 8 1 4 xp 8) Vow + 455 & Uo. + 4 U0 
1 bil 
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fiir die V,;, bekommt, nach deren Integration sich die V); aus einer der beiden Gleichuugen (40) 
bestimmen lassen. Fiir den Spezialfall der freien Auflagerung an den Tragerenden erhilt man 
z. B. mit dem Ansatz U);(&) = M, sin ka € nach kurzer Rechnung 

Wk 1 1 

oy int obs) YA (obs — obs)? + ma macoli) 
Entsprechend (39) kénnen freie Schwingungen mit den verschiedenen Mp1» We gleichzeitig auf- 
treten, so daB die Gesamtlésung lautet 


No(,T) me [Uox(&) (Ay cos Om, t+ By, sin w, 17+ C, cos ayyt-+ Dy sin Wx 7)], 
=1 


Bo (Et) = [Von (§) (Ay cos 7+ By, sin oy, 7+ Cy, cos Wyo t+ Dj sin Wy, 7)]- 
k=1 


c) Erzwungene Schwingungen ohne statische Vorbelastung. Die Differential- 
gleichungen (26) sollen der besseren Ubersicht wegen fiir den Fall 1 (e,= 0) integriert werden. 
Der im folgenden dargestellte Rechnungsgang wird aber auch fiir ¢)-- 0 im wesentlichen der 
gleiche bleiben. 

Die in den Lésungen (35) des ungestérten Problems auftretenden Konstanten A,, B,, C, 
und D,, kénnten bestimmt werden, falls fiir t—0 bestimmte Bedingungen fiir die Anfangsaus- 
lenkung und die Anfangsgeschwindigkeit vorgegeben waren, z.B. durch 


fh = No(é, 0) = 5 4. Uox(€), £1(&) = no(é; 0) = oe B, Uy, (€), | 


zs (41) 
Fale) = Vol, = 3 Gy Vou(S), — BalE) = Hol, 0) = ¥ wan Di Vor(é). | 
=1 
Im vorliegenden Falle sind die Anfangsbedingungen und damit auch die Funktionen (41) will- 
kiirlich wahlbar. 
Setzen wir zur Abkiirzung 
Ay, cos (ox + 9%) T+ By sin (Wy, + 90) T= Yn (T)> | 
C;, cos (wa, + qa) t+ D, sin (wax + Ga) T= Yn (7); | 


worin q alle ganzen Zahlen einschlieBlich Null bedeuten kann, so erhalten wir mit (35) aus (26) 
(fiir n= 1) und wegen 


(42) 


1] i 
Yq,k COS AT = = (Pq—1,.e 1 Pq+1,k)> Pq, COS AT = = (Yq—1.e Wq+1,k) 


folgende Gleichungen zur Bestimmung von 7, und 0: 


iy + in = — DV (arpa Pere) Ie | 

ie (43) 

dy + zo Stipes =>, [Uo (Pa.e + %+1,4)] | 
k=1 


Die Lésungen der homogenen Gleichungen entsprechen genau den Lésungen des ungestérten 
Problems, jedoch mit anderen willkiirlichen Konstanten A,, ..., D,, die gleich Null gesetzt 
werden. Somit interessiert nur je eine partikulare Liésung der inhomogenen Gleichung (43), 
fiir die wir ; 
ny (&, t) = 3 [U_1,14(&) p_a,n (4) + U4 1,14 (4) 41,4(7)1) | 
ee (44) 
9, (Et) = (Varn (€) p—a,e (0) + V4r,14(€) 941,27) | 
k=1 


ansetzen. Der erste Index der Koordinaten-Funktionen bezieht sich auf die jeweils mit ihnen 
verkniipfte Zeitfunktion (also auf q), wahrend der zweite Index die zu integrierende Gleichung 
von (26) angibt. Mit diesen Ansatzen und mit 


Pq, (7) = — (Won + 9%)” qx (T) » aq, (T) = — (Wan + 1%)? Pa, 4 (T) 


150 Reckling: Die Stabilitat erzwungener harmonischer Schwingungen usw. Ingenieur-Archiv 


erhalten wir aus (43) 4 k Bedingungsgleichungen fiir die Koordinaten-Funktionen 


IV (4, F &)? 


Ux1,1% 2 Uxi,14 = Vip Oo) % F& + uz, Sit, § + | 
“2 . 
+ u3;, cos ¥,€ + U4, sin ¥, € = F,(é), (45) 
J 


VeLin— An2 Vierin— 4 (sn FO)? Ver i e= 044 C0) A, S + Vak Sin A, + 
-}- V3 k cos A, € + Vgk sin A; € = F,(&) > 


worin die auf der rechten Seite stehenden Koeffizienten u;;,, v;; bekannte, von den Randbedin- 
gungen und den Eigenwerten A, , und 7, abhangige Ausdriicke von endlicher GréBe sind. Es 
sind also 4k Randwertaufgaben zu lésen, bei denen alle Koordinaten-Funktionen die Rand- 
bedingungen zu erfiillen haben. Thre Liésungen setzen sich zusammen aus den mit den Kon- 
stanten c=;, bzw. d+;, behafteten endlichen Lésungen U+;,; bzw. Vz,; der homogenen Glei- 
chungen und den Partikularlésungen Uz;,, bzw. Vexp: 


t= 


4 = 
Uriin— > rin Usik UE are | 
: (46) 


4 — 
Veewiee lag epee | 


=o 
Die Konstanten c= ;, bzw. d= ;, sind so zu bestimmen, daB die jeweils vier Randbedingungen (23) 
erfiillt sind, fiir die wir 
Rey (Cz) = 0 bzw. S,t (V1, 14) — 05 (r=, «6 SA)S 
schreiben wollen. Setzen wir die Lésungen (46) hierin ein, so erhalten wir die inhomogenen 


Gleichungssysteme 


4 4 
D) Wit Ch SSR (Uer,) Dawe De dap re ele (47) 
ee Sil 


zur Bestimmung der c+ ;,, bzw. d= ;,, welche stets dann genau eine Lésung haben, wenn ihre 
Determinanten 


Ave) == 0, baw. Al ea5,) = Or (48) 
von Null verschieden sind. Dann errechnen sich die Konstanten aus 
= (Cae ah = [Deore sant 
Cz ik A (yrik) a bzw. d= 5, = A erin) 3 (ee Il, De oF 4) 5 (49) 


wo die Zahlerdeterminanten C+ ;;, und D- ;,im allgemeinen von Null verschieden sind. Falls die 
rechten Seiten von (47) Null werden, wenn also die Partikularlésungen der inhomogenen Glei- 
chungen fiir sich allein allen Randbedingungen geniigen, wie beim beiderseits aufgelagerten 
Trager, werden die Zahlerdeterminanten null, da dann alle Elemente einer ihrer Spalten ver- 
schwinden; dann verschwinden alle c+ ;, und d= ;, und die Partikularliésung ist bereits die voll- 
stindige Lésung. 

Verschwinden dagegen die Determinanten (48) des Gleichungssystems, so haben die inho- 
mogenen Randwertaufgaben im allgemeinen keine eigentlichen Lésungen. Sie sind aber in diesem 


Falle doch lésbar, falls jede Lésung U= ;,(&) bzw. Vx in(é) der homogenen Gleichungen zu den 
rechten Seiten F(€) bzw. F,(&) orthogonal ist, wenn also gilt! 


fUsnO@ROdé=0, baw. f[ Veu(OR(Qdé=—0, (i= 1,... 4). 


£=0 g=0 


Diese Bedingungen sind jedoch sicherlich nicht erfiillt, wie sich im einzelnen zeigen 1aBt. Wenn 
also die Determinanten (48) gegen Null gehen, haben die inhomogenen Differentialgleichungen 
keine eigentlichen Lésungen, so daf nach (49) die c= ;, bzw. dx ;, und damit auch die Lésungen (46) 
iiber alle Grenzen wachsen. 

Das Verschwinden der Determinanten (48) bedeutet gleichzeitig, da® die zugehérigen homo- 
genen Differentialgleichungen (45) Lésungen haben, die den Randbedingungen angepaBt sind. 
Diese sind aber wiederum bis auf die bei U+1,1;, bzw. V1, 1; stehenden Konstanten den Differen- 


‘ Vgl. z. B. Kamke, Differentialgleichungen, Liésungsmethoden, S. 188. Leipzig 1943. 
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tialgleichungen (29) gleich. Die Bedingung, da®B die Lisungen beider Differentialgleichungen 
denselben Randbedingungen geniigen miissen, ist aber nur dann zu erfiillen, wenn die Menge der 
Konstanten in den sich entsprechenden Gleichungen fiir beliebige k itibereinstimmt, d.h. wenn 


Win = (@az = a)?, bzw. w3,—= (ws, a)? , 


ist. Mit anderen Worten: Wenn die Erregerfrequenz die Werte 


a= |Wak 25 Wp n | (50) 
annimmt, wachsen die Lisungen (46) iiber alle Grenzen. Dariiber hinaus kénnen sie auch noch 
gegen unendlich gehen, falls das ihre Partikularlésungen 


aun x? nove ; 

U= ie = x2 vi — (4, = a) (Urk Co} 5 + Uz, St eae 

“? } — 5 tneaas. 
< xp! —(w4, Sa ax)? (Us i COS VES ae U4, SIN Vj, é) ° 

I 
V. = ) FI 
ERP Tf 48 —4 (0,4, =o)? (14 Co} Ay FP Mx Sit Ay €) + 
1 


: 
Ay, + 4 %3 22 — 4 (o,), 7 4)? (vax cos A, E + v4, sin Ay é) 


tun; d.h. wir erhalten auch noch Unendlichkeitsstellen fiir 


61 Oak ae ik als [Map = PE 241. 


eon Fa Y/4 aE x2 : (51) 


Zur Berechnung der zu ¢j proportionalen Stérung erhalten wir in entsprechender Fort- 
setzung des Verfahrens durch Einsetzen von (44) in die Gleichungen (26) fiir n= 2 


jp le — i DUNNE aie koe V 44k) Por Vein eGo ee V 21. te Os kl s | 


k=1 


(52) 


i i = ov oo — a> (Uf 1,14+ U4 1,14) Port U1 ep—2,n+ U1 149+ 2, 4] | 
k=1 

Bei der weiteren Rechnung ist zu beachten, daB in den auf der rechten Seite stehenden Zeit- 

funktionen ¢,(t) und w,;,(t) die Frequenzen w,;,, bzw. wy; auftreten. Da die homogenen Lisungen 

von (52) den Lésungen der ungestérten Gleichungen (25) entsprechen und in ihnen daher die- 

selben Frequenzen ,,, bzw. @;,, auftreten — was dem ,,Resonanzfall*: entspricht —, wird fiir 

die Lésungen 


75 — ee Pst. (Ax sin yh, T — B, cos Wy T) ais (Gh Vole a= U_»2%P—2,4 + U,. 2,2kP+ 2, hI | 
La 


00 ras rf 

a >) [Vou t (CG sin wa, t — D, cos Wax T) + Vor, Yor +- V—2,2%~—2,4-+ Ve 2,20 Y+ 2,4] | 
k=1 

angesetzt. Durch Einsetzen in (52) ergeben sich wegen (69) und (70) die Bedingungsgleichungen 


(53) 


i! o uv a ) 

UW = — 4 (Via t+ Vi i118) » | 
fisbetyrs ahi (54) 

& Vi, — 24Ves) = FZ Mak Me (U1,1% + (Ui 1,1k) | 


fiir die Koordinaten-Funktionen, wobei jetzt zur Abkiirzung nur dic fiir eine Naherung bis ¢; 


wesentlichen Gleichungen hingeschrieben wurden. 
Wir haben also ganz entsprechend wie beim ersten Rechnungsschritt zwei Randwertaufgaben 


zu lésen und erhalten als allgemeine Lésungen 
Ue — 64S ak E27 e5, € + Cah 4 Usk p ° } (55) 
Van = fiz Co} 2 22 + fox Sin 2 x, F 4 KIS. SULT R Aten 


in denen die Konstanten e;,, f;, wiederum durch EHinsetzen in die Randbedingungen ermittelt 
werden kénnen und die aus (46) entstandenen Partikularlésungen U),,, und V,,,, fiir alle Erreger- 
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frequenzen gegen unendlich gehen, fiir die die Funktionen U-,,;;, und V+i,1, das tun. Dariiber 
hinaus werden sie — wie sich zeigen 148t — auch noch unendlich fir 


O&O = Wok >» Wdk > | Max FH 4, x3]. (56) 


Die Randbedingungen machten die Auflésung von zwei inhomogenen Gleichungssystemen 


4 4 
Spe as bzw. DE idk = — Orel Var) (r= 2, 3; A), 


i=1 i=l 
erforderlich, welche dann eine Lésung haben, wenn die Determinanten 
A(l,ix) =14+ 4(¢,4 C,) + 12 ¢,C,, 
A (E, 54) = 4 Hq {1 x9 [(Cy+ Cy) (2 42 Co] 2 x — Gin 2 x) + Gin 2 xy (1+ 4 x2 C,C,)] — oj 2 x,} 


von Null verschieden sind. Beide Determinanten sind aber stets gréBer als Null, da die Kon- 
stanten C,, C, und x, positiv sind. Die Konstanten e;, bzw. f;, und damit die Lésungen der 
homogenen Gleichungen gehen also nur dann gegen unendlich, falls die R,4.(U2%p) baw. S,1.(Vexp) 
das tun; diese aus den Partikularlésungen U,;,, bzw. V2; entstandenen Ausdriicke wachsen 
aber fiir die durch (50), (51) und (56) angegebenen Frequenzen iiber alle Grenzen; damit sind 
simtliche Unendlichkeitsstellen der U,, und V,; angegeben. 

In dieser Weise setzt sich das Verfahren weiter fort. Fir die Erérterung des Grundsatzlichen 
soll daher hier die Rechnung abgebrochen werden. Die Lésungen sind dann 


ees Dy [(Uo. +e? Ux) Port & (U_1,1% P_1,4+ 9+ 1,14 Y41,—) + 
k=1 + ¢? U,, t (Ay sin 4, T— B;, cos w,,7)] + O-E?) , (57) 
O= by (Vex = ei Vor) Wok & (Vis OIL en Vi 11k P+ 1,k) =a | 
- +e Voy t (Cy sin ogy tT — Dy c08 4x t)] + O(€8) j 
Durch eine fiir das Weitere sehr wesentliche Umordnung der Reihenglieder werden darin 
alle bei den Konstanten A,, ... D, stehenden Ausdriicke zusammengefaBt : 
n= 3 tAvell (G0) + Beall (Got) + Goer) + Deel ED] + 08), | a 


P= S [Aw MEL (G2) + Be HE (Er) + Ge HE (E.7) + Di HE ED)1+ O14), | 
k=l 


wobei die so entstehenden Funktionen gl), ni (it, j= 1, 2) selbst wieder unendliche Reihen 


aus Produkten von Zeit- und Koordinaten-Funktionen sind, deren Ausrechnung wir oben bei 
n= 2 abgebrochen haben: 


ght (€,7) = [Ux (6) + 6? Us, (€)] COS Wyx~T + Ef Uz, (€)T sin @,7 + Ole) , | 
ghia (6,0) = [Uon (6) + 2 Uen (&)] sin oy, 7 — €2 Us, (&) 7 cos nt + O(c) , 

Bei (57) = ey [U_a,1x (8) 08 (ae— 4) t + U4.1, 1% (6) c08 (ana) 7] +0 (€) . 

gh} (6,7) = & [U_1, 14 (€) sin (a, — &) t+ U 41,18 (€) sin (@ax + 4) 7] + O(E}) . | 
Entsprechend lauten die ni) , in denen nur @,j;, mit @g;, vertauscht ist und die entsprechenden 
Koordinaten-Funktionen Vy; usw. an Stelle von U,;, usw. treten. 


(59) 


7. Die Untersuchung des Grundbewegungszustandes auf Stabilitét, a) Aufstellung der 
Stabilit atskriterien. Der Grundbewegungszustand ist stabil, falls die Differentialgleichungen 
des Nachbarzustandes keine mit der Zeit anwachsenden Lésungen haben. Die Lésungen (58) 
scheinen aber in jedem Fall — d.h. fiir alle Erregerfrequenzen — instabil zu sein, da beispiels- 
weise bei einigen der in ihren Kinzellésungen (59) auftretenden Glieder die Zeit t multiplikativ 
neben den trigonometrischen Funktionen steht und bei Weiterfithrung der Rechnung sogar 
Glieder mit 7, 7?,... vorkommen. Das ist aber offensichtlich physikalisch sinnlos. Wir diirfen 
daraus auch in dieser Weise noch nicht auf eine Instabilitat der Gesamtlésungen schlieBen, da 
in ihnen unendliche Reihen solcher Funktionen auftreten, iiber deren Verhalten noch nichts 
ausgesagt werden kann. 
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Zur Beantwortung der Stabilitatsfrage werden die in der Arbeit von W.Kucharski} abge- 
leiteten Satze zugrundegelegt und fiir den vorliegenden Fall zweier gekoppelter Differential- 
gleichungen erweitert: 

Da alle in (58) auftretenden Funktionen gl? und hi) linear voneinander unabhangig sind. 
k6nnen jeweils zwei von ihnen hinsichtlich ihrer Zeitabhangigkeit zu einem Fundamentalsystem 
fiir eine lineare homogene Differentialgleichung zweiter Ordnung hinsichtlich der Zeit zusammen- 
gefaBt werden, wobei z.B. fiir die Funktion gf) folgende drei Kombinationen méglich sind: 


Coy eT Pale gl), a main gO) 
oder gf} = mit gf), sf} mit gf?) (60) 
oder gf) mit gi?) , ol) eerie 


Jede der sechs gewohnlichen Differentialgleichungen, fiir welche diese Kombinationen ein Funda- 
mentalsystem bilden, ware in der weiter unten dargestellten Weise auf ihre Stabilitat zu unter- 
suchen. Dasselbe gilt auch fiir die ni?) . Die spater durchgefiihrte Berechnung der Substitutions- 
Koeffizienten fiithrt nun aber fiir samtliche Funktionen mit ungleichen Indices j auf Wider- 
spriiche, so das als zusammengehérige Grundfunktionen nur die mit gleichen Indices j be- 


hafteten — d.h. nur die jeweils ersten Kombinationen (60) — mit den zugehérigen Differential- 
gleichungen. 
B+ r)2g)=0 und A+ s)AW=—0 (j=1,2) (61) 


in Frage kommen. Die Lisungen (58) werden dann 
n= Si lel) + el - = 5 TAD + AQ. 
Jedes Paar der Funktionen gj) und Aj) ist jetzt fiir sich zu untersuchen. Bei Ersatz von t 
durch t+ p mit p = 5 entstehen die Integrale gj) = g()(r-+-p) und hid) = hi) (t+ p), die sich 
bei Wahl der richtigen Substitutions-Koeffizienten durch die linearen Substitutionen 
Ries nega eaaile 
ere |. Get (62) 
A = A) We + BY MQ | 
aus den Fundamentalsystemen g/j), hY) aufbauen lassen. 
Die Lésungen von (61) sind periodisch zweiter Art 
gl) = mii) gi) , hi) — nld) ni) (j= 1, 2), (63) 
wobei die Zahlen mJ), n¥ aus den quadratischen Gleichungen 
(mip)? — (a 4 + af») mY +1 = 0. 
(nf) — (BUR, + BYR, ,) + 10 
folgen. Diese vereinfachen sich noch weiter, wenn wir fordern, daB eine der Partikularlésungen gf) 


eine gerade, die andere eine ungerade Funktion von T ist: 


A)=K(—), = — Al), (64) 
und daB die Lisungen die Anfangsbedingungen 
sh (=e, (Ea). =O. AO =O. LO =e Ea) (65) 
erfiillen. Dann ist 
Lay ti (66) 


Beweis: Die Substitutions-Koeffizienten kénnen durch Partikularlésungen fiir p ausgedriickt 


werden: 
) i A) ai oh (py et) oo) Oh p\ ncaa ap 
Ce cad ? Bi (Pp) = cli) OD ? Bl) (p) = ef, OY, 8) (P) Sk 22, k 


1 Vel. FuBnote 3 auf S. 137. fi 
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Aus (62) gewinnen wir 


2) a | laid al) 
BO) Bun Bure  P) eyes ca 
od) ali) (j) (j) y : 
BE Xow] | Mat, w %2"2, bk 
wobei die im Nenner stehende Determinante gleich ,,eins‘¢ wird. Setzen wir hier t= — p, so 


wird mit (64) und (65) 
£0 (— p) = ef) = af, (0) = ef, af 


Andererseits ist aber 


womit die Behauptung (66) folgt. Die allgemeine Lésung 
oli) = gli j) oli) we» ai) (ok!) + ob «)p 43) 
By = By Pye rea (ee By) 8 (7) (67) 


ist stabil, sofern der Realteil des charakteristischen Exponenten [uy Null ist. Das Kriterium 


fiir die Stabilitat der k-ten Lisung gi) lautet dann 


jo.) <1 Cy. (68) 


Entsprechend folgt fiir die Stabilitat der k-ten Lisung ni 


|PA,| <1 (j= 1,2). (68’) 


Zur Bestimmung der Substitutions-Koeffizienten a) ,, BY) , werden die in (59) auftreten- 


den Koordinaten-Funktionen nach den Eigenfunktionen wy n (€), Von (&) der freien Schwingungen 
entwickelt, z. B. in der Form 


Ure= 3 0) Von, Uae = Sw) Uins Vor > oH) V,,, usw: (69) 
n=1 ea — 
; —— 
mit A= | [ U2x,Urnd& (n=k) ce le ad (n = k) i ee 
Den tol) 5 = a : bi ; : 
lo  % c co 8 «6 (n+ k) | f Use Vom & (n ) 


Mit diesen Entwicklungen kénnen wir dann die Funktionen (59) als Reihen aus Produkten von 
Koordinaten- und Zeitfunktionen schreiben: 


= DS Upn () GY (x) , h—= J) Von(§) H(t) (i, j= 1,2). (71) 
i = im — 1 

Bei Ersatz von t durch t-++ p entstehen nur neue Zeitfunktionen Ge. H,,;, wihrend die Koor- 

dinaten-Funktionen unverandert bleiben, d.h. dann kénnen die Zeitfunktionen gl) bzw. hii) 

mittels der Substitution (62) durch die Gleichungen 


J Uon Gi= aly, 3 Use C+ af), 3 Von OY. | 
a | (Aspe) 


3 Von HO= AO), 3 Von HO, +B, 3 Vou HY 
ah n=1 
mit den gt!) baw. hi/) verbunden werden. Da die Koordinaten-Funktionen der freien Schwingungen 


linear voneinander unabhangig sind, kénnen diese Gleichungen aber nur dann bestehen, wenn 
fiir jedes n 


B= of N+ OMG | as 
seta peel Tuam iii ae th 
gilt, wobei die Wahl von n freigestellt ist, da die Substitutions- Koeffizienten unabhangig von n 


sind. Wenn wir n —k wahlen, so erhalten wir die Funktionen Gy) und Ht }) durch Gleichsetzen 
der einzelnen Gleichungen von (59) und (71) sowie Vergleich der bei Uy, bzw. Vo, stehenden 
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Funktionen mit Beriicksichtigung der Entwicklungen (69). Die Substitutions-Koeffizienten er- 
geben sich, indem man zunachst in den aley und HY? t durch t-+ p ersetzt und so GY), HY) bildet, 
sodann die so gefundenen Ausdriicke auf den linken Seiten von (72) fiir n =k einsetzt und 
schlieBlich die bei gleichen Zeitfunktionen stehenden Funktionen gleich Null setzt. 


Die Ausrechnung ergibt fiir die Naherung bis zu Gliedern mit e} die Stabilitatskriterien 


| of) | =| cos on P + ef P a) sin ws, p| <1, | 2) ,|=|cosoup| <1, i. 
| PY. | =| eos ae P+ ef p bY) sin og, p| <1, |PP) .| =| cos. P| <1. | ) 


Fiir j= 1 kann aus (73) direkt der Verlauf der Grenzkurven zwischen stabilen und instabilen 
Lésungsgebieten in einer (g,«)-Ebene angegeben werden durch die beiden Gleichungen 


+ 1 — cos @,, p + 1—cosw,,p 
a=]/ ? Che _, (74) 


pas) sin Oy 1, P po sin 41, P 


wobei das bei ,,1‘* stehende Vorzeichen so zu wahlen ist, daB ¢, reell wird; auBerdem sind auch 
nur positive Werte der Wurzel physikalisch sinnvoll. 


Fiir j—=2 sind die Stabilitatskriterien fiir alle Frequenzen erfillt. Es kommt daher fiir die 
Stabilitatsuntersuchung auf das Gleiche heraus, wenn man von vornherein alle Lésungsanteile 
mit j—2 in den Lésungen (58) fortlaBt und entsprechend Gleichung (41) folgende Anfangs- 
bedingungen vorschreibt: 

a) Fir die Lésung 7 braucht nur eine Anfangsstérung hinsichtlich 7 zu erfolgen, d.h. es kann 
C, = D, = 0 gesetzt werden. 

b) Fiir die Lésung @ braucht nur eine Anfangsstérung hinsichtlich ? zu erfolgen, d.h. es 
kann A, = B,= 0 gesetzt werden. 


b) Ermittlung der kritischen Erregerfrequenzen. Wenn man nur die kritischen Fre- 
quenzen und nicht den Verlauf der Grenzkurven der Instabilitatsgebiete ermitteln will, braucht 
man die fiir die zweite Naherung bis zu Gliedern mit ¢? einschlieBlich maBgebenden Koeffi- 
zienten a) und b{*) nicht vollstandig auszurechnen, sondern kann die kritischen Frequenzen 
sofort angeben: 

Da die Instabilitatsgebiete in einer (¢,,a)-Ebene von den Punkten ¢,—0 ihren Ausgang 
nehmen, tritt Instabilitat nach (74) dann auf, wenn @,; p=ra bzw. wa, p=ra (r= 1, 2, 3...) 
ist, bzw. wenn mit p—=2z2/« die Erregerfrequenzen die kritischen Werte «= 2 @),/r, 2 @ax/T 
annehmen. 

Zusiatzlich zu diesen Instabilitatsgebieten ,,erster Art‘+ kénnen aber noch weitere Gebiete 
,zweiter Art auftreten, die ihren Ausgang von denjenigen Erregerfrequenzen nehmen, fiir 
die die Entwicklungs-Koeffizienten a(*) und b) bzw. nach (70) die Funktionen U,)(&, a) und 
V,,(€ «) gegen unendlich gehen. Diese kritischen Frequenzen sind aber bereits durch (50), (51) 
und (56) angegeben. Es 148t sich ferner zeigen, daB bei Ausfithrung der Quadraturen (70) keine 
zusatzlichen kritischen Frequenzen hinzutreten. 

Zusammenfassend sollen noch einmal samtliche in zweiter Naherung ermittelten kritischen 
Erregerfrequenzen « aufgeschrieben werden, von denen im allgemeinen Einspannungsfall Insta- 
bilitatsgebiete ausgehen: 

ia cack ie eg tae alaee i Gah 


r , Ts 


2, Art: Ki Day Oa mls [Oak = VEXy|, [Mae F ME, | 


AR oo (75) 
on =e ayer 25 Me 


; | yn 4%, %5| « | 


8, Anwendungsbeispiele mit zahlenmafiger Durchrechnung fiir konkrete Einzelfalle. a) Der 
durch gleiche Endmomente belastete frei aufliegende Trager. Die Rechnung soll auf 
Frequenzbereiche beschrankt bleiben, fiir die (w/a, ,)?<1 ist. Die Gleichungen (27) geben die 


Lisungen des ungestérten Problems an mit den normierten Figenfunktionen 


Upon (€) = Vou (é) = /2sinkaé. 
AEs 
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Die Eigenfrequenzen der freien Schwingungen sind nach (32) 
Ae eel see ere ee 
Won = (kt)? xy, dan = = V(b)? + 473 


Fiir die in (54) auftretenden Koordinaten-Funktionen erhalten wir 


Ursin = Usi1n sinkré, Veiik = VF1,14 8D ka, | 


(76) 
Uy, = Us, sin ka, Vi, = vy, sin kaé | 
mit den Koeffizienten 
fe (k m)* ‘ Ph ye : (hate a | 
i V2 5, —(ax F 2)” V2 4, —@on F %) (76’) 


ke 0)? ka)? 
ee a ae V44,1k)> vy = EE (24.4% 4 U+1,1k) . | 


Die einfache Gestalt der Eigenfunktionen soll dazu ausgenutzt werden, um bis zu Stérungs- 
gliedern mit ¢4 einschlieBlich zu rechnen. Man erhalt nach langerer Rechnung, wenn nur die 
fiir die Rechnung bis zu Gliedern mit ¢4 wesentlichen Anteile aufgeschrieben werden, 


4 = SY [sin kav & + uy, t (Ay sin @y,t — B, cos wo, 7)], | 
a (77) 


0,= » [sin ka é- 4,0 (C, sin wa, tT — D; cos Wax T)] 


{ 
J 
mit den aus (76’) rekurrent zu bestimmenden Koeffizienten 


1 { = ae 1k fae 
Ty = | ID 2 2 Hy Pe 5 : uz 
ap ian | V2 Ua (v2aan + Oha,18) + + Us 


Oy, —% Op, +% 


l 
i) 
\ 


} 
| (7) 
Ea el 5 My fg (ket)! ( 21,1 k He Ne ae 
1 oa | eae iat wtaas) + eee wouik — SHE) at) | 
Wir verfahren hier, wie am Schlu® von Abschnitt 7a) angedeutet wurde, indem wir 
a) fiir die Lisung 7 die Koeffizienten C,— D,= 0 
und b) fiir die Lisung @ die Koeffizienten A,— B,— 0 


setzen. (Man kann auch fiir die Naherung bis zu Gliedern mit ef einschlieBlich zeigen, da damit 
alle Instabilitatsgebiete erfaBt werden.) Dann wird im Fall a) 


n= 3 {sin kav é [2 gon + €% Ug (Ay sin 5,7 — B, cos @ 5,7) + | (78) 
+ ef ug, (Ax Sin Wy, T — By, cos wy, 7)]} + Oleh). J 
Ordnen wir diese Lésungen entsprechend (58) um, so wird 
el (€, t) = sin kw [7/2 cos ayn 7 + (€2 ugy + ef Ugy) T Sin Coy t| + O(é), | 
6 (é, x) = sin ko [2 sin wy. 7 — (€2 uae + ef ax) T cos yn 7] + O(8). | 


Eine Entwicklung der in den Lisungen auftretenden Funktionen nach Eigenfunktionen der 
freien Schwingungen eriibrigt sich hier, da die Randbedingungen bereits die Zusammenfassung 
der in den Lésungen allein auftretenden Funktionen sin kaé mit verschiedenen Koeffizienten 
gestatten. Nach Bildung von gi) kénnen wir daher jetzt direkt die Substitutionen (62) anwenden 


und erhalten so durch Nullsetzen der bei zweckmaBig gewahlten Zeitfunktionen stehenden Aus- 
driicke die folgenden Stabilitatskriterien: 


(79) 


2 
oe = |cos cnn p(1 — stp? rue & is Sin Wp% P (Urn + ef Ugx)| <1. (80) 


Im Fall b) ergibt sich entsprechend 


2 
Pi =) COS dT Pp (1 eT Ps 4 ap tis SiN Wak P (Voz + €} M4,)| SI. (80’) 
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Pir eine Niherung bis zu Gliedern mit €? ergeben sich die Grenzkurven zwischen 
stabilen und instabilen Lésungsgebieten aus 


i 1 2 5 7 = 1p) 
pal \ cos (,, Pp) V2 oa eae he 1 — cos w,;, p) V2 


P Ug), $10 Wy, Pp P V2), 810 M4, P i) 
Da der Fall der freien Auflagerung nach Abschnitt 6c) eine Sonderstellung einnimmt, kommen 
hier nicht alle der in (75) fiir den allgemeinen Einspannungsfall aufgeschriebenen kritischen 
Frequenzen vor; insbesondere fallen die Frequenzen |@4, += @5m| mit (m-=Ek) fort, wahrend die 
tbrigen in |); — wa;| ttbergehen: Als kritische Frequenz haben wir also nur 
2 Oy, 2041 
Oe = ? as, P= 23s subi) — |@Mox + wax]. (82) 


Die letzten Frequenzen folgen auch direkt aus den Unendlichkeitsstellen der Koeffizienten (76’) 

und (77’). Man erkennt, das bei der Naherung bis 4 keine neuen kritischen Frequenzen hinzus 

treten. Treiben wir die Naherung noch weiter, so kommen weitere Instabilitatsgebiete hinzu, 

die von den Frequenzen « = 2 lWsn + 
r 


F wax| (r= 2,3, ...) ausgehen. Diese werden jedoch 
sicherlich von untergeordneter Bedeutung sein. 


& = My [Myp 


(————-__ Naherung bis et einschl. ; 


——-—-— Niaherung bis st einschl.; Instabilitatsgebiete der Naiherung 


bis et schraffiert.) 


Tragerlinge l=3.30m. 


OY 


= 


1 Wl \\ | WW / | // 
ne eels WA) i 
4 se es 4) | 
ar VY \\ 4) if! 
y 44 i 
ibe 6 | I 
4 | | ly 
az YY 4 | 
i | ' 
| j ! | 
| / 
07 | 
0-064 
o 40 | 15 | 20 Go | 30 ad a 
NE fy ee a ee ee ee eee 
Oy 4 S097 5 6 wy 7 Loy 8 g 10 wgyr3y «17 10 2ayy 
] | l l co sek) 
700 150 200 250 JOO 
Abb. 5. 1. Beispiel: Frei aufliegender I-Triger mit Momenten M cos wt. Grenzkurven zwischen stabilen und instabilen 


Lésungsgebieten. 


Als Anwendungsbeispiel wahlen wir einen Trager, der bereits von verschiedenen Ver- 
fassern in ihren Beispielen benutzt wurde!. Der Tragerquerschnitt ist in Abb.5 dargestellt. Die 
Tragerlange betrage = 3,30 m. Die Querschnittswerte sind: 


Bi Ooo0 tm, By 99 tm t= 363 Cm. 2.06 tm. 
Ferner ist : 
Wh= 6,6, @ = 0,00563 tsec*m™, M,,—= 17tm. 
Damit werden die dimensionslosen Tragerkonstanten 
#,== 0,381, x= 1,349, U.= 6,90, b= 0,943 . 
Fiir die Frequenzen gilt die Umrechnung a = 0,04125 . Die Rechnung soll fiir k=1 durch- 


gefiihrt werden. Die dimensionslosen Eigenfrequenzen der freien Schwingungen sind: w,,=- 3,76 
und wg,— 6,5. Fiir den vorliegenden konkreten Fall ergeben die Abschatzungen von Abschnitt 4b), 
daB die Naherung M, ~ M, = const fiir den untersuchten Frequenzbereich ausreicht. 


1 Vel. FuBnote 2 auf fos diate 
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Das Ergebnis der Rechnung ist in Abb.5 dargestellt; die Erregerfrequenzen a sind dabei auf 
die erste Biegeeigenfrequenz @,, mittels 6 = a/@), bezogen. Fiir d> 3,46 gibt es kein Instabili- 
tatsgebiet mehr, wahrend sich fiir 0 <6<1 die Instabilitatsgebiete haufen; dabei werden diese 
Gebiete aber immer schmaler, so da&B Abb.5 die wesentlichen Instabilitatsgebiete erfaBt. Die 
ausgezogenen Kurven gelten fiir die Naherung bis zu Gliedern mit e4, die gestrichelten fiir die 
Naherung bis ¢2. Fiir kleinere ¢, (etwa ¢,<0,2) stimmen die Kurven fiir beide Naherungen gut 
iiberein. Bemerkenswert sind die Frequenzbereiche bei a=2a,, und a=20@4,: Wahrend die 
Naherung bis ¢? an diesen Stellen ausgepragte Instabilitatsgebiete ergibt, reduzieren sich diese 
bei der Naherung bis ¢4 auf ganz schmale Gebiete, die praktisch durch eine einzige Kurve dar- 
gestellt werden kénnen. 

b) Der durch Endmomente belastete, am linken Ende frei aufliegende, am 
rechten Ende vollstandig eingespannte flanschlose Trager. Jetzt wollen wir nur 
noch («/@3y)?<1 voraussetzen und die Naherungsmethoden von Abschnitt 5 anwenden, indem 
wir das gemaB (16) von x und der Erregerfrequenz  abhangige Schnittmoment M, = — B, v’’ 
in geeigneter Weise durch ein nur noch von w abhangiges Moment ersetzen. 

Im vorliegenden Falle haben wir als Grundschwingung 


ie) 


4 1? 1 . 0% UW og, Toe So 
v(%,t) = pos My cosmt : (@ } sin = spp, zs in 
oe poe 


SS i ay oe ) = wel! 
a [sh papas er 


und zur Abkiirzung 
1 


racragn 


einfithren, geht v(x,t) iiber in 


41? 3 
M, cos @ t|(B(o) — 1) sin + os (1 — »)|= va(x) cost. 


CED reer 


Ersetzen wir nun v(x,t) durch die Durchsenkung @ (x,t) = 04 (x) cos wt mit der Mafigabe, daB vq (x) 


durch ein konstantes Moment M, hervorgerufen und Vmax == Vd max = Va (1/2) sein soll, so er- 
halten wir das nur noch von der Erregerfrequenz @ abhaingige Naherungsmoment 


M,= M, c + 3 (Blo) — i) (83) 


Wenn wir die zweite Naherungsmethode anwenden, so ergibt sich als Naherungsmoment 


M,=M,)/1 + 5 (Bo) — 1. (83') 


Schreiben wir ohne Riicksicht darauf, welche dieser Naherungen wir anwenden wollen 


M, = C(o) M, > 
so kénnen wir alle fritheren Rechnungen beibehalten und miissen nur am Schluf der Rechnungen 
in den Gleichungen fiir die Stabilitats-Kriterium gema$ (20) ¢, durch & = e, C(w) ersetzen. 
Fiir den flanschlosen Trager vereinfacht sich die zweite Gleichung (22) durch Fortfall des 
Gliedes + GIN 


Die Randbedingungen lauten mit 1/C, = C, = 0 nach (23) (0) = (1) = 7'(0) = n'(1) = 0. 


Die normierten Eigenfunktionen der freien Biegeschwingungen sind 


ee 2 Ak (Sin Ak sin Ak E —sin Ak Sin Ar €) 
|/7x[sin? 2x (Gin 2 Ak — 2 An) — Gin? Ax (sin 27, —24n)] 
Die erste Frequenzgleichung (33) vereinfacht sich zu 6(A,) = 0 mit den Wurzeln A, = 3,9265, ... 
In der Differentialgleichung (29) fiir die Vj; fallt wegen der verkiirzten Ausgangsdifferential- 
gleichung — Fehlen des Gliedes 1/4 J1V — das erste Glied fort, so da® die normierten Eigen- 


funktionen der freien Torsionsschwingungen 


Vou (&) = 72 sin % & (85) 


(84) 
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lauten. Die dimensionslosen Eigenfrequenzen der freien Schwingungen sind 

Opn An ie Odk= MHz = kr xp. (86) 
Damit ist das »ungestérte‘ Problem gelist. Das in Abschnitt 6c) dargestellte Berechnungsver- 
fahren liefert uns dann die Koordinaten-Funktionen U,, (€,0) und V,,(&,a) sowie schlieBlich nach 


langerer Rechnung durch Ausfithrung der Quadraturen (70) die fiir eine Naherung bis zu Gliedern 
mit ¢? einschlieBlich interessierenden Koeffizienten 


alt) — Pu Kilt) [ots 24 (Ke(A,) + K,(Ax) 04 tg oy) + 
+ oferta (Kel) + Ky(a) 042 ct o42)] + 


1 1 i zat all 
K ; 
+ K,(Ax) pao ual ee) ac rete (% = ork Mh 


(87) 
bik) — Hy Ma Ph Oat | i 7j, S (04) 0» ae 
2k 4x2 |i? — (a, —%)? |B (0_,) Cee Ae 
1 1 | VE (Opn) Cpe ‘| 
xi Pi, — (Wg, +4)? |BlOry) (Fi—O4,) 2 I) 
mit 
eae 2 <2 ee 1 ) 
2 *y i rte He ‘ ee = Lines, ee) 


Dabei sind K,(/,) usw. endliche Ausdriicke, die ihrer Lange wegen nicht aufgeschrieben wurden, 
wahrend die Funktionen Q(0.,), ©S(e4x) usw. durch (33’) definiert sind. Die Grenzkurven 
zwischen stabilen und instabilen Lisungsgebieten kénnen damit aus (74) errechnet werden. 
Die frither allgemein ermittelten kritischen Erregerfre quenzen ,,zweiter Art‘‘— vel. Gl. (75) — 
lassen sich fiir das spezielle Beispiel auch direkt aus den Unendlichkeitsstellen der Koeffizienten (87) 
entnehmen: Die alt) gehen fiir ctg or,— 0, also fir o+,= +ma(m= 0,1, 2,...) gegen 


unendlich, was 4= |), — mxzx,| zur Folge hat. Wegen max_= Wam nach (86) treten also 
kritische Erregerfrequenzen a = |@s, —- @am| auf, wie sie auch friiher angegeben wurden. Die- 
selben Frequenzen ergeben sich auch fiir O(o+,) = 0, d.h. fiir oz, —=A,, aus den Unendlich- 


keitsstellen der Koeffizienten b$). Au®erdem haben die Koeffizienten noch weitere Unendlich- 
keitsstellen, die zu den entsprechenden iibrigen in (75) angegebenen kritischen Frequenzen 
fiihren; die letzten Frequenzen von (75) kommen hier iibrigens wegen der vereinfachten Aus- 
gangsdifferentialgleichung fiir # nicht vor. 

Es wird folgendes Anwendungsbeispiel gewahlt: Als freie Tragerlange wird 1—1,0m 
angenommen. Fiir einen Tragerquerschnitt von der Héhe = 10cm und der Breite = 1 cm 
gelten folgende Querschnittswerte : 


Be=-15tm?, By, 0115 tm, == 8.42 cm", C= 0,2535 tm’. 
Ferner ist 90 = 0,8- 10-8 tsec2m~? und M,,= 0,946tm. Damit werden die dimensionslosen 
Tragerkonstanten 
4, = 0,2901, x)= 12,04, wy=— 4,553, w,—94,1. 


Fir die Frequenzen gilt die Umrechnung « = 0,0196 m. Auch hier soll die Rechnung fiir k = 1 
durchgefiihrt werden. Die dimensionslosen Eigenfrequenzen der freien Schwingungen betragen 
mit A,= 3,9265 nach (86) a, = 4,48, wa, = 37,80. Die Abschatzungen von Abschnitt 4b) 
ergeben hier, da erst die Quotienten (w/w, ,)? gegeniiber ,,eins‘* vernachlassigt werden diirfen. 
Wir wollen uns der Naherungsgleichung (83’) bedienen, um den Einflu8 des von w abhangigen 
Biegemomentes zu beriicksichtigen. 

Nach Berechnung von a!!! und 6!) erhalten wir die in Abb.6 dargestellten Instabilitats- 
gebiete in der (¢, «)-Ebene. Es zeigt sich, daB die von den Frequenzen |@p, + %A,| und 
|wa, = %,92| ausgehenden Bereiche so schmal sind, dafs sie durch eine einzige zur €,-Achse 
parallele Gerade dargestellt werden kénnen und daher praktisch nicht von Bedeutung sind. 

Fiir die gleichen Tragerabmessungen wurde die Rechnung zum Vergleich auch fiir den frei 
aufgelagerten Fall durchgefiihrt. Die Ergebnisse sind in Abb.7 dargestellt, wo die beiden Auf- 
lagerungsfalle miteinander verglichen sind, indem die Grenzkurven der Instabilitatsgebiete fiir 
ein bestimmtes Gebiet jeweils iiber dem Verhaltnis der Erregerfrequenz o zu derjenigen Frequenz 
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aufgetragen sind, von welcher das betreffende Gebiet seinen Ausgang nimmt. Die Instabilitats- 
gebiete fiir beide Auflagerungsfalle stimmen dann recht gut miteinander iiberein. Fiir eine 
Abschatzung wird man wohl unbedenklich von dieser Tatsache Gebrauch machen kénnen: Man 


h-focm 


sels 
b=70m 
Tragerlange L=100cm. 


Resonang der Grundbewegung 


iy l oe 
40 | 45 50 
| a ee 
2 Ces 
$041 Oy 47-057 \%74-Og7| — Wgz HY; +Oqy OgjtOp7 HpAq-57 HA 7107 (ek) 
7500 7150 2000 2250 2500 


Abb. 6. 2. Beispiel: Links frei aufliegender, rechts eingespannter Trager, belastet durch Momente M, cos wt. 
(Instabilitatsgebiete schraffiert.) 


braucht dann nur die @,;, bzw. wg, fiir eine beliebige Einspannungsart zu berechnen und die 
Instabilitatsbereiche vom frei aufliegenden Trager zu iibernehmen. Allerdings gilt dieser Ver- 
gleich nur fiir diejenigen Instabilitatsgebiete, die ‘von den in (75) angegebenen kritischen Fre- 
quenzen ,,erster Art“ sowie 
von |@a~— @o),| ihren 
Ausgang nehmen; das sind | 
aber auch die wesent- | 
lichsten Bereiche. | 

c) Der nur durch 
Axialkrafte belastete 
freiaufliegendeTrager. 
Gleichung (12) ist jetzt die 
grundlegende Differential- 
gleichung, in der wir der 
Einfachheit halber 


P(t)==-Pycosiw t 
annehmen. Sie geht dann 
mit 


é = pee) 
Rainy bt) iSx ae 
aan os a Oe — 
aa) a7 Og7+ aD : 
$47 7+ 67 ss ple /By (88) 
Abb. 7. Vergleich zwischen einigen Instabilitatsgebieten des eingespannten und frei P Q ; 
aufliegenden Tragers vom 2. Beispiel. 2 
Links freie Auflagerung, rechts Einspannung, Py — Pl 
———— Beiderseits freie Auflagerung. O) = Prr k2 72 By 
iiber in 
1 tet ae e(kan)* cose tty BO (89) 


mit den Randbedingungen 
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Wir nehmen wieder einen Stérungsansatz 


(517) = mo(S.t) + & tn (E,t)e” » (90) 
mit dem sich die n Differentialgleichungen 
Mn + oY = — (ka)? cosat: ny 4 pearl aces) (91) 


ergeben, in denen wiederum die rechten Seiten als Liésungen der jeweils vorhergehenden Gleichung 
fir n—1 bekannt sind. (In der ersten — ungestérten — Gleichung fiir n= 0 ist die rechte Seite 
Null.) Wir wollen die Naherung bis zu Gliedern mit ¢? einschlieBlich fiihren. Die Rechnung 
geht ganz analog wie beim Kipproblem vor sich, und wir erhalten als Lésung des ungestérten 


Problems 
k=1 


Mo(E.t) = DS) Uon(€) Gox(t) = yo: [sin kat € (A; cos wm, 7 + By, sin w, T)] (92) 
k=1 


mit @, — (k)?— Higenfrequenz der Biegeschwingung. 


Upon (€) = 2 sin ka€& sind die normierten Eigenfunktionen der Biegeschwingung. Durch 
Einsetzen in die zweite Gleichung (91) fiir n— 1 folgt 


+ 4 2° . 
ft mY = Fe D, (Mein kw (pant gas): 
il 
Der Liésungsansatz 


1 (U_11, Gaye t Usir G+1,8] 


fiihrt zu 


Uz 11k — U=+1,1% SID kaé 


mit 

OF O% 
—— SS u ; ——— — . 
a (2 (2a, —a)” aa «2 (20, +4) 


Aus der dritten Gleichung von (91) fiir n—2 folgt mit den fiir die Naherung bis e? wesentlichen 
Gliedern aus dem Lésungsansatz 


J 


fo ©) 


N= DS) (Us, t (A; sin w, tT — B, cos wo; 7)] 
ik 


fiir die Koordinaten-Funktionen 
of 


2V2a0(4 wok — a?) ; 


(93) 


U2, = Uz, sin ka & mit U,, = 
sodaB sich als Gesamtlésung bis zu Gliedern mit ¢? 


N= 21 {sin kaé [V2 vox + & (U—1,1k V1, + U4 1,1b P41,k) + EF Un T (Ay sin @, T— B, cos ox) |} 
we (94) 


ergibt. Ordnen wir diese wieder entsprechend (58) um, d.h. setzen wir 


| = 2 [Ag 81 + Be Se] » 


so kénnen wir die Stabilitatsuntersuchung genau so wie beim Kipproblem durchfithren: Wir 
bilden die Funktionen g,;, die sich durch eine Substitution entsprechend (62) mit geeigneten 0;;, 
aus den g,; aufbauen lassen. Durch Vergleich der Zeitfunktionen cos w,1 finden wir fiir die 
Substitutions-Koeffizienten 0,,,, und damit das Stabilitatskriterium 


aa el cos p+ ep tein a P ell. (95) 
Setzen wir die Erregerfrequenz « durch den Faktor 6 mit der Biegeeigenfrequenz w in die Beziehung 
ibt sich fi 93) uy —= ——“# dvcites witep =" = ana 

« = da,, so ergibt sich fiir u,;, aus (93) U2; Tees und weiter mit p—=— So 
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| ,4,~|— 1 fiir die Grenzkurven der Stabilitatsgebiete der Ausdruck 


(= 1 + cos Llp 6(62 — 4) 
e= (96) 
ect 
Se 


Die Instabilitatsgebiete in der (c,«)-Ebene nehmen jetzt ihren Ausgang von den Stellen O == 2/ih, 
d.h. von den kritischen Erregerfrequenzen « = 2 @,/n. Fir den beliebig eingespannten Trager 
treten bei der zweiten Niaherung bis ¢? noch die Frequenzen ,,zweiter Art |@, = @,,| hinzu, 


05 


G3 


S 
* 
—— oS 


G2 


‘ Jag 
0 
a OS | 10 15 20 
72 : : 
7 7 Zz 
UO G 6b G 7 é 


Abb. 8. 3. Beispiel: Durch P, cos mt belasteter frei aufliegender Trager. (Instabilitatsgebiete schraffiert.) 


wie sich zeigen la8t. Die Grenzkurven (96) haben fiir alle Tragerabmessungen und fiir alle k 
Giiltigkeit, wahrend die Grenzkurven des Kipproblems jeweils immer nur fiir einen ganz be- 
stimmten konkreten Fall gelten, da die dort auftretenden Entwicklungs-Koeffizienten von den 
Tragerabmessungen und von k abhangig sind, wahrend das bei 6 in (96) nicht der Fall ist. Hier 
brauchen wir daher im Gegensatz zu den beiden ersten Beispielen kein konkretes Zahlenbeispiel 
anzunehmen. Die aus (96) errechneten Grenzkurven sind in Abb.8 aufgetragen. Die Ergebnisse 
gelten nur fiir solche Erregerfrequenzen, die wesentlich kleiner sind als die erste Eigenfrequenz 
der Longitudinalschwingungen. Die Instabilitatsgebiete stimmen mit den von Mettler! mit 


Hilfe einer Mathieuschen Differentialgleichung berechneten Bereichen schon bei der Naherung 
bis €? genau iiberein. 


(Eingegangen am 2. November 1951.) 


Anschrift des Verfassers: Dr.-Ing. K. A. Reckling, Berlin-Charlottenburg, Technische Universitat, 
Lehrstuhl fiir Mechanik (Prof. Dr.-Ing. I. Szab.). 
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Uber die Einschaltung zusatzlicher Punkte beim Verfahren von Multhopp. 
Von J. Weissinger. 


Die Naherungsverfahren von Trefftz-Glauert } und Multhopp? zur Lisung der Prandtlschen 
Tragfliigelgleichung 


(8) = L(8) [a (8) — 7(8)] (1) 
‘mit der abkiirzenden Bezeichnung 
Fg) — 2 f) 
£9) ee 
6 


fitr eine beliebige Funktion f(#)) haben das eine gemeinsam, da® die gesuchte (dimensionslose) 
Zirkulationsverteilung y(#) als ein Sinus-Polynom P(#) m-ten Grades angesetzt und dieses 
durch strenge Lésung der m Gleichungen 


P(@,) = L(8,) [«(8,) — P(9,)] (y= 1,2,...,m) (2) 


oder kiirzer 
Po L, (a, — =) (yizas Ly 2056, 12) 


bestimmt wird, wobei die 9, m im Intervall (0,2) gelegene Punkte sind. Ein Unterschied besteht 
darin, da bei Trefftz-Glauert auf Giund des Ansatzes 


P(sy= S acsin wo 


n=! 


die Koeffizienten a,, bei Multhopp die Funktionswerte y,— y(0,) in dem Interpolationspolynom 


= wei Dy mH sin wd (3) 


als Unbekannte angesehen werden. Wahrend bei Treffiz-Glauert die Punkte @, beliebig gewahit 
werden kénnen, kommen die grofen rechnerischen Vorziige des Mulihoppschen Ansatzes nur 
dann zum Tragen, wenn die #,, das Intervall (0,7) aquidistant unterteilen: 


1,= (n= 1,2,..., m). 

Haben also die Vorgaben L(#), «(#) nur in einem kleinen Bereich des Fliigels einen ,,unregel- 
maBigen‘’ Verlauf, so kann man das bei Trefftz-Glauert durch Einschaltung einiger zusatzlicher 
Punkte im kritischen Bereich beriicksichtigen, wahrend bei der tiblichen Form des Multhopp- 
schen Verfahrens eine Hinzufiigung von Punkten nur in den Intervallmitten méglich ist und 
zwar gleichzeitig in allen Intervallen, so daf also der Wunsch nach Einschaltung eines Punktes 
mindestens eine Verdoppelung der Gesamtzahl der Punkte und damit der Unbekannten, also 
mindestens eine Vervierfachung der Rechenarbeit zur Folge hat. 

Zweck dieser Note ist der Hinweis, daB dieser einzige Nachteil des Multhoppschen Verfahrens 
insofern nur scheinbar ist, als man jedenfalls einen beliebig gelegenen Punkt hinzufiigen, den 
Grad des Interpolationspolynoms also um 1 vergréSern kann, ohne den Rechenaufwand merk- 
lich zu erhibhen. 

Sind namlich y,, yo, --+» Ym> Ys beliebige Zahlen, #, ein beliebiger, von &,—=na/(m- 1) 
(n= 1,..., m) verschiedener Punkt aus dem Intervall (0,2), P(#) der Ausdruck (3) und 
P,= P(#,), so gilt fiir das Sinus-Polynom (m- 1)-ten Grades 


les 
7 (9) =) | eG + 1)# 


~ sin (m —- 1) 


1 H.Glauert, Die Grundlagen der Tragfliigel- und Luftschraubentheorie. Berlin 1929. 
2 H. Multhopp, Luftf.-Forschg. 15 (1938), S. 153. 
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offenbar 
y (8s) == YS y (%,) 4h (v NB go och, m) 
und ; 3 ee 
= a ys — Ps m + 1)sin m + 1) 
CS ass sin (m + 1) Os 2sin B i 


Wie bei Multhopp erhalt man daher fiir den induzierten Anstellwinkel «;() = y (0) 


Xiy — y (3,) = b,,¥,— = by nn (v= 1,2,..., m) 


und iiberdies 
m+ 1 


a= KO) = P+ .— Peas ee Oe 


Man kann also die Werte y, (v= 1, 2,..., m) aus dem iiblichen Gleichungssystem von 
Multhopp 
Vy 7: L, (a, a Oi») (y = Ike Ze eo on m) 


unbeeinfluBt durch die Hinzunahme von *, bestimmen und dann y, aus der Gleichung 


m+1 ep 
Mes L, (a, ae Os) ae L, a, 2 (vs B,) Jsin Os P, D 


also 
= il , m +1 = 4 
Ys 1 m+ 1 Oe T 9 gin Bs ae P,) (4) 
L; | 2sin ®s 
ausrechnen. Eine einfache Umformung liefert 
pepe I S “eyelet mn hee a—Ss 
rag P,— P, (emir aa L 1—yp)sin wd, sin pd, 25 hays 


mit 


dsnsin On | sin® ers (n grade), ‘ 
oe (m + 1) (cos 8s — cos On)?” ear. | cog heual es (n ungrade) . 
Setzt man noch 
“arma baz, thes (6) 


so kann man fiir (4) schreiben 


1 m 
Vs = Whe (a. + oD i 7.) ° (7) 
3 (eal 


und der induzierte Anstellwinkel an der Stelle @, ist 
ais b. 5s —  OsnYn + (8) 
i 


Ist insbesondere #, = (2k + 1)a/(2 m—+ 2) ein Punkt des (2 m- 1)-punktigen Verfahrens, 
so werden alle a,, = - und damit die b,,, baw. b,, gleich den von Multhopp tabulierten Koeffi- 


zienten b,,, baw. b,, (v= 2 k + 1) des (2 m-+ 1)-punktigen Verfahrens. Der nach (7) berechnete 
Wert y, ist dann also gerade der Wert, der sich beim ersten Iterationsschritt des (2 m- 1)- 
punktigen Verfahrens ergibt, wenn man nach der Multhoppschen Rechenvorschrift als nullte 
Naherung die Lésung des m-punktigen Verfahrens beniitzt. 

Nach Multhopps Vorschlag wird eine beliebige Zirkulationsverteilung zweckmafig als Summe 
einer symmetrischen und einer antisymmetrischen Verteilung berechnet. Zu den entsprechenden 
Gleichungen von Multhopp tritt dann statt (7) die Gleichung 


m+1 
1 2 
ik Fasee By = ep ais Oia ae B Sh m+ 
n= aaa 


s s 
2 
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im symmetrischen und 


m—1l 

i oy 
een. Ie a, + > Con Yn r) C= by) Dain 

: n—1 


im antisymmetrischen Fall. 


Zusammenfassung, Es empfiehlt sich, den Wert y, der Zirkulationsverteilung an einer 
Zwischenstelle #, nicht durch Einsetzen von #, in das Multhoppsche Interpolationspolynom (3), 
sondern nach (7) auszurechnen, da auf diese Weise die Prandilsche Tragfliigelgleichung (durch 


Erhéhung des Grades des Naherungspolynoms um 1) nicht nur an m, sondern an m~-- | Stellen 
exakt befriedigt wird. 


(Eingegangen am 3. November 1951.) 


Anschrift des Verfassers: Prof. Dr. J. Weissinger, Hamburg 26, Wichernsgarten 1. 
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Die Auftriebsverteilung von Tragfliigeln mit Tiefensprung. 
Von J. Weissinger. 


1, Einleitung. Da bei unstetiger Anstellwinkelverteilung die tiblichen Methoden zur Lésung 
der Prandtlschen Tragfligelgleichung zu langsam konvergieren, hat Multhopp' ein Verfahren 
zur ,,Abspaltung“‘ eines Sprunges angegeben, bei dem die gesuchte Zirkulationsverteilung y dar- 
gestellt wird als Summe einer bekannten, nur von GréSe und Lage des Sprunges abhangigen 
Funktion yr und einer unbekannten Verteilung y*, die einer Prandtlschen Gleichung mit 
stetiger Anstellwinkelverteilung geniigt, also wie tiblich berechnet werden kann. Der Versuch 
einer schematischen Ubertragung dieser Ahspaltungsmethode auf den Fall eines Sprunges in der 
Tiefenverteilung scheitert daran, daB die Tiefenverteilung in der Prandtlschen Gleichung hin- 
sichtlich y nicht additiv wie die Anstellwinkelverteilung, sondern multiplikativ auftritt, so daB 
die abzuspaltende Funktion yz notwendig vom Werte der unbekannten Funktion y an der 
Sprungstelle abhangt. Im Folgenden soll gezeigt werden, daB die Ubertragung des Multhopp- 
schen Verfahrens, das insbesondere in seinen Bezeichnungen als bekannt vorausgesetzt wird, 
dennoch méglich ist und zu einer einfachen Rechenmethode — auch im Falle eines gleichzeitigen 
Sprunges von Anstellwinkel- und Tiefenverteilung — fiihrt ?. 


2, Das Rechenvyerfahren, Die mit dem Produkt hV der Spannweite b und der Anblasgeschwin- 
digkeit V dimensionslos gemachte Zirkulationsverteilung geniigt der Integrodifferentialgleichung 


y(n) = L(m) [o(n) — a(n), kurz: y = L(a — a,) (1) 


mit 


1 
__ a x(m) U7) me dy dn’ 
L(y) = Sey ee ai(7) oe | dn’ y—1"’ 

a 
in der »= 2 y/h die dimensionslose Spannweitenkoordinate, c’,.(7) die Profilkonstante, 1(7) 

die Tiefen- und a (7) die Anstellwinkelverteilung bedeuten. 
Die Funktionen L(y), «(y) sollen stetig sein aufer evtl. an der Stelle 7 = 7,, wo sie einen 

Sprung haben diirfen, dessen GréBe durch 


D LGje 20) LG 0) 5. fe. 
Gere 0s o = a(n, + 0) — a(n, — 0) (2) 


gekennzeichnet sei. Da y(n) iiberall stetig sein muB, gilt an der Sprungstelle die Gleichung 


Ys = (ns) = L(ys-r 9) [a(ns + 0) — a(n + 0)] = L(y. — 0) [a(ns — 0) — ai (ns — 9)], 
aus der man fiir den Sprung des induzierten Anstellwinkels 


ai(ns + 0) — a(y,s— 0) =o+T7y, (3) 


errechnet. 
Wahlt man in der Zerlegung 
(nq) = yun) + ¥*(n) (4) 
die Funktion y; so, da8 der zugehérige induzierte Anstellwinkel a, an der Stelle 7, den Sprung 
o- Ty, macht, so ist der von y* induzierte Anstellwinkel a} stetig, und y* kann daher mit ge- 
niigender Genauigkeit nach dem einfachen Multhoppschen Verfahren (ohne Abspaltung) be- 


rechnet werden. 


1 H. Multhopp, Luftf. Forsch. 15 (1938), S. 153. 

* VeranlaBt zur Veréffentlichung dieses Rechenverfahrens, das im letzten Kriegsjahr wahrend meiner 
Zugehorigkeit zur DVL-Adlershof entstanden ist, wurde ich durch eine Arbeit von F. B. Hildebrand (N. A. 
C, A. Techn, Notes No, 925), die mir kiirzlich zuganglich wurde. Hier wird das gleiche Problem nach einer 
Methode der kleinsten Quadrate behandelt, die in ihrer numerischen Durchfiihrung aber einen merklich 
eeeen Arbeitsaufwand erfordert als die im folgenden beschriebene Abwandlung des Multhoppschen 
Verfahrens. 
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Nach Multhopp hat nun die stetige Zirkulationsverteilung 


ae 2 
5 8 
io) — = (cos 8— cos @,) In my + #, sin (} = are cos n) 
$ ager 


die Eigenschaft, da8B ihr induzierter Anstellwinkel im Intervall — 1 < 1 <5 gleich 0, im Inter- 
vall ns< 74 <1 gleich 1 ist. Setzt man daher 


yi(q) = (0 + tys) f(m) » (5) 


so gilt, wie verlangt, 


: eo Rite a Sn <n, 

«lo tay, fir m <n St, (©) 
und y* geniigt der Gleichung 

yt=L-(a*—at), ao =a—ar—% . (2) 


Das zugehirige Multhoppsche Gleichungssystem lautet 
Appa bs (ay — b,,y7 + >’ 5,,y%,) (yee Canetee, m), (8) 
ial 


reicht aber, falls nicht gerade 7, mit einem der Multhoppschen Punkte 7, zusammenfallt, zur 
Bestimmung der m Unbekannten y¥ noch nicht aus, weil a* den unbekannten Wert y, enthalt. 
Man mu8 daher die Prandtlsche Gleichung auch noch an der Stelle 1, heranziehen und erhalt 
als zugehérige lineare Gleichung ? 


yt = Ljat— byt De! bn 78] (9) 
1 
oder unter Beachtung von (4), (5), (7) 
Y= Ya = De} oe Oils = bs Biel ie 2) Se) aie bn 7 > 


d.h. 
1 i 
ye| pe + Baa af] = a — sarah beef t 3 ban PE (10 
s n=1 


Fiir die unstetigen Funktionen L, a, «;,ist dabei der Grenzwert von links oder rechts zu nehmen, 
aber natiirlich fiir alle Funktionen von derselben Seite. Nimmt man z. B. den Grenzwert von 
links, so lautet (10) 


Tl ween b,(1—tf.)| = an 0) bis fo Den Yn - 


Ist tf, 41, so kann man statt y, leicht y¥ als Unbekannte einfithren. Denn nach (4) und (5) 
gilt an der Stelle 7, 


Y¥s= Vet (O+TYs) fs: 


folglich 
7 Ons ottys 
ye Ee Ce ce eee aE (11) 
Das vollstandige Gleichungssystem zur Bestimmung von y* und y lautet dann 
1 ia o a Tye fr << f * — 1 Y 12 
(b+ 7) B= a — ity L—tf ie ey OY es (v r) pee a5uMbe'S) ( ) 
mit 
0 fiir =S Il = Ny << "s ’ ( 
Oi Ty = %:1(N,) = Cae bys os Z| 13) 
=e CUE) oe le 


Fiir » = s mu8 wieder darauf geachtet werden, daB der Grenzwert fiir die unstetigen Funk- 
tionen L, a, «;1 von der gleichen Seite her genommen wird. 


1 Vel. J. Weissinger, Ing.-Arch. 20 (1952), S. 163. 


168 Weissinger: Die Auftriebsverteilung von Tragfliigeln mit Tiefensprung. Ingenieur-Archiv 


3. Bemerkungen. a) Sind mehrere Sprungstellen 7,4, ..-, sr vorhanden, so muff man yy 
als Summe von Ausdriicken der Form (5) wahlen, also ansetzen 


Y= P+ & (ot tse) Jo 
o= 
und zu den Multhoppschen Gleichungen noch die (9) entsprechenden Gleichungen an den Stellen 


so hinzunehmen. 

b) Normalerweise ist L(y) eine symmetrische Funktion. Dann ist y(y) zugleich mit a(m) 
symmetrisch bzw. antisymmetrisch, und man kann die von Multhopp fiir diese Falle angegebenen 
Vereinfachungen benutzen. Ist z. B. o(7) symmetrisch mit Sprungstellen bei 7 = 4, >0 und 


n= — 7, 80 nimmt man statt f(#) die Funktion 
oye atc eee O+ Ds 
9 sin —> cos —p 
F (9) = — | (cos & — cos #,) In | (cos #-+ cos B,) n —_,—-+ 29, sin 8 
I ae |\O9—@s| 0 — Os 
2 cos —5 
Dann bleiben alle obigen Formeln richtig, wenn man f durch F und b,, durch B,,, ersetzt und y, n 
eae! : 
ae von 1 bis 5 (m- 1) laufen 1aBt. 
L(0) 2 
40 
08 
06 
OY 
° M=15 
G2 a nach Hildebrand 
ee Li 4 
0 G2 O4 06 G8 10 0 Q2 OY G6 OF 40 
Abb. 1, Auftriebsverteilung des unverwundenen Abb. 2. Auftriebsverteilung des Fliigels II, a) mit konstantem 
liigels I. Anstellwinkel, b) mit Sprung des Anstellwinkels. 


c) Im Gegensatz zu den Gleichungen von Multhopp sind die Gleichungen (12) nicht immer 
durch Iteration lésbar, oder das Iterationsverfahren kann zu langsam konvergieren. Man iiber- 
zeuge sich also vorher, ob das Zeilensummenkriterium eine gute Konvergenz sicherstellt und 
wahle sonst ein anderes Auflésungsverfahren, etwa das von Gauf. 


4, Zahlenergebnisse!, Von Hildebrand? ist der Fall 


yi 9 OES | 
L(y) = Le I : 
a? aol! | 

also ein unverwundener rechteckiger Fliigel der Streckung (16/15) c’,. mit einem rechteckigen 
Ausschnitt in der Mitte behandelt worden. Das Ergebnis ist in Abb. 1 eingetragen, in der die 
ausgezogene Kurve die Auftriebsverteilung darstellt, die sich nach dem oben geschilderten Ver- 
fahren mit m= 15 (9 lineare Gleichungen) ergibt. Man sieht, daB beide Methoden praktisch das 
gleiche Ergebnis liefern; die Unterschiede iiberschreiten nirgends 1%. Zum Vergleich ist auch 
das Ergebnis der Rechnung mit m= 7 (5 lineare Gleichungen) eingezeichnet; auch diese Ge- 
nauigkeit diirfte praktisch schon ausreichen, zumal wenn man nach Liésung des Gleichungs- 
systems noch einige Punkte dazwischen schaltet . Kine iibertriebene Genauigkeit in der Lésung 
der Prandtlschen Gleichung hat praktisch schon deshalb wenig Sinn, weil in dieser Gleichung 


O(n) == 1 (di hig —=0)- (1) 


1 Fir die Durchfiihrung der Rechnungen habe ich Frl. Dr. H. Rohde zu danken. 
2 Vel. FuBnote 1 von S. 166. 
8 Vel. FuBnote 2 von S. 167. 
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Pfeilungs- und Randeinfliisse, die durch Spriinge in der Tiefenverteilung entstehen, nicht erfaBt 
sind. 


Als zweites Beispiel wurde der Fliigel 


\> Stim |) |P<e0, 6; 
i (11) 


L(n) = 
ees fir 0,6 <|y| <1 


betrachtet und zwar mit den Anstellwinkelverteilungen 
| (1 fir |n| < 0,6, 
oo) SUNS aia Reale yt ee Sein 

Auch hier stimmt das Ergebnis der Rechnung fiir m = 7 praktisch mit dem Resultat fiir m= 15 
itberein, wie aus Abb. 2 zu entnehmen ist. 

Sind ¢,, cy, zwei beliebige Konstante, so ergibt sich die zu 

c fiir <a 0500 
a(n) {2 Tl 
le, fir 0,6 <|y| <1 


gehérige Auftriehsverteilung y durch Superposition der beiden berechneten Verteilungen y, 
und ¥;, als 


Y= €,Ya-+ (C)-- ¢) Yo - 


5. Zusammenfassung, Das Verfahren von Multhopp zur Berechnung der Auftriebsverteilung 
wird auf Tragfliigel mit Spriingen in der Tiefenverteilung ausgedehnt. Der Rechenaufwand 
ist nicht viel gréBer als in dem von Multhopp behandelten Fall eines Sprunges in der Anstell- 
winkelverteilung. 


(Eingegangen am 3. November 1951.) 
Anschrift des Verfassers: Prof. Dr. J. Weissinger, Hamburg 26, Wichernsgarten 1. 
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Zur Stabilitaét erzwungener Schwingungen elastischer Kérper 
mit geschwindigkeitsproportionaler Dimpfung”. 


Von G. Grammel., 


1, Einleitung. In zwei grundlegenden Arbeiten hat E. Mettler! eine Theorie der Stabilitat 
erzwungener Schwingungen elastischer Kérper aufgestellt, dabei jedoch keine Dampfungs- 
einfliisse beriicksichtigt. F'. Weidenhammer ® hat fiir den besonderen Fall des axial pulsierend 
belasteten Stabes ebenfalls Stabilitatsuntersuchungen angestellt, wobei er die Mettlerschen Er- 
gebnisse auf einem zweiten Weg bestatigt und auch eine geschwindigkeitsproportionale Dampfung 
beriicksichtigt hat. Im folgenden wird nun gezeigt, wie man fiir jeden beliebigen elastischen 
Kérper bei harmonisch pulsierender Belastung das Stabilitaétsverhalten mit Beriicksichtigung 
einer geschwindigkeitsproportionalen Dampfung berechnen kann. Den Ausgang hierzu bilden 
diejenigen Bewegungsgleichungen des Kérpers, die man wie in M. I. ohne Beriicksichtigung der 
Dampfung erhalt. Diese bilden ein System von partiellen Differentialgleichungen, das nur Ab- 
leitungen gerader Ordnung enthalt. In diesen Differentialgleichungen werden die geschwindigkeits- 
proportionalen Daimpfungsglieder additiv hinzugefiigt. Jetzt hat man ein System von Diffe- 
rentialgleichungen, das die gedampfte Bewegung des wirklichen Kérpers beschreibt, aber auch 
Ableitungen ungerader Ordnung enthalt. Durch eine Transformation mit einem Exponential- 
faktor gewinnt man daraus ein neues Differentialgleichungssystem, das wieder nur Ableitungen 
gerader Ordnung enthalt. Es wird gezeigt, das dieses System die Bewegungsgleichungen eines 
Ersatzkérpers darstellt, der denselben geometrischen Bedingungen (Form und Lagerung) und 
denselben Belastungen unterliegt wie der wirkliche Kérper, aber an Stelle der geschwindigkeits- 
proportionalen Dampfung eine ausschlagproportionale Dampfung besitzt. Fiir diesen Ersatz- 
kérper entwickelt man dann streng nach M. II. die mathematischen Stabilitatsbedingungen, die 
durch Riicktransformation in diejenigen des wirklichen Kérpers verwandelt werden. Damit 
ist das Problem grundsatzlich gelést, und man mu nur noch aus den Stabilitatsbedingungen 
die zahlenmaBigen Werte berechnen. Fiir diese Berechnung wird ein Weg gezeigt, auf dem man 
in einem (@, €)-Koordinatensystem, wo w die Frequenz der Belastung und ¢ ein MaB fiir die 
Amplitude der Belastung ist, punktweise die Grenzkurven der gesuchten Instabilitatsbereiche 
des wirklichen Kérpers erhalt. 


2, Die Ausgangsgleichungen mit Beriicksichtigung der Dimpfung. Die Mettlersche Theorie 
geht aus von der Elastizitatstheorie endlicher Verschiebungen *. In ihr treten an Stelle der 
klassischen Dehnungen ¢ und y und der klassischen Spannungen o und t die Verzerrungsgriéfen 


3 
_ Oui, Ouj | > Up, OUp nee 
8ij Ox; | Oxi | } (i, = 1 2, 3) 
ies 


Ox, Oxj 


mit g;;—= gji, und die SpannungsgréBen k;;, welche mit den g;;im Bereich des Hookeschen Ge- 
setzes durch die Elastizitatsgleichungen 


kj 6 gi; (i + j) 
3 
= Spy ales > 
hom (8 De e| 


(tj == 3) 


(Sa 


mit dem Schubmodul G und der Poissonschen Konstante » — 1/m verbunden sind. Hierbei ist 
(%1, %,%3) ein kérperfestes Koordinatensystem, das im unbelasteten Zustand rechtwinkelig ist; 
u; (t= 1, 2, 3) sind die elastischen Verschiebungen eines Kérperpunktes in Richtung der un- 
verformten Koordinatenachsen. Vernachlassigt man die in den Verschiebungsableitungen 


* Dissertation Miinchen 1951. Referent: Prof. Dr. L. Féppl, Korreferent: Prof. Dr. R. Sauer. Der 
Verfasser ist Herrn Prof. Dr. L. Féppl fiir seine stets freundliche Hilfe zu herzlichem Dank verpflichtet. 

* E. Mettler, Ing.-Arch. 16 (1947), S. 135, weiterhin mit M. I. bezeichnet; Ing.-Arch, 17 (1949), S. 418 
weiterhin mit M. II. bezeichnet. ‘ i 

* F. Weidenhammer, Ing.-Arch. 19 (1951), S. 162. 

° R. Kappus, Z. angew. Math. Mech. 19 (1939), S. 271 u. 344. 
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quadratischen Glieder bei den Verzerrungsgrifen gij und den SpannungsgréBen k;;, so erhalt 
man mit den iiblichen Bezeichnungen der klassischen Theorie 

PEN ra 4. | 

ec bey ee ery 228). 

O;= k;;, tii kj; 


In der Elastizitatstheorie endlicher Verschiebungen hat man mit X, als Volumkraften und Pi 
als Oberflachenkraften in Richtung der unverformten Koordinatenachsen 1, 2, 3 als elastisches 
Potential (Formanderungsarbeit) des Kérpers 


ie : 
prea emetic : 
sak la oe ee 


als Potential der auReren Krafte 
3 3 
Voss fife Mimav— ff) 3 prado (2) 


und als kinetische Energie ohne Dampfung 
Pay ea 
eS ae ( ou dV, 3 
| | | é Di at 2 
wobei sich die Raumintegrale iiber das gesamte Volumen, das Oberflachenintegral iiber die ge- 
samte Oberflache des Kérpers erstrecken und dV ein Volumelement, dO ein Oberflachenelement 


des Kérpers ist. Nach dem Hamiltonschen Prinzip berechnet sich damit die eintretende Be- 
wegung des Kérperpunktes aus der Variationsgleichung 


ty 
6f/(T—A—V)d=0. 
ty 
Das Stabilitatsverhalten des Kérpers untersucht man nach dem Prinzip der kleinen Schwin- 
gungen, indem man dem bekannten Grundbewegungszustand u() (x,, %, %3, t) eine infinitesimale 
Nachbarbewegung u; (x, %5, x3, t) linear iiberlagert, so daB 
u; = ul) + 4; (i= 1 2.3) (4) 
gilt. Instabilitat tritt dann ein, wenn eines der u; mit wachsender Zeit t beliebig grof8 werden 
kann. Setzt man die u; aus (4) in die Ausdriicke (1), (2) und (3) ein, so erhalt man 
A= 4 2 4s AO, 
V=vO+4+VOQ+0, 
Oe eee ey 
wobei die oberen Zeiger (0), (1) und (2) jeweils Glieder nullten, ersten und zweiten Grades in 


den wu; kennzeichnen. 
Fiir die Grundbewegung gilt dann die Variationsgleichung 


ty 
6 { (TO — AO — VO) de— 0. (5) 
ty 
Ferner ist 
f(r AC) di -= 0, 
ty 


weil dieses Integral die erste Variation von oh A — V)=0 ist, wenn man fir die Ou; 
insbesondere die u; wahlt, und weil der Grundzustand ein méglicher Bewegungszustand ist. 
Die A, VO und T( sind vollkommen frei von allen u;, so daB also als Variationsgleichung 
fiir die wu; lediglich 
te 
6 | (T@ — A) d= 0 (6) 
ty 


iibrig bleibt. 
13* 
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Mettler zeigt, da unter den gemachten Voraussetzungen 


Al) — Abtass to TA es 


mit 
Autass = +] | 2 KY gl dV, (7) 
gl ue a ee 
JJ Jd om 
und 


sili? as 4 u 


ist. Man erkennt, daB Aj,),,, und 7T'(?) die klassisch berechnete Formanderungsarbeit und kine- 


tische Energie der Zusatzbewegungen wu; ist, wenn man sie am Kérper ohne die Grundverformungen | 


ul) anbringt. 

Bildet man von den Variationsgleichungen (5) und (6) die Euler-Lagrangeschen Differential- 
gleichungen, so weiB man aus der analytischen Mechanik', da diese Differentialgleichungen 
die Bewegungsgleichungen eines Kérperelementes darstellen. Kine geschwindigkeitsproportionale 
Dampfung erfaBt man also sinngema48 dadurch, da8 man bei diesen Differentialgleichungen 
ein Dampfungsglied 6 du(/dt bzw. B du,/dt (r= 1, 2, 3) additiv hinzufiigt. Mit diesem Damp- 
fungsglied erhalt man [vgl. M. I. Gl. (25) u. (26)] als Gleichgewichtsbedingungen, hier auf die 
Volumeinheit bezogen, aus der Variationsgleichung (6) die Differentialgleichungen 


3 
an) Se, ( ar, On 
ui (0) tries ur Our = 
= Ab | ee, Fe tae erage i ey SN 24, (a 


J 


mit den SS 


Say (1, 4; +> Pig = 0 (r= az, os (11) 


j= jas 
wobei nt der Oberflachennormalenvektor ist. Es gilt hier 


(1) — G o(} i ah 
kt Gg) (u J): 


Mes o(a = Eee | (12) 


oc 


dl wlan a 
wahrend die KK”) die klassischen Spannungen o und t des Grundzustandes sind. Entsprechende 
Gleichgewichtsbedingungen erhalt man aus der Variationsgleichung (5) fiir die Krafte der Grund- 
bewegungen ul) ohne die Zusatzbewegungen u,;. Diese Gleichgewichtsbedingungen mit Be- 
riicksichtigung der Dampfung sind nun keine Euler-Lagrangeschen Differentialgleichungen 
mehr, da sie jetzt auch Ableitungen ungerader Ordnung enthalten. Durch die Transformation 


cama 
U; (1, Xp, X%3, t) == @ 2° U; (x, Xo, 3, t) (u = 1, 2, 3) (14) 
geht (10) tiber in 
3 3 
ok) 3 
r rT] | : (0) OU, a 0?U, B? eo 
= Ox; | A, Ox; kj; Ox; O Ot2 40 U. Cc ile 2, 3), (15) 


und die Randbedingungen (11) in 


(1) 0U, 
| > k,j cos (in, x;) + > Pj ane (n=sle2a)e (16) 
—— j= 


a Vel. |e Joos, Lehrb. d. Theor, Physik. S. 98, Leipzig 1945. 
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Diese transformierten Differentialgleichungen (15) haben nur noch gerade Ableitungen der U; 
nach der Zeit, kénnen also wieder als Euler-Lagrangesche Differentialgleichungen eines neuen 
Variationsproblemes 


t 


6 { (T*— A*) dt:=0 (17) 


th 


ra f[fES eels [f[ider os 


7 
A * * 
A i Arias == Ma 


mit 


und, mit Beriicksichtigung der den Elastizitatsgleichungen (12) und (13) analogen Glei- 


chungen, 
3 
* G eri ats COP oUj 2 | 4a y dU; OU; 7 
3 3 
* 1 0) 0U, 0U 
= == ki? Bis epee! VE C 
A egy) ees J (5 ax, Ox; dt} (19b) 


der U; gedeutet werden. Es gilt jetzt statt (13) 


qin Casi CU; 
vj Ox; 0x; 


(i, 7= 1, 2, 3), 


wahrend die Elastizitatsgleichungen (12) formal identisch bleiben. Die Richtigkeit dieser Be- 
hauptung beweist man durch Ausfiihren der Variation (17). 

Wahrend (19a) und (19b) mit (7) und (8) formal identisch sind, hat sich der Ausdruck der 
kinetischen Energie (18) gegeniiber (9) durch Hinzutritt des zweiten Gliedes der rechten Seite 
geandert. 

Betrachtet man die U; als die Zusatzverschiebungen eines Ersatzkérpers, der denselben 
geometrischen Bedingungen (Form und Lagerung) wie der wirkliche Kérper unterworfen ist, 
und nimmt man bei diesem Ersatzkérper fiir die Zusatzverschiebungen eine ausschlagproportio- 
nale Dampfung, gekennzeichnet durch den Dampfungsfaktor (— f?/4 0”), an, so fiihrt dies auf 
das Variationsproblem (17) bis (19). Es sind also die U; die Zusatzverschiebungen eines solchen 
Ersatzkérpers, der unter denselben Belastungen steht wie der wirkliche Kérper, der aber fiir 
die Zusatzverschiebungen U; durchweg eine durch (— f?/4 07) gekennzeichnete ausschlagpro- 
portionale Dampfung besitzt. Man erkennt ferner, daB eine Gleichung (14) entsprechende 
Transformation der Gleichgewichtsbedingungen der Grundbewegungen uf) entstanden aus 
der Variationsgleichung (5), zu einem analogen Ergebnis fiihrt. Diese Gleichgewichtsbedingungen 
werden in der weiteren Rechnung jedoch explizit nicht bendtigt, lediglich die Erkenntnis, daB 
man zur Berechnung des Stabilitatsverhaltens des wirklichen Kérpers von einem Ersatzkérper 
ausgehen kann, der nicht nur fiir die Zusatzbewegungen die oben entwickelte ausschlagpropor- 
tionale Dampfung besitzt, sondern dieser in seinen sdmtlichen Bewegungen unterliegt. 

Die Problemstellung lautet jetzt also folgendermafen: Fiir einen Ersatzkérper, der denselben 
geometrischen Bedingungen und denselben Belastungen wie der wirkliche Kérper unterliegt, 
jedoch an Stelle einer geschwindigkeitsproportionalen Dampfung jetzt eine ausschlagproportio- 
nale Dampfung hat, sind die Stabilitatsbedingungen aufzustellen. Durch eine der Trans- 
formation (14) analoge Riicktransformation erhalt man dann aus ihnen die Stabilitatsbedingungen 
des wirklichen Kérpers. ; 

Jetzt lassen sich auch die GréBen Ky der Gleichungen (15) und (19b) deuten; sie sind die 
klassischen Spannungen dieses Ersatzkérpers fiir die Grundbewegungen u(), also mit den 
GréBen k) aus (8) nicht mehr identisch. 

Ren char muf werden, daB diese Uberlegungen auf der angefiihrten Arbeit von I’. Weiden- 


hammer aufbauen. 


3, Liésung des Variationsproblemes, Die Aufstellung der Stabilitatsgleichungen kann man mit 
kleinen Anderungen, die durch das Dampfungsglied bedingt sind, nach der von Mettler in M. II. 
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dargelegten Theorie durchfithren, wenn man fiir die 4uBeren Lasten wie dort die besondere An- 
nahme macht 

X,;=AXj+ eX; cos wt, | 

Pi=Api + epi’ cos at | 
wobei X; die Kriafte und p; die Driicke bedeuten, und worin die GréBen Xj, X;’, p;, pi’ reine 
Ortsfunktionen sein sollen. Weil die ki die klassisch berechenbaren Spannungen des Grund- 
zustandes des Ersatzkérpers sind, gilt bei ihnen das Superpositionsgesetz, und sie haben demnach 
die Form 


(i= 1, 2,3), 


KS) = Aol; + eoj) cos wt (i, 7 = 1, 2, 3), (20) 


wo oj; eine Funktion von Xj und pj, oj/ von Xj’, pj’ und @ ist. 
Das Variationsproblem wird auf direktem Wege mit dem Ansatz 


U, = Uy Wy Uy wy + U1, n—2Wn—2» 
U, = Ugg Wy + Ups Wz + + °* + Up n—1Wn—1> (21) 
U; = Ugg Ws} UggW_ + *** + Us n Wn | 


gelést, worin die u;; beliebig wahlbare, jedoch die Randbedingungen (11) erfiillende reine Orts- 
funktionen, die w; unbekannte reine Funktionen der Zeit t sind. Aus (18) und (19) erhalt man 
mit (21) 


L= T*¥— A* = yy (A fmf Wm— Dem WeWn— € C08 Wt+ Es, WeW,) - (22) 
ai 
Hierin ist 
D n= Bem A Cla ioe laa 


3 
1 
Am=z | { fo > U;spUim dV, (23) 
t= : 


3 
B EEC Oui f | Ouj f\ (OWim Ouj m | Ay Quif Oujm 
fm colo le I Ox; Ox; fe ax, i are On, Ox; dV, (24) 


und 


3 3 
Co. == als O:; OUy f OUpm GY 
maa fff S00 See es) 
i,j=1 pel 
ee : S7 Ouy f 8 
pi Se of, Uy Um an 
i 2 ie abs Ox; Ox; 26) 
Pil B= 


- 3 
H,,= 1 i eda: 
eel pres OE (27) 


In (23) bis (27) sind diejenigen u;;, welche im Ansatz (21) nicht vorkommen, gleich Null zu setzen. 
Aym, Drm und E;,, sind in den Zeigern f und m symmetrisch, auBerdem sind A;, und D,,, fiir 
jeden Kérper eigene feste Zahlwerte, waihrend Ey, vermége o;; Funktionen von sind. Die 


Ausdriicke 7 
A = DI Asm Wp Wry und D= mh DjmWs Wm 
f,m=1 f,m=1 
werden durch eine Transformation in Normalkoordinaten 
Ww; = @;40-+ Gig Ve + ie ~ Gin Uv, (== Ike Be On cn n) (28) 


iibergefiihrt in 
Ae, und D=») wv}. 
j=u =a 
Die oF erhalt man dann bekanntlich als Wurzeln der OF aus der Determinantengleichung 


|| Pm — OF Afm|| =O (f,m=1,2,...,n); (29) 
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sie sind die Eigenfrequenzen des nur mit 4X; und Ap; belasteten Ersatzkérpers. Mit der Trans- 
formation (28) erhalt man an Stelle der Epm (26) 


LS O; ¢ Oj mE; ;j (fom == Weedews ass nt) , (30) 


i,j=l 
welche wieder in den Zeigern f und m symmetrisch sind. Damit hat man jetzt fiir L (22) den Aus- 


druck 


L= 3 (uf — wf 7) — Ecos at: 3) Femvstms 
' f=1 fe 
und es sind die Euler-Lagrangeschen Differentialgleichungen des Variationsproblemes 6 | Lv dt==0 
also 
Up pee pe Come eS F102) (fail 25 ansjr0) 
m=1 


oder mit der Substitution T= wt 


d?v x 
(he _s CoFUF BOSE cee ra Mnias ” (fia recess 1h) s (31) 


Die 2n? Partikularlésungen dieses Systemes Mathieuscher Differentialgleichungen kann man in 
der Form 


Vir) op (t), Ve () =e ep (—2) (r= 1,2,--5m) (82) 


ansetzen, wo rf, (t) bzw. yy,(—T) periodische Funktionen der Periode 27 sind. Da nur die 
charakteristischen Exponenten 0, mit wachsender ,,Zeit‘*‘ tT unendlich groB werden kénnen, so 
sind diese maBgebend fiir das Stabilitatsverhalten. Beim Ersatzkérper liegt fiir ¢ = 0 in den 
Punkten 

|@r + ox| 
= = (k= 2. spe 0, sanz) (33) 
Instabilitat vor. Die Grenzkurven seiner Instabilitatsbereiche gehen also in einer (q, ¢)-Ebene 
von diesen Punkten w, der w-Achse aus. Wie schon aus friiheren Untersuchungen ! bekannt ist, 
liegen die Instabilitaétsbereiche des wirklichen Kérpers im Innern der Instabilitatsbereiche des 
Ersatzkérpers, enthalten jedoch die w-Achse nicht mehr. Man nennt die Instabilitatsbereiche 
bei @,) mit r= k Instabilitatsbereiche erster Art, die bei @) mit r +k Instabilitatsbereiche 
zweiter Art und p die Ordnung des Instabilitatsbereiches. 

In die weitere Rechnung fiihrt man mit 


QO 


0, = th, (34) 
die dimensionslosen Groen 
om Cae a ea ers 
oy tk eae a se oe e) 


ein, entwickelt die Fy,, in Potentenzreihen nach (@— @p) 


Fer, = Fm + (@ — 6p) Fm = (@ = @y)* Fyn -b * > (36) 
setzt 
FY oJ—?=6?, (j=0,1,2,...) (37) 
und entwickelt die periodischen Funktionen gy, in komplexe Fourier-Reihen 
ae ae ; 
Pfr (t) — OS Cfrn ent . 
n=—o 


Die Lésungen (32) lauten jetzt also 


++ © ; +9 : 
Vi, (x) = Dy Cry A (hy + n)t - Vi (r) = ay Ane i(h,+n)t . (38) 


y=>—o if iH 


Es geniigt, die Lésungen V;, zu untersuchen; die Lésungen V;/, bringen dann nichts Neues mehr. 


Setzt man V,;, (38) in die Differentialgleichungen (31) ein, so hat man mit (34) bis (37) durch Null- 


1 Vgl. E. Mettler, Mitt. d. GHH.-Konz., Bd. 8, Heft 1 (1940), S. 1. 
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setzen der Faktoren der e‘(',+)* die Gleichungen 
(a+ 79)? — 9] | 
ee = (Cm, req — 1 Gm, 1 + 1) [Gfn + Gfn (4 — 1) + Gm (Q- 1)? +.) Re 
re (fale Re, 798 = Ota Neer eeeerel | 
Man entwickelt die c;,,, q und q, welche Funktionen von ¢ sind, in Potenzreihen nach é: 
a oO +e on a e ae set 
q=lteq)+ &q%+-::, (40) 
fe agen ayonre es 


Setzt man diese Reihen in die Gleichungen (39) ein, so erhalt man durch Vergleich der Koeffi-: 
zienten von é fiir die ef) die Rekursionsformel 


ee n)? me qi] ¢ of) = es KY) (Gps OSE 3) (41) 
mit 
Kh, = 0, (41,0) 
K}), = — 2 (a +m) (Q® + 4g) + een vie RPE is) a 
ine fe s 
Ky) a ee ee) ea aula 2°) a De ere 

el B=0 
1% rae : 
Sa Gy (oe eee cl r.7 i y+ (41, j) 

m= 


DS ~ G aa Ved —t—1 : 
aR oe, ine (Cee p+ ef Ey Che eT oo) 


Die GroBen ee. stellen die Koeffizienten der ¢7 in der Reihe (q — 1) GY? + (q— 1)? Ge osee 
dar. Mit den Abe eth GréBen (35) gilt fiir die Punkte @, (33) 


qd =- qk = P (r wie oder Tr == k; P = 0, ganz). 


Fiir einen bestimmten Instabilitatsbereich sind also k, r, p fest vorgegebene, f, m, 7 laufende 
Zahlen. Nach (41) und (41,0) sind die ey) nur dann von Null verschieden, wenn 


CS (42) 
ist. Dies tritt ein fiir die Zeigerpaare 
(— ene a0) (=) mite ep (43) 


Fiir alle weiteren Zeigerpaare sei (42) verboten. Dannist fiir alle fund 7, welche die Bedingung (43) 
nicht erfiillen, 


p,=0- 
Setzt man 
cl ye —], 


und dividiert man die Lésungen V,, mit Cry,» So werden alle 
= 0 (>0), 
und die 
of) =e) (7 20) 


bleiben unbestimmte Integrationskonstanten. Nach diesen Festlegungen kann man aus den 
Gleichungen (41) fiir die Zeigerpaare f und 7, welche die Bedingung (43) nicht erfiillen, der 
Reihe nach die oy) berechnen, wahrend man fiir die Zeigerpaare f und 7, die (43) erfiillen, 


Ki) = a= (u Pree (44) 
K4) ven egos 3 (45) 


ieee 
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erhalt. Durch geeignete Wahl der q (und q(J) lassen sich diese beiden Gleichungen stets erfiillen, 
und Mettler nennt sie die beiden Hauptgleichungen. Bei ausfithrlicher Schreibweise erhalt man 
fiir diese Hauptgleichungen (44) und (45) 


| ae Gi< Pp); (46) 
Pegg tae wt oC nF ep)» (47) 
BG) + x» cli—p) (=p) (48) 
und 

| Ey (i< Pp), (49) 

Gt qk (0) _| 1 Ge 
2 ge pW) =) ae “(t DOCG) GS): (50) 

Pieters (enema 

| EQ) +. ci—P) aaa an one (51) 


Hierin ist x eine Konstante, die von qf und den Gt” abhangt; die BU) und E%) sind fir j <p 
Konstante, die von den Gh) (j= 1, 2,...) und gq; abhangen, fiir 7 > p Polynome von 1/c, c(, 


ce), ...,cG—P—). Zu bemerken ist, daB insbesondere fiir p= 1 gilt x = a (EP und da fiir 


‘Dee tammer bO—= EO) —_.0 ist. 
Die Reihe q (¢) = 1+ ¢ q+ ¢? q@) + - - - stellt in einer (g, ¢)-Ebene bzw. wegen (35) in einer 
(@, €)-Ebene eine Kurve dar, und in dieser Ebene ist fiir jeden Punkt (q, ¢) aus 


| 


er (e) = gt e gq) + & 24 


q()=—i 


ein bestimmter Wert der charakteristischen Exponenten 0, gegeben. 

Damit ist das Variationsproblem grundsatzlich gelést, und man kann beim Ersatzkérper fiir 
jeden Punkt einer (@, ¢)-Ebene angeben, ob er im Stabilitats- oder Instabilitatsbereich liegt. 
Denn nach den Lésungen (32) der Differentialgleichungen des Variationsproblems (17) kann dann 
instabiles Verhalten des Ersatzkérpers fiir ein Wertepaar (, €) eintreten, wenn fiir dieses einer 
der charakteristischen Exponenten 9, einen nicht verschwindenden Realteil besitzt. Fiir diesen 
kann dann V;, oder V;* mit wachsender ,,Zeit“‘ t beliebig groB werden. 

Nach der Transformation (14), dem Ansatz (21) und der Transformation (28) hat man als 
Zusammenhang zwischen den u; und den Zeitfunktionen vy 

pe ES 
ee OS! wr aye, (b= 1, 2,8); 
f,j=1 
wo wieder die in (21) nicht vorkommenden u;; gleich Null zu setzen sind. Da nun die u, in den 
vy und die vein den Vy, und V/, linear und homogen sind, so kann man den Exponentialfaktor 
exp (— {1/2 o@) in den Lésungen V;, und Vj, (32) zu deren Exponentialfaktoren exp (0,7) baw. 
exp (— 0,7) hereinmultiplizieren, und es herrscht demnach fiir den wirklichen Kérper erst dann 
Instabilitat, wenn exp (+ 0,— 6/2 0m) tv mit wachsender ,,Zeit’ + unendlich gro werden kann. 
Nach dem oben Gesagiten ist dies erfiillt, wenn 
B 2 

R(e)|— 52, >0 (52) 
ist, wo R(o,) den Realteil von o, bedeutet. Mit (34) und den dimensionslosen G1éfen (35) er- 
halt man fiir diese Bedingung (52) 


Jl >e, (53) 
wo J(q) den Imaginarteil von q bedeutet. 

4, Die Grenzkurven der Instabilititsbereiche. Nach der Ungleichung (53) kann fiir den 
betrachteten Kérper Instabilitat auftreten, wenn q= q,-+ €q")-} e?q-4- ++ einen Imaginar- 
teil besitzt, dessen Betrag gréfer als C ist. Fir € = 0, was f = 0 (keine Dampfung) entspricht, 
tritt demnach Instabilitat fiir die Werte (q, ¢) auf, fiir die q reell und q komplex ist. Setzt man 
die willkiirliche Integrationskonstante c() als komplexe Zahl 


cl — R(cosp + isin@) 
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an, so erhalt man nach (47) und (50) mit R= }q,/q, und beliebigem p = 0 
J (244 pg”) = 0, 
q@) =/Zsing. 
J (2 qr gq?) |/% sing 


Demnach ist hier q(®) reell, g) komplex, und fiir € = 0 kann Instabilitat auftreten. Nach M. II 
sind die Koeffizienten g{j) der Grenzkurven qc der Instabilitatsbereiche des Ersatzkérpers ge- 


geben durch 


: (J) A 
gy) a (i< Pp), (54) 
24% P 
(p) “ J 
(pP)__ ES = (55) 
qe) = (j= P)> 
Caan p.  pigeds 
; EW) : 
qi) — © (j >Pp)- (56) 
24q%P 
Tn (56) ist ED) der Wert von E(), wenn man in ihm cj) = ef) mit 
Oe ye (57) 
dk 
; 1 3/2 a 77(j) 5 
c(i—P) ee ees whe OE ——s il 2 Saatts 58 
i a aE GaP th pte.) (58) 


setzt. Wegen beider méglichen Vorzeichen in (55), (57) und (58) erhalt man damit fiir jeden 
Instabilitatsbereich des Ersatzkérpers, gekennzeichnet durch die zwei Zeiger k und r, zwei 


Grenzkurven 
ic=lted+erqQh+---. (59) 
Um zu den Grenzkurven des Kérpers mit € + 0 zu gelangen, bildet man eine den betreffenden 
Instabilitatsbereich des Ersatzkérpers iiberdeckende Kurvenschar und bestimmt auf jeder dieser 
Kurven der Schar denjenigen Punkt (qo, €9), in welchem fiir ¢ > ¢) Instabilitat, fiir ¢ < e) Sta- 
bilitat herrscht. Die Gesamtheit dieser Punkte bildet die Grenzkurve des entsprechenden 
Instabilitatsbereiches des wirklichen Kérpers und liegt noch ganz im Stabilitatsbereich des 
wirklichen Kérpers. Diejenige Kurvenschar 


G = IMGs 2) = 8 gee ice eb) Peg ee (60) 
erfiillt in 0< m< a diese Bedingung, in welcher die Koeffizienten qi) fiir j + p mit den q{) 
von (54) und (56) identisch sind, in q'?) jedoch 


() — ye (cosp + ising) (61) 
ist, man also : 
E®) ve, 
(Pp) — 
q cos 
2%? | plage 


hat. Aus der Forderung q\) — q{) fiir j ++ p und mit c() aus (61) erhalt man die ci —?) nach (51) 
zu 


Gasp ep 2. )e (62) 


Fiihrt man diese Werte c(i— P) von (61) und (62) in die q(/) (7 = p) von (47) und (48) ein, so erhalt 
man fiir 0 << y< a bei q einen nicht verschwindenden Imaginiarteil. Die Kurven f (9, ) (60) 
liegen also alle im Instabilitatsbereich des Ersatzkérpers. Ferner gilt ¢(°) = c!) fiir p = 0 oder a, 


und durch Bildung des Grenzwertes erhalt man 
Gap) — AUS 
Cima) . 0) =&g 
C 
Es schwenken somit die Kurven (60) fiir cl) > cl) d. h. fiir p — 0 oder z stetig in eine der beiden 


qc-Kurven (59) iitber. LaBt man also @ stetig von 0 nach a gehen, so schwenkt die Kurve (60) 
stetig durch den Instabilitaétshereich des Ersatzkérpers. 
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Nach (53) gilt fiir Stabilitat |J (q)| <€, fiir Instabilitat |J (2)| >¢. Nimmt man fir |J (q)| 
einen festen Wert an und variiert ¢ stetig, so liegt bei |J (4)| = ¢ der Grenzpunkt, bei dem gerade 
tia Stabilitat herrscht. Setzt man also c@) (E> 0) aus (61) und (62) in (47) und (48) ein, so 

ommt 


q(%, 2) = ge g® 4 e7 g@) os, 


und man erhalt den gesuchten Wert ¢, der Grenzpunkte (qo, &,) fiir jede durch ein festes p aus 
0< 9) < % definierte Kurve (60) aus der Gleichung 


IJ (4 [Po €ol)| = 6. (63) 


Hieraus flieBen méglicherweise fiir ein gy mehrere Werte ¢. Derjenige Wert ¢, ist zu nehmen, 
fiir den |J (q [p. €o|)| ein Minimum besitzt. Diesen Wert ¢, und den entsprechenden Wert Po 
setzt man in den Koeffizienten q'?) von (60) ein und erhalt damit 


(Pp) 

eek ga 3 oes eet) a (64) 
24hP PV9r % 

Damit sind die Grenzpunkte (qo, ¢9) bestimmt, und man kann die Grenzkurven des betreffenden 

Instabilitatsbereiches des wirklichen Kérpers mit Dampfung punktweise in ein (q, ¢)-Diagramm 

einzeichnen. 

Die Gleichung (63) hat keine Wurzel ¢, = 0. Es gibt deshalb in jedem Instabilitatsbereich 
einen Minimalwert é€ 9 nin > 0, und es kann fiir keinen Wert q in diesem Instabilitatsbereich fiir 
E9< €gmin Instabilitat auftreten. Diesen Wert &9;, nennt man den Schwellwert, und man be- 
rechnet ihn, indem man q)= 2/2 setzt. 


5. Beispiel. Im Anschlu8 an die Theorie sollen die Kippschwingungen eines Stabes mit recht- 
eckigem Querschnitt berechnet werden (Abb. 1). Mettler hat in seinen Arbeiten ebenfalls dieses 
Beispiel gewahlt, und es sollen hier zu Vergleichs- 
zwecken dieselben Abmessungen und Material- 
konstanten verwendet werden. Um die gewohn- 
ten Gleichungsformen zu erhalten, gilt wie bei 
hei Crea = Xs, Ky, Se, 

Die Abmessungen des Stabes sind: Lange 
l= 100cm, Héhe h=10cm, Breite b= 1 cm. 
Ferner sei der Elastizitatsmodul E = 2,1 - 10° ; =i 
kg/em2, ere seh ahaa Coe 8,1 - 105 kg/em?, Abb. 1. Bezeichnungen fiir das Beispiel. 
die Masse 0 = 8,0- 10-* kgsec?/cm*, also die Masse der Langeneinheit «4 — 8,0 - 10~> kgsec?/cm?. 
Als Tragheitsmomente berechnet man mit den iiblichen Bezeichnungen J,— 0,8333 cm4, 
J,= 83,33 cm*, Jz —=3,130 cm! und als Tragheitsarm des Querschnittes i, = 2,901 cm. Der Stab 
hat an beiden Enden Gabellagerungen, so daf er sich dort nur in x- und y-Richtung bewegen 
kann. Die Belastung besteht aus Endmomenten 


Mr; =AM,-+ €M, cos ot, 


welche den Stab in stabilem Zustand in y-Richtung ausbiegen. Fiir die Rechnung wahlt man 
1 
4 My)=~ Mi. M,= M,,. 
wo M,— 66170 cmkg das Kippmoment des Stabes ist. 
Nach M. I. macht man hier fiir die U; den Ansatz 
Ua = te we OY, 


wo v (x, t) und w (x, t) die zusitzlichen Ausbiegungen der Stabachse in der y- und z-Richtung, 
und # (x,t) der Verdrehwinkel der Querschnitte um die Stabachse sind. Als Spannungen des 


Grundzustandes hat man 
_ My _ OM 
ie i und [ [rae aa 


mit M als Biegemoment. Das Doppelintegral exstreckt sich ittber den Stabquerschnitt. 
Mit Beriicksichtigung der ausschlagproportionalen Dampfung erhalt man fiir den Ersatzstab 
entsprechend den Gleichungen (18), (19) und (20) [vg]. M. II., Gl. (15) bis (18)] 
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1 


1 
x Ow)\2 wo [OD\2) [ee te 2 “2 99 
ie 5 ites mus ) Cie ae (w? + 15 0) dx 
0 0 
1 ! 
» __ Edy [ (ew\2 |, GJa f foo? 
Aitaes = 2 [(e) dx | 2 ie dx , 


0 


A*,, = [ea o (M9) dn, 


eke 
oO sin Te 
M=AM,-++ —— cos wt’ =f Rk Way ake) 
ashy 
ky 


wo >’ bedeutet, da nur iiber ungerade k aufsummiert wird. Die @, sind die k-ten Kigenfrequen- 
zen der Biegeschwingungen des Ersatzstabes in y-Richtung, und es gilt 


Oe km\4 E Jz Bp? 
oty=(T] emg (65) 
Nach M. I. erhalt man hier also fiir die mike c wake Soda (15) 
O4w o2 2 F2 " 
Fla woke) a Fe Berge! ale 


O20 oy g FP? 


aro 0?w : 
CJizat Moab an lp i reek 


Als Randbedingungen (16) hat man fiir Gabellagerung 


ae ae OG (bie Ge=—) ciel oe == i. 
Ox 


Nach M. II hat man fiir den Ansatz (21) 


w= wysin =p xt oysin xf ay _asin 2 a, | 
oo (66) 
9 = wy sin x4 wy sin 22 x4 +--+ a, sin @ ay, | 


Nach Mettler kann man namlich bei diesem Ansatz die geraden Koordinatenfunktionen sin (kz x/1) 
mit geradem k weglassen, da die gesamte weitere Rechnung in zwei analoge Halften mit ge- 
radem k und ungeradem k zerfallt, welche beide dasselbe Ergebnis liefern. Durch Hereinnahme 
des Dimpfungsgliedes andert sich an diesen Uberlegungen nichts. 


Fiir die Ausdriicke D und A erhalt man 


n—l 
EP eat ko? F22 Be) F? 2B 
Paani k || EJs(=) saat” wi + 2A Mou, uy. 1 + \|6Sq hap a Weis 


Der Transformation (28) entspricht 


2 
w, = (vj, COSA,—Vy_ + 1 SIN Oy) 
Jui | 
| (67) 
Wy 4 1 == (% sin a, + vy 4 1 COs Oo) | 
ip)ul 
mit 
tg D) 0, = 21 M, Up : 
kip 
EI] ) —G Ja 
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jeweils fiir k — 1, 3,5,...,n—1. Die Determinante (29) wird damit 


ka\4 Oy Re 
aGey Eon] aa 
fetes _|=0. (68) 
kn\? Pe Gaia 
f amlisy fealty PPE nieg 


| 


Ferner ist, wenn man mit 4A den mit ¢ multiplizierten Anteil von A*,, bezeichnet, 


1 


a ; n—l1 S as 

= 4 /a\2 ie _ bx ka ee 

A=eM, cos ot (7 De ke? essa fn in oe oe ae dx. (69) 
ik=1 on 


| es) 


Das logarithmische Dekrement der Schwingungen ist hier 6 = {2/2 0m, und zwar bleibt 
fiir Metall! stets 6< 10-%. Fiir die nachher berechneten Instabilitatsbereiche ergibt dies 
B< 10-4kg sec/em*. Ein derart kleines f andert die Eigenfrequenzen (65) und die Wurzeln der 
co; in (68) nur um vernachlassigbar kleine Werte, so das man fiir die @,, und die m,die Werte von 
Mettler tibernehmen kann, namlich 


Oyy= 1460 sec}, 3 y= 13 140 sec 


und 
@, = 126,3 sec, w,= 1294 sec=, 


@, = 1929 sec, @,— S180 secu.. 


Fiir die Instabilitatsbereiche w)= @, und @)= @,-+ @,ist bei k > 3 


so daB man nach M. II. 


fo) sin —- % Py 
= a i | Z(w) sin = x 
aap 
Mk y | 
mit 
] 
Z(@) 


2 
My 


setzen darf. Rechnet manim Ansatz (66) bis n = 4 (Mettler hat gezeigt, da dann die Reihen (40) 
bei 6B = 0 schon gut konvergieren), so erhalt man statt (69) 


2 
ze £08 Mt (WwW, W. + 9 w; w,) + 


16 M. 1 1 3 81 
+ Z(@)— *€ cos ots Wy Wy — zx Wy Wy — FH Wy We rae ty tt). 
und die Transformation (67) wird zahlenmahig 
oe = (0,9994 v, + 0,0380 »,) , 
Ve 


2 


W, = ——=(— 0,0380 v, + 0,9994 v9) , 
ip Vl 

w, = = (0,9992 v, + 0,0389 »,), 
\ul 

w= —“= (— 0;0389 2, + 0,9992 v,) . 
ip ul 


Daraus berechnen sich die Fs, nach der Vorschrift (30), und aus ihnen fiir jeden Instabilitats- 
bereich die ey nach (36) und (37). 


1 Vel. Bennewitz u. Roétger, Physik. Z. 37 (1936), S. 578. 
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a) Instabilitatsbereich erster Art, zweiter Ordnung. Es gilt hier 
Wo = @, = 126,3 sec™ 1, 
k=r=1, p=2, g=¢=—q=1, 


also 
c(%), = 1. ci) c= 0 J = 1) ° 


cy) = (G20). 
Aus den Gleichungen (44) und (45) ergeben sich sukzessiv mit den aus (41) berechneten 
| mal 5 60, firey anes 
1 6 
“nin GPa fy oa ed at 
mal +-G0), fir y=-+ 1, 
I l ea) (0 Poe D 
| mal + GO, e+ g@ 2 GO), e() fir n=—3, 
1 l eo a y+ ¢a 0 ee 
oe mal CO). ¢)-— 9 GO) (oO) oy) ae eee 
oe pe as 2 mi T' 2 aie ceed Bey lbs 
1 
| mal — G0) fir 7=-+1 
die q) und q() zu 
qo. 
q) == 0 J 
g°) = — 0,4474 — 0,3906 o() 
1 
g) = — 0,4474 — 0,1958 (< ae zi) 
7) = — 0,3906 c), 
: iL 
aie ow On = 
g®) 0,1958 ¢ (1 =) 


Damit kommt fiir die qg-Kurven (59) 
qc = 1— (0,4474 + 0,3906) 2 + 0 - <3, 
fiir die Berechnung der ¢, nach (63) 
0,3906 sin gy - e6 + 0-25 —C=0 
und fiir die Berechnung der gq, nach (64) 
J = 1 — (0,4474 + 0,3906 cos gy) e8 + 0+ &2. 


Hier ist C/8 = 494,9 cm4/kgsec. In Abb. 2 sind die Grenzkurven fiir 6 = 10-4 und 10-1 kgsec/cm* 
sowie die gco-Kurven, also die Grenzkurven fiir B= 0; eingetragen. 


b) Instabilitatsbereich zweiter Art, erster Ordnung. Es gilt hier 
Wo = @,-+- @.= 2055 sec“, 


[han bee ee 5 ese q. = 90,0614, gq, = 0,9386, 
und man erhalt 


9) = — 2,088 - 10-2 e , 

g) = — 2,088 - 10-* e + 3,190 - 10-1 1, 

gq) = — (0,1872 e + 2,088 e) — 4,382) 10-2, 

Oe (121.1 van ee = + 2,088 () + 0,1872 c()* +. 1.749) 10-2, 
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Berechnet man q@) und q@) und setzt die Werte c() und ci —») aus (61) und (62) ein, so ist fiir 
y = 16°, was e = 0,5 entspricht, is (q°))| <10~*. Die Rechnung zeigt, daf man jeweils die 
Glieder mit ¢? vernachlassigen darf, und es kommt somit fiir die qc-Kurven nach (59) 

qc = 1 + 0,1632 ¢ — 0,1253 ¢?, 
fiir die Berechnung der ¢) nach (63) 


0,02088 sin @, + &)-+ 1.368 Se eee ee 


— cos 2g 1 —cos 29} ° 
€ 
O¢ 
b=10% 
G3 
p-0 
02 
G7 = 
b=" 
« 
0 1__ 
| O85 G90 O95 700 r G90 Q9S 400 105 . 
eee ee ee a an eee See ee l : 
105 710 5 720 105 tab 3sek™' 7800 7900 2000 20853 2100 220088K * 
Abb. 2. Instabilitétsbereiche bei m) = @.- Abb. 3. Instabilititsbereiche bei wy) = w 1+. 


und fiir die Berechnung der q, nach (64) 
GU 0,1632 cos gg* ey—0, 1259 FF. 
Hier ist C/6 = 30,415. In Abb. 3 sind wieder die Grenzkurven fiir 6 = 10-4 und 10-5 kgsec/em* 


sowie die gc-Kurven eingetragen. 


(Eingegangen am 14. November 1951.) 


Anschrift des Verfassers: Dr. Giinther Grammel, i. Fa. Schnellpressenfabrik Koenig & Bauer A. G. 
(13a) Wiirzburg 7. 
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Zur Synthese der Schiebepaar-Getriebe. 
Von K, Hain, 


1, Einleitung. Die Getriebe mit Schiebepaaren haben die gleichen willkommenen Figen- 
schaften wie die Getriebe, die nur Zylinderpaare enthalten. Ein Schiebepaar besteht bekanntlich 
aus einem Gleitschuh, in welchem sich der Gleitschieber ohne Drehen hin- und herbewegen kann. 
Da heide Teile sich gegenseitig in Flachen beriihren (niederes Elementenpaar), kénnen solche 
Schiebepaare ebenso wie die Zylinderpaare wesentlich starker beansprucht werden als Paarungen 
mit Linien- oder nur Punktberiihrung (héhere Elementenpaare), wie diese z. B. bei Kurven- 
getrieben angewandt werden. Wegen der geringeren Abnutzungen dndert sich das Bewegungs- 
gesetz der Schiebepaar-Getriebe praktisch itberhaupt nicht, und es kénnen beliebig groBe Krafte 
und Momente beihohen Geschwindig- 
keiten iitbertragen werden. 

Die Synthese der Schiebepaar- 
Getriebe, d.h. deren Entwurf fiir 
vorgeschriebene Bedingungen, wurde 
von R. Beyer1, auf den Burmester- 
schen Grundlagen aufbauend, im 
eigentlichen Sinne begriindet. Im 
vorliegenden Aufsatz sollen nun 
auch hierfiir die vom Verfasser ent- 
wickelten Punktlagenreduktionen 
verwendet werden. Mit ihrer Hilfe 


= IRS 

a Sy 

7 FOS Vy 
Ai > we 


Ag 
Pe “2/ itich Na Ve 


| 


Abb. 1. Konstruktion einer Kurbelschleife fiir Zwei- Abb. 2. Konstruktion einer Kurbelschleife fiir g,,= yp». 
und Drei-Winkelzuordnungen. 


kénnen unter bestimmten Voraussetzungen héhere Anforderungen erfiillt werden als dies bisher 
méglich war. Es sollen zunachst nur die Getriebe der Viergelenkkette behandelt werden. 


2. Winkelzuordnungen, Winkelzuordnungen sind dann zu erfiillen, wenn die Bewegungen 
zweier im Gestell gelagerter Glieder einander zugeordnet werden sollen 2. So sollen in Abb. 1 die 
Winkel ¢, und 3 des Gliedes A) A den Winkeln y,, und y,, des Gliedes B, B, das eine Gleit- 


fiihrung B haben soll, entsprechen. Man tragt im Gestellpunkt A, den Winkel 5% undin By 
den Winkel = Wy, an Ay By an. Der Schnittpunkt A, der freien Schenkel beider Winkel bestimmt ” 


die Kurbellange A,A. Mit Hilfe der gegebenen Winkel g,. und y,3; hat man auch die Kurbel- 
lagen A) A, und A, A;. Nunmehr verschwenkt man die Punkte A, und.4, um By um die Winkel y,5 
und 7,3in dem Sinne, wie sich das Glied B, B von den Stellungen 2 und 3 in die Stellung 1 be- 
wegen soll, und man erhalt die Punkte A,, und A,,, wobei gema® den vorbereitenden MaB- 
nahmen A,, mit A, als Punktlagenreduktion zusammenfiallt. Man braucht nun nur noch die 
Gleitfiihrung B als die Gerade A, A, auszufiihren und hat die gestellte Aufgabe mit Hilfe der 
so entstandenen Kurbelschleife erfillt. 

Ist noch ein drittes Winkelpaar p,— 4 zu verwirklichen, so liegen nach den gleichen 
Konstruktionen die drei Punkte A,= A,,, A, und A,, im allgemeinen nicht auf einer Geraden. 


* R. Beyer, Zur Synthese ebener und raumlicher Kurbeltriebe. VDI-Forschungsheft 394, Berlin 1939, 
2 K. Hain, Instrumentenkunde 63 (1943), S. 170. 
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Da aber in vielen Fallen nicht die absoluten Winkel y oder w, sondern nur deren Verhiltnis- 
werte vorgeschrieben zu werden brauchen, kann man die Konstruktion so lange wiederholen, 
bis eine Gerade durch die drei Punkte zustande kommt. 

Kine Kurbelschleife entsteht auch dann, wenn zwei zuzuordnende Winkel g und yw gleich 
groB sind (Abb. 2). Der Schnittpunkt A, der freien Schenkel der in 4, und B, an A, By an- 
getragenen Winkel _ Py. und + Ye liegt im Unendlichen, da diese parallel zueinander sind. 


Die Kurbel 4) 4 wird hier ersetzt durch eine Gleitfithrung g,, die in A, auf dem freien 


Ahb. 3. Punktlagenreduktion an der Geradschubkurbel. 


Schenkel von Yj2.senkrecht steht. Die Lagen g, und g, liegen symmetrisch zum Gestell, wahrend 


die Lagen g, und g, mit g, die Winkel g,, und g,, einschlieBen. Verdreht man die Gerade g, um B, 
um den Winkel y,.—=q .in dem Sinne, wie sich die Schwinge By, B von der Lage 2 in die Lage 1 
drehen soll, so fallt die so gedrehte Lage {omit g, zusammen; esist also eine ,,Geradenreduktion“ 
entstanden. Die Drehungen der Geraden g, und g, um By um die Winkel y,3 und y,,im gleichen 
Sinne ergeben die Lagen g,, und g,,. 
An die drei Geraden: g,== 8.4. £315 a 
&,, kann man den beriihrenden Kreis pega — Nee + 
k mit seinem Mittelpunkt B, zeich- i | | 
nen. Mit diesem Punkt sind die Ab- | 
messungen der Kurbelschleife be- 
stimmt, wenn man die Geradfithrung ve 
& parallel zur Geraden g, legt. Wenn 
sich also die Schwinge B,B um die Ay _| £2 

: s | 
vorgeschriebenen Winkel yw dreht, loops We 
wird der Gleithebel A)g, um die Rit oa ee | | 

. 4. Punktlagenreduktion 

ebenfalls vorgeschriebenen Winkel p an der Kreuzschleife. - 
verschwenkt. Fiir dieses Beispiel der 
Kurbelschleife kann also eine Lagenreduktion nur herbeigefiihrt werden, wenn zwei zuzuordnende 


Winkel » und y einander gleich sind. 


se 


UZ >| 


3. Weg-Winkelzuordnungen, Mit Hilfe des Geradschubkurbel-Getriebes kénnen vorgeschrie- 
bene Winkel der Antriebskurbel vorgeschriebenen Wegen des Gleitsteines zugeordnet werden. 
In Abb. 3 sollen den drei Winkeln ,9, (13, Py, die geraden Bewegungsabschnitte nj., m3, Ny4 eines 
Gleitsteines entsprechen. Fiir den Gleitstein B ist die Richtung und der Richtungssinn von B, 
nach B, vorgeschrieben, nicht aber eine Punktlage B,, B,, B, oder B, beziiglich des Drehpunktes A) 
der Kurbel A, A. Man zeichnet beliebig in der Zeichenebene beispielsweise die Strecke A, A;—= Nog 
und konstruiert iiber dieser Strecke als Grundlinie das gleichschenklige Dreieck mit dem Spitzen- 
winkel g3 und dem als Kurbeldrehpunkt dienenden Spitzenpunkt A,. Mit Hilfe der gegebenen 
Kurbelwinkel ¢,, und y,, kann man nun auch die Kurbellagen 4, A, und A) A, aufzeichnen. Man 
]aBt nun die Punkte A,, A, und A, um die Strecken ny, M43, M,in dem Sinne wandern, wie sich 
der Gleitstein von B,, B,, B,nach B, bewegen soll. Die Richtung dieser Bewegungen ist durch 
die Strecke A,A, gegeben. Man erhalt dabei die Punkte A,,, A3,, Ay,, wobei A, mit A,, als 
Punktlagenreduktion zusammenfallt. Der Kreis durch die drei Punkte A,, A,;= Asis Aix 

14 
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ergibt als Mittelpunkt den Gleitstein-Gelenkpunkt B,, von dem aus die Gleitfithrung in der 
vorgegebenen Richtung, also parallel zu A,A,, ausgefiithrt werden kann. Es ist natiirlich zu- 
lassig, fiir jeden der sechs méglichen Bewegungsabschnitte n mit dem dazugehorigen Winkel 
eine Punktlagenreduktion zu erzwingen. Auch doppelte Punktlagenreduktionen sind hier még- 


Abb. 5. Punktlagenzuordnungen fiir eine Geradschubkurbel, wenn 
die Verbindungsstrecke A, A, senkrecht auf dem Strahl B, X, steht. 


lich, wenn voneinem der zuzuordnenden Werte 
n oder y nur deren Verhiltnis vorgeschrieben 
ist. Es sind dann also Vierwinkel-Vierweg- 
Zuordnungen erfiillbar. 

Fiir die Geradschubkurbel von Abb. 3 
tritt ein Sonderfall ein, wenn die Relativ- 
punkte A auf einer Geraden liegen. In Abb. 4 
soll fir die Lagen 3 und 4 eine Punktlagen- 
reduktion zustande kommen. Es ist A, 4,==N4,. 
La®t man die Punkte A, bis A, in Richtung 
A, A,um die Strecken n wandern, so fallt A,, 
mit A,, in einen Punkt zusammen. Da die 
drei Punkte A,, A,,, As, = Ay, auf einer Ge- 
raden liegen, stellt diese Gerade einen Gleit- 
schieber g, dar, in dem ein im Punkt A an- 
gelenkter Gleitstein gleitet. Das Getriebe ist 
damit eine Kreuzschleife geworden. 


4, Erzeugung gegebener Kurven, Wenn 
eine gegebene Kurve als Koppelkurve z. B. 
eines Geradschubkurbel-Getriebes erzeugt 
werden soll, kann man eine Anzahl von Punk- 
ten auf dieser Kurve annehmen und das Ge- 


triebe so konstruieren, da diese Punkte von einem Koppelpunkt des Getriebes genau durch- 
laufen werden. Zwischen diesen Punkten wird dann die Koppelkurve von der verlangten Kurve 
um so weniger abweichen, je mehr man genau zu durchlaufende Punkte annehmen kann. 


Abb. 6. Punktlagenzuordnungen fiir eine Kurbelschleife, wenn die 
Punkte By, By, Bo; auf einer Geraden liegen. 


In Abb. 5 sollen die fiinf Punkte E, bis 
E, auf der gegebenen Kurve nach praktischen 
Erwagungen ausgewahlt worden sein. Man 


Abb. 7. Punktlagenzuordnungen fiir eine 
Winkelschleife. 


bringt die Mittelsenkrechten e,; und e,, der Verbindungsstrecken zweier beliebiger Punktpaare 
FE, E; und E,E, zum Schnitt By und erhalt die Winkel y, und Wea2- Mit By als Ursprung 
zeichnet man dann ein beliebig in der Zeichenebene liegendes Strahlenbiischel mit den 


* K. Hain, Maschinenbau-Betrieb. Beilage Getriebetechnik 11 (1943), S. 29. 
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: 1 1 
Winkeln < X,)B, X, = > Vs und) = Xb, X.== 9 Ysa Mit beliebigem, aber jeweils gleich 


groBem Halbmesser schlagt man Kreisbigen um E, und E,, welche die Strahlen B,X, und 
B,X; in den Punkten A, und A, schneiden. Dies wiederholt man so lange, bis die Verbindungs- 
strecke A, A, senkrecht auf dem Strahl B,X, steht. Den Punkt A, findet man als Schnitt- 
punkt der Geraden A, A, mit dem Kreisbogen um E, und E, A, als Halbmesser. Die Punkte 
A, und A, liegen auf der Geraden 
A, A, symmetrisch zu A, und A, 
und haben von E, und E; die 
Entfernung E,A,. Macht man: 
fi E A, By /\ E, A,B, und 
A E, A, Byy= A E, A; By, s0findet 
man die noch fehlende Gelenk- 
punktlage B, als Mittelpunkt des 
Kreises durch die drei Punkte B,, 
By, Bos. Wiirde man die Dreiecks- 
konstruktionen fiir die Punkte F,, 
A, und E,, A; durchfiithren, so 
wiirden als Punktlagenreduktio- | 

nen die Punkte B, mit B,; und | AG | II 
By; mit By, zusammenfallen. Der | 

Gestellpunkt B, ist wegen der be- 
schricbenen, vorbereitenden MaB- “*-& Tonfucs Puntulgenserdngngen, for sng Cornea 
nahmen gleichzeitig der Pol P,; 

und P,, der Koppel AB. Es ist klar, daB man auch noch sechs E-Punkte auf der gegebenen 
Kurve annehmen kann, wenn man die Lage des X-Strahlenbiischels so lange verandert, bis vier 
B,-Punkte auf einem Kreis liegen. AuBerdem gibt es fiinf Méglichkeiten zur Paarung der finf 


E-Punkte, so daB auch dadurch noch v7) 

bestimmte Anspriiche beziiglich des ys A é 

Platzbedarfes und der Ubertragungs- f- Nas ff 

verhaltnisse erfiillt werden kénnen. wets ea | ee 
Fiihrt man in Abb.6 die Konstruk- | ay | 


tion wie in Abb.5 durch und findet am ; 
SchluB, daB die drei Punkte By, By, Nika 
By; auf einer Geraden liegen, so erhalt | 
man eine Kurbelschleife mit der Gleit- 


fiihrung A,B, welche parallel zu der | re 

Geraden durch die drei erwahnten i 

Punkte liegt. Man kann fiir das be- fae 

liebig in der Zeichenebene liegende AAS dis 


X-Strahlenbiischel den Halbmesser so Aaa 
lange verandern, bis die drei By-Punkte / 


auf einer Geraden liegen. ys 
In Abb.7 entsteht eine Winkel- sf 
schleife, wenn sowohl die Punkte A, und Vi, 


A, auf einer Geraden senkrecht zu By Xo ae 
Re nee tee fare i arte, Uaeecasnes nega. fo iaee Conn shebh abel 
einer Geraden liegen. Bei der Forderung 

nach einer Winkelschleife mu8 man fiir jede Lage des X-Strahlenbiischels den Halbmesser E, A, 
so lange verandern, bis die Gerade A, A, senkrecht auf By) X, steht. Daraufhin muf man die 
Lage des X-Strahlenbiischels unter Beibehaltung seines Ursprunges B, verandern, bis auch die 
drei B,-Punkte auf einer Geraden liegen. Es leuchtet ein, daf fiir eine Winkelschleife nach der 
angegebenen Konstruktion im Héchstfalle fiinf Punktlagen auf einer gegebenen Kurve erfiillt 
werden kénnen. 

Ein Sonderfall bei den Punktlagenzuordnungen entsteht, wenn die Verbindungsgeraden 
zweier Punktlagenpaare, die zu Punktlagenreduktionen fiihren sollen, einander parallel sind. 
Dies trifft beispielsweise in Abb. 8 fiir die Punktlagenpaare E,— F, und E,— E,zu. Der Punkt B, 
als Schnittpunkt der hier parallelen Mittelsenkrechten auf E, E; und E, E, liegt im Unendlichen. 


14* 
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also bewegt sich der Punkt B auf einem Kreisbogen mit unendlich groBem Halbmesser, d.h. es 
muf eine Gleitfiihrung vorgesehen werden, die parallel zu den Strecken E,E, und EF, E, ist. 
Uber die Lage dieser Gleitfithrung ist dabei noch nichts ausgesagt. Das X-Strahlenbiischel, 
das seinen Ursprung im Punkt B, hat, besteht, da dieser Punkt im Unendlichen liegt, aus den 


1 I : 
parallelen Geraden X,, X,, X,, wobei die Abstande X)X,— = ™s> 9, CD, > Mea sind, wenn 


ny; —= E, E, und n,,— E, E, festgelegt werden. Im itbrigen kénnen die X-Geraden mit ihrer Rich- 
tung senkrecht zu FE, E,; und E, E, beliebig in der Zeichenebene liegen. Man schligt um E, und E, 
Kreisbégen mit beliebig groBem, aber gleichem Halbmesser, welche die Geraden X, und X, in 
den Punkten A, und A, schneiden. Die Mittelsenkrechte a,, auf der Strecke A, A, schneidet 
dann die Gerade X,im Punkt A,, und mit A, A, —A,A,ist die Kurbellénge der im Gestellpunkt Ay 
gelagerten Kurbel gegeben. Die Punkte A, bis A, findet man als Schnittpunkte des Kurbel- 
kreises mit den Kreisbigen um E,, E,, E; mit E,A, als Halbmesser. Die Punkte A; und A, 
miissen dabei symmetrisch zu den Punkten A, und A, beziiglich der Geraden X, liegen. Fiir die 
Konstruktion des Getriebes ist es gar nicht notwendig, ihre Lage zu ermitteln. Nunmehr be- 
stimmt man zwei beliebige Pole, die eine gleiche Indexziffer haben, z. B. die Pole P,, und P,3, 
die beide die gleiche Indexziffer 1 haben. Den Pol P,, findet man als Schnittpunkt der Mittel- 
senkrechten a,, auf A, A, mit der Mittelsenkrechten e,, auf E, E, und den Pol P,, als Schnittpunkt 
der Mittelsenkrechten a,, auf A,A, mit der Mittelsenkrechten e,, auf E,E,. Damit sind auch die 
Drehwinkel 6,,— x E,P,,E, =< A,P,,A, und 6,:=<X E,P,,E,=<X A, P,,A; bekannt. Andie 


durch P,, und P,, gehenden, senkrecht zu E,E,; stehenden Geraden hy, und h,, tragt man in dem 
Sinne, wie E, zu E, und E, liegt, die Winkel = Oy. und = 6,3 an, deren freie Schenkel sich in B,, 


der Punktlage 1 des Punktes B schneiden. Mit B, ist die Koppellange A,B, und die Lage der 
Gleitfiihrung gegeben, die, wie bereits erwahnt, parallel zu E, FE; liegen muB. 

Man kann auch eine Geradschubkurbel fiir sechs Punktlagen konstruieren, wenn die Ver- 
bindungen dreier Punktlagenpaare einander parallel sind, wie in Abb. 9 die Strecken E, E,= 1g, 
E, E;= no5, E, E,=:nz,. Das X-Strahlenbiischel hat die zu E, E, senkrecht stehenden Geraden Xp, 


X,, X,, X, mit den Abstanden : Thies : Nox» : N34. Man schlagt um E,, E,, E, mit beliebig 


groBem, aber gleichem Halbmesser Kreisbigen, welche die X-Geraden in A), Aj, Aj schneiden. 
Diesen Halbmesser verandert man so lange, bis der Mittelpunkt A) des Kreises durch die drei 
Punkte A{, Aj, Aj auf der X)-Geraden liegt, d. h. bis die Kurve « aller Punkt Aj die Gerade X, 
im Punkte A, schneidet. Damit hat man die Kurbelpunkte A,, A,, A, und die zu ihnen sym- 
metrisch beziiglich der Geraden X, liegenden Punkte A,, A;, A,. Es geniigen auch hier wieder 
zwei beliebige Pole mit einer gleichen Indexziffer, z. B. die Pole P,; und P;,, die sich als Schnitt- 
punkte der Mittelsenkrechten a,; und e,;, sowie a;, und es, ergeben. Die Geraden hy; und hz, 


stehen wieder senkrecht auf E, E,. An diese tragt man die Winkel 2 Oy, und 2 Os, in dem Sinne 


an, wie E; zu E, und zu E, liegt. Die freien Schenkel schneiden sich im Punkte B,, womit auch 
alle Abmessungen der die Aufgabe erfiillenden Geradschubkurbel in der Getriebestellung 5 ge- 
geben sind. ) 


(Eingegangen am 14. November 1951.) 
Anschrift des Verfassers: Ing. Kurt Hain, Braunschweig, JuliusstraBe 31 f. 
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Uber Hauptschwingungen mit endlichen Schwingweiten*. 
Von Th, Péschl, 


1, Erklérung des Begriffes der Hauptschwingungen fiir endliche Schwingweiten, Die Haupt- 
schwingungen eines Systems von Punkten oder starren Kérpern kénnen anschaulich durch die 
Bedingung gekennzeichnet werden, daB alle zu dem System gehirigen Punkte zu gleichen Zeiten 
durch ihre Gleichgewichtslagen hindurchgehen und auch zu gleichen Zeiten ihre gréBten Aus- 
weichungen aus den Gleichgewichtslagen erreichen. Auf diesem Gedanken beruhen die elemen- 
taren Methoden der Berechnung der Eigenfrequenzen fiir die kleinen Schwingungen der mate- 
riellen Systeme. Die gleiche Definition ist auch anwendbar, wenn die Schwingungen nicht auf 
kleine Ausweichungen von den Gleichgewichtslagen beschrinkt werden. Es ist z. B. auch in 
diesem Fall jeder Lage der einen Masse eines Doppelpendels eine bestimmte Lage — genauer 
gesagt: zwei! — der zweiten Masse durch die Bedingung zugeordnet, da die beiden Massen zu 
gleichen Zeiten durch ihre lotrechten Gleichgewichtslagen hindurchgehen und zu gleichen 
Zeiten ihre gréBten Ausweichungen erreichen. 

Es besteht jedoch eine wesentliche Schwierigkeit, wenn man die Frage zu beantworten hat, 
wie tatsachlich die einander zugeordneten Lagen der beiden Massen des Doppelpendels berechnet 
werden sollen. Es ist wohl bekannt, da®B das Problem des Doppelpendels exakt nur fiir kleine 
Werte der Ausweichungen lésbar ist. Ebenso ist bekannt, da die Berechnung der Schwingdauer 
des einfachen Pendels schon auf elliptische Funktionen fiihrt, und es kann als sicher gelten, 
daf{§ die Erweiterung auf die Berechnung der Schwingdauern des Doppelpendels fiir endliche 
Ausschlage nicht einfacher werden kann. Andererseits kann man einen Naherungswert fiir die 
Schwingdauer des einfachen Pendels angeben, ohne elliptische Funktionen zu verwenden; so 
hat z. B. Lord Rayleigh auf direktem Wege eine zweite Naherung in der wohl bekannten Form 
gegeben: 

2 Diets 1 Sa 

. Poe al eal, (1) 
worin s die endliche Amplitude und / die Lange des einfachen Pendels bedeutet. Fiir kleine 
Schwingweiten verschwindet die Korrektion. 

Es diirfte zweckmaBig sein, das Adjektiv ,,klein“* in der Schwingungslehre nur im Sinn der 
Beschrankung auf lineare Ausdriicke in den Ansatzen und in den Ergebnissen zu verwenden. 

Da die Lésungen der Bewegungsgleichungen des Doppelpendels und ahnlicher Probleme der 
Schwingungstechnik fiir endliche Schwingweiten nicht bekannt sind, wiirden zunachst nur die 
Methoden der Reihenentwicklung iibrig bleiben, die jedoch z. B. im Falle des Doppelpendels 
recht umstandlich werden. Fiir die praktische Auflésung derartiger Probleme kann man jedoch 
einen Kunstgriff anwenden, der im folgenden dargelegt werden soll. 


2, Einfiihrung des Prinzips der kleinsten Wirkung in der Jacobischen Form, Fiir die Integration 
der Bewegungsgleichungen eines Doppelpendels ist es wesentlich, daf in diesen die Zeit t explizite 
nicht vorkommt und daher mittels des Energieintegrals eliminiert werden kann. Durch diese 
Elimination wird die Ordnung des Problems um zwei Einheiten herabgesetzt. Die tatsachliche 
Ausfithrung dieser Elimination erfolgt am einfachsten schon in dem zugehorigen Variations- 
problem, wobei statt des Hamiltonschen Prinzips das Prinzip der kleinsten Wirkung zur An- 
wendung kommt, das in der von Jacobi angegebenen Form bekanntlich so lautet: 


2 
6f 2 (h— V) ds=0, (2) 


wobei V die potentielle Energie und h die Energickonstante bedeutet. Fiir das Doppelpendel 
erhalt man dann unmittelbar eine von t unabhangige Differentialgleichung zweiter Ordnung 
in den beiden Lagenkoordinaten des Doppelpendels (x und y oder p und y), welche gleichen 
Werten der Zeit t entsprechen. Nach der Elimination von ¢ sind diese beiden Koordinaten 


* Erweiterte Bearbeitung eines beim IV. KongreB der Italienischen Mathematiker-Vereinigung in 
Taormina gebaltenen Vortrages. 
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Funktionen voneinander und entsprechen einander in dem oben erklarten Sinn. Jedes Werte- 
paar (x und y oder gy und y) erfillt daher vollstandig die notwendigen Bedingungen fiir die 
Lisung des Problems. Wenn einander zugeordnete Werte dieser Koordinaten den Geschwindig- 
keiten Null entsprechen sollen, also die Grenzlagen der Schwingungen bedeuten, so kann die 
Energiekonstante h sogleich angegeben und die Integration bis zu diesen Grenzlagen ausgefihrt 
werden. Fiir die beiden Hauptschwingungen eines Doppelpendels bilden diese Grenzlagen zwei 
Kurven, von denen jede durch eine Differentialgleichung erster Ordnung bestimmt ist. Diese 
Kurven werden im folgenden kurz als Grenzkurven der Schwingungen bezeichnet und in 
einem Schaubild dargestellt. 

3. Anwendung auf kleine Schwingungen. Zur Erlauterung der angegebenen Methode mégen 
hier zunichst die Bewegungsgleichungen fiir die kleinen Schwingungen des Doppelpendels an- 
gefiihrt werden; der Einfachheit halber mégen die beiden Pendelfaden gleiche Langen (/) und die 
Pendel gleiche Massen (m) haben. Werden als Koordinaten nach Abb. 1 die Abstande von der 
Lotrechten x, y eingefiihrt, so lauten diese Bewegungsgleichungen bekanntlich: 


(3) 


Die beiden letzten Glieder stellen die Kopplungsglieder dar. Fir die Integration dieses Systems 
hat man bekanntlich zu setzen 


x= Acos (wt-+ a), y = Boos (wt-+ a) 
und erhalt fiir @m die algebraische 
Gleichung vierten Grades 


Thre Wurzeln sind 


O12—= (2 - 12) = ? (5) 


und die Verhaltnisse zwischen den 
Amplituden in den beiden Schwin- 
gungsarten haben die Werte 


B = 
q=g=l2+1, | 


(6) 


g =I Co = vi = (7 z ee 
Abb. 1. Grundschwingung und Abb. 2. Schaubild der Schwingungs- 8 
(erste) Oberschwingung des ausschlage fiir das linearisierte Di eid H j Sees 
Doppelpendels. Problem. 1€ eiden auptsc wingungen 


unterschieden sich nicht nur durch 
die Werte der Frequenzen (die Eigenwerte des Problems), sondern auch durch die Werte der 
Verhaltnisse der Amplituden in den beiden Schwingungen. 

In iibersichtlicher Weise kénnen diese Ergebnisse durch ein Schaubild nach Abb. 2 dar- 
gestellt werden, wobei entsprechende Werte der Veranderlichen x und y als Cartesische Koor- 
dinaten auf zwei rechtwinkligen Achsen aufgetragen werden. Die beiden geraden Linien mit den 
Gleichungen 


ee meg Sian (7) 
geben dann die Gesamtheit der Werte von x und y fiir entsprechende Lagen der beiden Massen 
des Doppelpendels in den beiden Schwingungstypen an. 


4, Die Differentialgleichung fiir zugeordnete Lagen des Doppelpendels. Diese beiden Geraden 
und die vollstandige Lésung kann man auch unmittelbar aus dem Jacobischen Prinzip entnehmen, 
das im vorliegenden Falle die Form annimmt 


s[ 20 d= [yh le xt — tar}! ylty2dx—0. (8) 
1 1 
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Wenn der Integrand mit F bezeichnet wird, so lautet die Lagrangesche Gleichung dieses Varia- 
tionsproblems 


dfar\ are 
de\ay’} — ay — ) 
und nimmt im vorliegenden Falle die Form an 


Ble—f(e—eyt ty \fly— £0a2—yyte—naty9=0. (19 


Von dieser Gleichung bendtigen wir nur die beiden Lisungen fiir den singularen Punkt x — 0, 


y= 0 mit den Anfangspunkten x0 und y—0 und den Umkehrwerten x—6, y=n. Die 


Energiekonstante h nimmt dann die Form an 


ae é. ie g2— Ey + : | , (11) 


und der Koeffizient von y’’ verschwindet. Somit bleibt dann die folgende Gleichung iibrig, in 
der wir wieder x, y an Stelle von &, ” schreiben: 


JS are erage (12) 
Diese Gleichung hat die Liésung 
y=t=c; (13) 
daher ist 
ne eS a 2 
Die Wurzeln sind 
a= leay2, (15) 
und dies sind dieselben Werte wie in Ziff. 2. 
Die Schwingdauer ergibt sich aus der Energiegleichung in der Form 
s (€, n) 
J» V2(b— V) J Y2[h—(3 x#*—2 xy + y*)] g/l 
0 0 


indie man y=—c,x bzw. y= c,% einzusetzen hat. Man erhalt daher 


r=42 42) es (17) 
und 


2 
2 4 o0* 


OO) Org a (2 == /2) 4. , (18) 


also wieder dieselben Werte wie in Ziff. 2. 


5. Die Differentialgleichung fiir allgemeine Form der potentiellen Energie, Fiir einen Schwinger 
mit derselben Form der kinetischen Energie T, aber von einer beliebigen Form der poten- 
tiellen Energie V — V(x, y) ergeben sich die Differentialgleichungen der Grenzkurven durch die 
Gleichung 

OV OV 
as een = i()h 19 
Ox oy ‘ (19) 
Die Gleichgewichtslagen x —0, y =0 des zweigliedrigen Schwingers stellen in allen Fallen einen 
singularen Punkt dar, und die Werte der Ableitungen y’ —dy/dx in diesem Punkte kénnen als 
bekannt angesehen werden, da die Lésungen fiir kleine Schwingungen eben durch die elementaren 
Methoden gegeben sind. 


6. Ausfiihrungen iitber das Doppelpendel. Fiir endliche Ausweichungen des Doppelpendels sind 
die in Ziff. 3 verwendeten Koordinaten x, y, also die Entfernungen von der Lotrechten durch 
den Aufhangepunkt, nicht geeignet. Es ist vorteilhafter, die Winkel p, p der beiden Pendelfaden 
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als Lagrangesche generalisierte Koordinaten einzufiihren. Die Ausdriicke fiir die kinetische und 
die potentielle Energie nehmen dann fiir m= 1 die Form an (Abb. 3) 


T = + Ply? + 2 cos (p—9) pp + 29"), V=—gl(cospt2cosg). (20) 


Das Energieintegral lautet T+ _V—=h. Nehmen wir als Anfangs-Koordinaten und Geschwindig- 
keiten, die sich auf die Punkte der Grenzkurve bezichen, die folgenden an: yp =u, p=, P=0; 
y —0, so nimmt die Konstante h den Wert an 


h=— gl(cosv-+ 2 cos u). (21) 


Nach der Elimination von t mittels des Prinzips der kleinsten Wirkung in der Jacobischen Form 
erhalt man die Bedingungen zwischen den Winkeln y und y in der Form des folgenden Variations- 
problems, in der y als unabhangig Veranderliche angesehen und zur Abkiirzung dy|/dp =y" 
gesetzt wird: 


2 2 


f) [pera = ofy/2 aa (cos p+ 2 cos y)| Jy? + 2 cos(y— g) yp’ + 2dp=0. (22) 
i 1 | 

Nunmehr hat auch der Ausdruck fiir ds eine etwas allgemeinere Form angenommen. Wir 
wollen daher direkt die Bedingungen angeben, die fiir die Grenzkurven bestehen, indem wir von 
T und V folgende verallgemeinerte Formen annehmen: 


poe (aw2+2biv+ ci), V=V(u,») ; (23) 


die a, b, c und V sind jetzt beliebige Funktionen von u und v. Fir 
die Punkte der Grenzkurven verschwinden nach Bildung der Euler- 
schen Gleichung die Glieder mit y2 (h — V), und mit der Bezeichnung 
v' = dv/du bleibt die Gleichung iibrig 


OV OV 
etn tain ped. 
Ou Ov 
, 
d (singo+ sin) > aV aV © (24) 
a= —b— 
I Ou Ov 
Abb. 3. Einfiihrung yon » und , é 
als generalisierte Koordinaten.. die unter der Annahme der verallgemeinerten Formen von T und V 


die Differentialgleichung fiir die Grenzkurven darstellt. Als Anfangs- 
werte von v’ sind wieder jene anzufiihren, die sich bei kleinen Schwingungen ergeben und wie 
zuvor als bekannt angesehen werden kénnen. 
Fir die Probleme der Eigenschwingungen des Doppelpendels fiir groBe Ausweichungen, fiir 
welche die oben angefiihrten Formen der kinetischen und potentiellen Energie T und V Geltung 
haben, erscheint v’ in der folgenden Form: 


sin (v —u)c 
oe (v — u) cos u 


(25) 


f Les 
sin u —-> sin v cos (v — u) 


Die angenaherte Integration dieser Gleichung fiir die bekannten Anfangswerte von v’ in den 
Ruhelagen u=0, v=0 kann nach einer der bekannten numerischen Methoden leicht durch- 
gefiihrt werden und liefert in Abb.4 als Ergebnis der beiden mit k, und k, bezeichneten Kurven, 
welche die zusammengehérigen Werte der Koordinaten der Grenzkurven des Doppelpendels 
angeben. In Abb. 5 und 6 sind einige Zwischenlagen des Doppelpendels fiir die beiden Eigen- 
schwingungen fitr die Werte der Winkel fiir die kleinere Frequenz von py =0° bis etwa gp = 70° 
und fiir die gréBere von py =0° bis p=120° angegeben. Diese Ergebnisse werden auch durch 
Versuche bestatigt, bei denen nur zu beachten ist, daB die Pendelfaden fiir gréBere Werte der 
Ausschliage schlaff werden kénnen, da sie bei solchen nicht nur Zug-, sondern auch Druckkriafte 
erfahren. 

Zur vollstandigen Lésung der hier behandelten Aufgabe gehért noch die Bestimmung der 
Schwingdauern fiir jede einzelne dieser Eigenschwingungen, d.h. fiir jedes Paar zugeordneter 
Lagen der beiden Punktmassen des Doppelpendels. Diese Bestimmung kann auch nur durch 
angendherte Ermittlung der auftretenden Integrale geschehen. 
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Da die Schwingdauer des einfachen Punktpendels 


Pe 8 x /- 
& 
und die der Grundschwingung des Doppelpendels fiir kleine Ausschlage 
= Laren 
T= 2a(1+772))/2 (26) 
ist, so schwingt das Doppelpendel langsamer als ein einfaches mit gleicher Pendellange /, und 
seine Schwingdauer verlangsamt sich weiter mit gréBer werdenden Ausschlagen. 


Z 
g 
Abb. 4. Grenzkurven y = y(g) fir das Abb. 5 und 6. Grenzlagen fiir die Grundschwingung 
Doppelpendel. und die Oberschwingung des Doppelpendels. 


7, Ansatz fiir das Dreifach-Punktpendel, Fiir mehrgliedrige Schwinger werden die Ansatze 
naturgema} verwickelter, bleiben aber grundsatzlich ausfiihrbar. Wir wollen uns hier auf das 
Dreifach-Punktpendel beschranken und nur die Ausgangspunkte der Rechnung angeben. Die 
wirkliche Durchfithrung kann daraus unschwer entnommen werden. 

Fir gleiche Massen der Punkte und gleiche Langen /lauten die Bewegungsgleichungen in den 
Lagekoordinaten x, y, z (Abb. 7) 


eo meies 


227229 13=6,29 
Z y 


y¥=—(—2a+3y a 
Be (che tcp | 


Die Eigenfrequenzen (der kleinen Schwingungen) er- 
geben sich als die Wurzeln der Fre quenzdeterminanten- 


gleichung (mit w?l/g = p) 


p3)- 2, 0 | 
2s p—3, 1 = p> — 3 p*+ 18p—6=0; 
0 1, p—1l a ree : 
Abb, 7. Kleine Schwingungen fiir das 
ihre Wurzeln haben (angenahert) die Werte bneatech Bankt peu acl. 


PaO a Wp 2 20 pe= 6,29; (28) 
und die Ausschlagverhaltnisse fiir die drei Hauptschwingungen sind 
| 41: Bee Cabs ape? sro 922, 
ee Carl 5 1,045, | (29) 
| As: B,:C, = 1: — 0,645 : 0,122. 
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Fiir die potentielle Energie erhalt man 


oe (30) 


bo| & 


id i ce eae 
Vi=e Te 2xy+sy ¥ 2 


und die kinetische Energie ist 


THF (PtHP MSH. GD 


Das zugehérige Variationsproblem in der Jacobischen Form lautet, wenn x als unabhangige 
Veranderliche angenommen und dy/dx = y’, dz/dx =z gesetzt wird, 


bf /2b—V) ds = 5f FRY) fy EFPLI de =O. (32) 


Die Grenzkurven sind wieder als die Eulerschen Gleichungen dieses Variationsproblems gegeben 
und sind zwei lineare Differentialgleichungen erster Ordnung in y’ und z’, die nur die Lésungen 
y’ = konst und z’ = konst haben. 

Fiir endliche Schwingungen ist es wieder zweckmabig, statt der Koordinaten x, y, z die drei 
Winkel y, y, y einzufiihren; die Ansitze und die Eulerschen Gleichungen werden dadurch er- 
heblich verwickelter, so daB auf die weitere Ausfiihrung verzichtet werden kann, zumal der 
Rechnungsvorgang sonst ganz so verlauft wie zuvor. 


8. Zusammenfassung und Folgerungen, Die Ordnung des Systems der Bewegungsgleichungen 
eines Systems mit mehreren Freiheitsgraden kann durch Elimination der Zeit mittels des Prinzips 
der kleinsten Wirkung — zweckmafig in der von Jacobi angegebenen Form — um zwei Ein- 
heiten erniedrigt werden. Fiir das Doppelpendel ergibt sich dadurch eine einzige Differential- 
gleichung zweiter Ordnung. Obwohl ihr Integral nicht bekannt ist, so kénnen doch die zuge- 
ordneten Lagen der beiden schwingenden Massen auch fiir beliebig groBe Ausschlage ermittelt 
werden; und zwar sind diese fiir die beiden Scharen der Hauptschwingungen durch je eine Diffe- 
rentialgleichung erster Ordnung gegeben, die angenahert durch eine der bekannten Naherungs- 
methoden gelést werden kénnen. 

Dieselbe Methode ist unmittelbar auch fiir das Doppel-Kérperpendel und — nach sinngemaBer 
Erweiterung — auch fiir mehrfache Pendel anwendbar und selbstverstandlich auch fiir alle 
schwingungsfahigen Systeme verwandter Art. 

Von besonderem Interesse diirfte jedoch die Tatsache sein, daB durch die angegebene Eli- 
mination der Zeit t auch das allgemeine Problem des Doppelpendels (fiir beliebige Anfangs- 
bedingungen), dessen Lisung noch nicht bekannt ist, — und dhnlicher Probleme — auf die 
Integration einer einzigen Differentialgleichung zweiter Ordnung zuriickgefiihrt und dadurch 
wesentlich vereinfacht wird. 

Die Behandlung der erzwungenen Schwingungen fiir Systeme mit endlichen Ausschlagen, 
wie sie z. B. bei dem System Glocke und Kléppel verwirklicht ist, erfordert jedoch eine besondere 
Untersuchung. 


(Eingegangen am 18. Dezember 1951.) 
Anschrift des Verfassers: Prof. Dr. Theodor Péschl, (17a) Karlsruhe i. B., Wendstr. 5, Fernruf 5414. 
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Schubspannungsverteilung und Energiedissipation 
bei turbulenten Grenzschichten. 


Von J. Rotta, 


1, Einleitung. Die Berechnung turbulenter Grenzschichten wird bei den meisten bjsher be- 
kannt gewordenenen Naherungsverfahren auf die Lésung von zwei simultanen, gewéhnlichen 
Differentialgleichungen zuriickgefithrt. Dieses sind die v. Kdrmdnsche Impulsgleichung und eine 
Differentialgleichung fiir die Anderung, die die Form der Grenzschichtprofile beim Fortschreiten 
in Strémungsrichtung erfahrt. 

Damit das Gleichungssystem lésbar ist, miissen noch weitere Zusammenhange gegeben sein. 
So muB vor allem, wenn als Gleichung fiir die Anderung der Profilform der Energiesatz eingefiihrt 
wird, die Abhangigkeit zwischen der auftretenden Dissipationsfunktion — das ist das iiber die 
Grenzschichtdicke erstreckte Integral der durch die Zahigkeit dissipierten Energie — und der 
Form des Geschwindigkeitsprofiles bekannt sein. Die Ermittlung dieses Zusammenhanges er- 
scheint zur Zeit als eine der dringendsten Aufgaben auf dem Gebiete der turbulenten Grenz- 
schichten. 

In einer unlangst erschienenen Arbeit! hat sich der Verfasser bemiiht, diese Abhangigkeit 
aus Versuchsergebnissen rein empirisch festzustellen. Da die Bestimmung des Dissipations- 
integrals aus Versuchsergebnissen von Natur aus mit gewissen Ungenauigkeiten behaftet ist, soll 
in der vorliegenden Arbeit der Versuch unternommen werden, den genannten Zusammenhang 
theoretisch zu klaren und durch erneute Sichtung vorliegenden Versuchsmaterials vermittels 
einer neuen Auswertungsmethode auf eine festere Grundlage zu stellen. 

Die Ausfiithrungen beschrénken sich auf stationadre, zweidimensionale Grenzschichten in 
inkompressiblen Medien. 


2. Uber die Einparametrigkeit der Grenzschichtprofile. Als auffallendes Ergebnis experimen- 
teller Untersuchungen an turbulenten Grenzschichten wurde gefunden, daB sich die Geschwin- 
digkeitsprofile, d. h. die Verteilung der Grundstrémungsgeschwindigkeit iiber den Wandabstand, 
mit guter Annaherung als einparametrige Familie darstellen lassen. Diese Kigenschaft wird bei 
den meisten Naherungsverfahren in Ermangelung besserer Anhaltspunkte als Arbeitshypothese 
eingefiihrt. 

Nun beeinflussen bei hinreichend groBen Reynoldszahlen die kinematische Zahigkeit y und 
die geometrische Beschaffenheit der Wand (Wandrauhigkeit) die Geschwindigkeitsverteilung nur 
in einer im Vergleich zur Grenzschichtdicke 6 sehr diinnen Schicht 0,,in Wandnahe; im iibrigen 
Bereich der Grenzschicht erscheinen die StrémungsgréBen dagegen bei geeigneter Normiecrung 
als von der Zahigkeit und Wandrauhigkeit unabhangig. Die Beachtung dieser Erkenntnis fiihrt 
dazu, nicht die ganzen Geschwindigkeitsprofile als einparametrig anzusehen, sondern nur den 
duBeren, von der Zahigkeit und Wandrauhigkeit unbeeinfluBten Teil als einparametrige Schar 
aufzufassen. Es ist dann mdglich, den Einflu8 der Reynolds-Zah] und der Wandrauhigkeit 
naherungsweise zu beriicksichtigen, ohne daB ein grundsatzlich komplizierteres Berechnungsver- 
fahren notwendig ware. Diese Uberlegungen wurden in der zitierten Arbeit des Verfassers 
behandelt. 

Der benétigte Zusammenhang zwischen dem erwahnten Dissipationsintegral und dem Form- 
parameter der Geschwindigkeitsprofile kann aber nur dann als eindeutig vorausgesetzt werden, 
wenn zu jedem Geschwindigkeitsprofil auch eine ganz bestimmte Verteilung der statistischen 
SchwankungsgréBen gehort. Mit anderen Worten, es muf vorausgesetzt werden, dafs sowohl die 
Geschwindigkeitsprofile als auch die Verteilungen der Geschwindigkeitsschwankungen und Korre- 
lationen zwischen ihnen als einparametrige Kurvenscharen angesehen werden kénnen. Diese 
Voraussetzung wird durch vorliegende Versuchsbeobachtungen gestiitzt, doch lassen sich die 
Bedingungen fiir ihre Giltigkeit auch theoretisch angeben. 

AuBer der Zulassigkeit der iiblichen Grenzschichtvereinfachungen ist die erste Bedingung fiir 
den eindeutigen Zusammenhang zwischen den Geschwindigkeitsprofilen und der Verteilung der 


1 J. Rotta, Uber die Theorie der turbulenten Grenzschichten. Mitteilungen aus dem Max-Planck-Institut 
fiir Strémungsforschung, Nr.1. Géttingen 1950; gekiirzte Darstellung Ing.-Arch. 19 (1951), S. 31. 
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Geschwindigkeitsschwankungen eine hinreichend grofBe Reynoldszahl, denn nur dann bleibt der 
Bereich, in welchem die Str6émungsverhaltnisse unmittelbar von der Zahigkeit beeinfluBt werden, 
klein gegen die Grenzschichtdicke. Die weitere Beurteilung, ob zu einer gewissen Verteilung der 
Grundstrémungsgeschwindigkeit U iiber dem Wandabstand y auch eindeutig bestimmte Ver- 
teilungen der Geschwindigkeitsschwankungen iiber y gehéren, erfordert die Beantwortung der 
Frage, wie weit der Turbulenzmechanismus nur durch die Ableitungen von U nach y gesteuert 
wird und wie weit auch die Verteilung von U iiber x (der Koordinate in Strémungsrichtung) und 
die senkrecht zur Wand gerichtete Geschwindigkeitskomponente V daran beteiligt sind. 

Zur Untersuchung des Turbulenzmechanismus kann man sich eine Anzahl von Gleichungen 
herleiten, indem man die Navier-Stokesschen Differentialgleichungen fiir die allgemeine, zeitlich 
und riumlich schwankende Fliissigkeitsbewegung mit einer oder mehreren Komponenten der 
Geschwindigkeit multipliziert und dann eine statistische Mittelung vornimmt. Solche Ent- 
wicklungen gehen bereits auf 0. Reynolds zuriick und wurden fiir Produkte von drei Geschwin- 
digkeitskomponenten von P.Y.Chou! angeschrieben. Ohne hier darauf einzugehen, welches 
physikalische Bild auf diese Weise von dem Turbulenzmechanismus entsteht, und welche Méglich- 
keiten sich zu seiner quantitativen Erfassung bieten”, midge folgendes dargelegt werden. 

Durch das Auftreten der substantiellen Ableitung der Geschwindigkeit in den Navier-Stokes- 
schen Gleichungen von der Form dp/dt+-» vy, wobei » der Geschwindigkeitsvektor ist, erscheinen 
in den genannten Gleichungen fiir den Turbulenzmechanismus Glieder, die einen konvektiven 
Einflu8 der Grundstrémung auf die Turbulenzbewegung beschreiben. Bei ausgebildeten, turbu- 
lenten Strémungen durch Leitungen gleichbleibenden Querschnitts fallen diese konvektiven 
Glieder ohne weiteres fort, nicht aber bei Grenzschichten, wenn auch die iiblichen Vereinfachungen 
eingefiithrt werden. Man hat daher von diesen Gliedern eine Einwirkung der Verteilung von U 
in x-Richtung und der V-Komponente auf die Turbulenzbewegung zu erwarten. 

Wir wollen uns die Gré®enverhaltnisse an Hand der Gleichung fiir die kinetische Energie 
der Schwankungen?, 
aU 


0Q 
iyi ue ae ays (2.1) 


dy 


aE aE\ 
ON oe aay le 


konvektives Glied 


klar machen, wobei E die kinetische Turbulenzenergie 
pst hee 
O 


(2,2) 


S die Dissipation 
___{ fau\2 du\2 | dv\2 | Ow\2 , {dw , dv\2 _ (du, Ow\? , /dv , du\? 2 
S=ol(ay) + 2(ze) +2(59) +2() +(e +a) tet ae) t (eta) 9 
Q die durch turbulente Vermischung in y-Richtung transportierte Energie (Energiediffusion) 


= —o(F+ F)+ (ASS + aye), (2.4) 


dy : 
t die Schubspannung, p die Schwankungen des Druckes und @ die Dichte des Strémungsmediums 
sind. Die Energie E, die Dissipation S und die Energiediffusion Q werden auf die Masseneinheit 
der Fliissigkeit bezogen. 


Nun hat die Komponente V bekanntlich die GréBenordnung U, dé/dx, ebenso ist 
OE/Ox ~ (0E/dy) (dd/dx) , 


wenn U, die Gréfe von U am AuSenrand der Grenzschicht ist. Hieraus folgt, daB UdE/dx 
und VoE/dy etwa von gleicher Gri®enordnung sind; das konvektive Glied in (2,1) hat insgesamt 
die GréSenordnung 
dE dé 
ayadae 

1 P. Y. Chou, Quarterly of Appl. Math. 3 (1945), S. 38. 

* Vel. hierzu J. Rotta, Z. Physik 129 (1951), $. 547; 131 (1951), S. 51. 
_ ® Gl. (3,4) der in FuBnote 1 von S.195 zitierten Arbeit. Wie friiher bezeichnen die groBen Buchstaben 
die zeitlich gemittelten Geschwindigkeiten, die kleinen Buchstaben bezeichnen dagegen die Geschwindig- 


keitsschwankungen, deren Mittelwerte gleich Null sind. Uberstreichungen der Schwankungsgréfen be- 
deuten zeitliche Mittelung. 


wo 
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Aus dem Vergleich dieses Wertes mit dem Diffusionsglied 
aQ Of. Wet w 
Q dy == dy 5 | oe ee 
folgt, daB das konvektive Glied in (2,1) nur dann klein gegen das Diffusionsglied (oder gegen 
eines der anderen Glieder in (2,1)) und somit vernachlassighar ist, wenn die Bedingung 


di v 
ib <i @.5 


iC 
erfiillt ist; denn es ist gréBenordnungsmabig 
0 ( a) 
v 
oy 


Die GréBe )v/ U, ist aber erfahrungsgem4B selbst klein gegen 1. 

AuBer daB die Reynoldszahl hinreichend gro8 ist, muB also weiter die Anderung der Grenz- 
schichtdicke in Strémungsrichtung dd/dx von zweiter Ordnung klein gegen 1 sein, wahrend 
fiir die Zulassigkeit der iiblichen Grenzschichtvereinfachungen nur erforderlich ist, da® déd/dx 
von erster Ordnung klein gegen 1 ist. Dies soll heiBen, um einen eindeutigen Zusammenhang 
zwischen dem Geschwindigkeitsprofil und der Verteilung der Schwankungsgeschwindigkeiten 
annehmen zu diirfen, muB dé/dx etwa eine Zehnerpotenz kleiner sein, als es fiir die Annahme 
der von L. Prandtl eingefiihrten Vereinfachungen bei gleichen Genauigkeitsanspriichen not- 
wendig ist. Fiir die praktische Grenzschichtrechnung vermittels Impuls- und Energiesatz sind 
die Anforderungen allerdings nicht ganz so scharf. Anstatt einer strengen Einparametrigkeit der 
Geschwindigkeitsprofile und Schwankungsverteilungen ist nur ein eindeutiger Zusammenhang 
zwischen dem Dissipationsintegral und dem Formparameter fiir das Geschwindigkeitsprofil 
notwendig. Beides sind integrierte GréBen, bei denen sich maBige Abweichungen in den nor- 
mierten Verteilungen noch nicht merkbar auswirken, so da diese Voraussetzung viel leichter 
erfiillbar ist. 

Bei den in der Regel vorliegenden Verhaltnissen bei turbulenten Grenzschichten ist dd/dx 
von der GréBGenordnung 0,01 bis 0,02, so da damit der experimentell beobachtete eindeutige 
Zusammenhang zwischen Schwankungs- und Grundgeschwindigkeitsverteilungen erklarlich ist. 
In der Nahe der turbulenten Ablésestelle ist dd/dx allerdings keineswegs von héherer Ordnung 
klein gegen 1, so daB hier bei Anwendung der einparametrigen Berechnungsverfahren keine Treff- 
sicherheit mehr gewahrleistet werden kann. 


3. Naherungsansatz fiir die Schubspannungsverteilung. Die je Zeit- und Langeneinheit dissi- 
pierte Energie steht in engem Zusammenhang mit der Schubspannungsverteilung. Die genaue 
Kenntnis dieser Verteilung wiirde zweifellos fiir die Berechnung turbulenter Grenzschichten einen 
bedeutenden Fortschritt bedeuten, und es ist daher schon mehrfach versucht worden, sie aus 
Versuchen oder theoretisch zu ermitteln. 

K. Fediaevsky macht in seiner Arbeit! einen Polynomansatz und fordert insbesondere, daB 
die Schubspannungsverteilung die aus der Grenzschichtgleichung und Kontinuitatsgleichung 
folgenden Wandbedingungen 

‘ or Op Or 

yo By Bx” yt oe 

erfiillt. Nun ist die Schubspannung aufer in der meistens sehr diinnen Schicht in Wandnahe, 

in der die Strémungsverhaltnisse unter dem vorherrschenden Kinflu®B der Zahigkeit stehen, 
durch die Beziehung 

T=— ow (3,2) 

gegeben. Sehen wir die im vorigen Abschnitt untersuchte Annahme, nach welcher zu jedem 

Geschwindigkeitsprofil eine bestimmte Verteilung der statistischen Schwankungen gehdrt, als 

giiltig an, so ist jedem Geschwindigkeitsprofil ein bestimmtes wv-Profil und mithin auch ein 

bestimmtes Schubspannungsprofil zuzuordnen. Diese Annahme ist aber nicht mit den Wand- 

bedingungen (3,1) vertraglich, auBer wir wiirden jedem drtlichen Druckgradienten ein bestimmtes 

Geschwindigkeitsprofil zuordnen. Letztere Methode, die dem Berechnungsverfahren von A, Buri? 


1 K, Fediaevsky, J. Aeron, Sciences 4 (1937), 8. 491. . 
2 4, Buri, Eine Berechnungsgrundlage fiir die turbulente Grenzschicht bei beschleunigter und ver- 
zogerter Grundstrémung. Dissertation. Ziirich 1931, 
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zugrunde liegt, hat sich jedoch nur beschrankt als brauchbar erwiesen. Im Abschnitt7 werden wir 
auf diese Frage naher eingehen. 

Wir werden die Bedingungen (3,1) also zunaichst auBer acht lassen, ohne ihre Exaktheit 
irgendwie in Zweifel ziehen zu wollen. Fiir das Geschwindigkeitsprofil hatte sich der Naherungs- 


ansatz! 
Ce A (1 — n) - Inn fiir = EMS. | (3,3) 


v* “ 


als brauchbar gezeigt. Dabei ist 

ee * ? (3,4) 
v*—Yz,/0 die mit der Wandschubspannung 7, gebildete Schubspannungsgeschwindigkeit, x eine 
universelle Konstante (% —0,4), A der Profilparameter und 6, die Dicke der von Zahigkeit und 
Wandrauhigkeit beeinfluBten Unterschicht. Dieser Ansatz stellt das Geschwindigkeitsprofil im 
Mittel ganz gut dar, ohne daB er mit der sich aus den Grenzschichtgleichungen ergebenden Reihen- 
entwicklung von der Wand aus genau itbereinstimmt. Bei der Suche nach einem Ansatz fiir die 
zugehirige Schubspannungsverteilung scheint deshalb die AuSerachtlassung der Wandbedin- 
gungen berechtigt zu sein. 

Wir machen wie Fediaevsky einen Potenzreihenansatz von der Form 


% WO oe Cw 

Se a [Oy 9 i O30? “eis ete = (3,5) 
nehmen aber statt der Wandbedingungen (3,1) einen linearen Zusammenhang zwischen a, und A 
i Pea fh. (3,6) 


Der lineare Zusammenhang zwischen a, und A ist theoretisch begriindet. Wollte man z. B. mit 
einem Differentialgleichungssystem fiir den Turbulenzmechanismus, wie es in den zitierten 
Arbeiten? entwickelt wurde, unter Vernachlassigung der konvektiven Glieder die zu der Ge- 
schwindigkeitsverteilung (3,3) gehérige Schubspannungsverteilung berechnen, so kame man 
tatsdchlich auf einen linearen Zusammenhang zwischen a, und A. Eine einfache quantitative 
Abschatzung lat sich mit dem Prandtlschen Mischungswegansatz 
er eee 
T= Ol PEen (3,7) 


machen. Verwendet man dabei fiir den Mischungsweg | den Ahnlichkeitsansatz nach v. Kérmdn, 


dU/dy 3.8) 


[= 4 |= 
*|\PU/dy? 


) 


so bekommt man a,=— 4 A, also a= 4. 

In (3,6) ist Ay diejenige GréBe von A, die fiir die Darstellung des Profiles fiir die turbulente 
Grenzschicht ohne Druckgefalle zu wahlen ist (A4,~0,9). Hierdurch wird erreicht, daB® bei 
dp/dx —0 sich a, 0 ergibt, wie es der Theorie und den Beobachtungen entspricht. Abgesehen 
von der Abschatzung mit dem Mischungswegansatz soll die Proportionalitatskonstante « auf em- 
pirische Weise bestimmt werden. 

Als weitere Annahme zur Bestimmung der Koeffizienten in (3,5) wahlen wir fiir 

rae hee o.! Op 
N= Vets ees (3,9) 
Diese Bedingungen begriinden sich damit, daB die quadratischen Mittelwerte der Geschwindig- 


keitsschwankungen, wu? usw., und damit auch wv am Auenrand der Grenzschicht asymptotisch 
auf Null abfallen. Versuchsbeobachtungen? zeigen weiter, da® auch der Korrelationskoeffizient 


wo/| ue vzam Aufenrand der Grenzschicht auf Null abfallt, so daB der Mittelwert des Produktes uv 


am AuBenrand von héherer Ordnung gegen Null geht als die Schwankungsquadrate u? und v2. 
Hiernach ware es berechtigt, fiir 71 auch noch die Bedingung 
te 0 
dy? aan (3,10) 
1 Gl. (6,14) der FuBnote 1 von S$. 195 zitierten Arbeit. 
2 FuBnote 2, von S. 196. 


° G. B.Schubauer und P.S. Klebanoff, Investigation of Separation of the Turbulent B 
NACA T.N. 2133 (1950). i N ventilate oa 
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vorzuschreiben. Diese Vorschrift scheint allerdings etwas zu streng zu sein, weil der Ansatz (3,3) 
wegen des scharfkantigen Uberganges in die freie Strémung eine etwas zu kleine Grenzschicht- 
dicke 6, bezogen auf die Verdrangungsdicke 


(0-2. (3,11) 
0 


oder auf die Impulsverlustdicke 


f U U 
3, le (1 u,) ay (3,12) 
0 


ergibt. Vergleichsweise soll deshalb die Schubspannungsverteilung ohne und mit Bedingung (3,10) 
untersucht werden, wobei der Ansatz (3,5) im ersten Fall ein Polynom dritten Grades ist, wahrend 
im zweiten Fall ein Polynom vierten 
Grades angesetzt wird. Die Koeffi- 


zienten 4), @3, a, lassen sich dann in @ 
Abhangigkeit von a, ausrechnen: 
dr e j 98 5S 
(1) aga 0 fir gl; \bersuche 
a= — (3 —=-2 a,), Ob cS 
a,—2-+ a,, 
== (( 
@ ve (3,13) gy 
t a: : 
(2) i OFir i: 
a, = — (6-++ 3 a,), 2 
a,=—8+3a, yo™ 
a, = — (3+ a) =- a” 
; ; 0 GOS G7 GIS G2 O25 O38 
In Abb. 1 sind die Schubspan- Abb. 1. Schubspannungsverteilung bei konstantem Druck. 


Rechnungsannahme; a, = 0; A = 0,9. Versuche nach F’. Schultz-Grunow. 


nungsverteilungen T/T, fiir die Grenz- 
schicht mit konstantem Druck 4z 
To 

5 


(a, = 0, A) = 0,9) 


itber der dimensionslosen Abszisse 
yv*/6,U, dargestellt und mit der 
aus Versuchen von F’. Schultz-Gru- 
now! errechneten Verteilung ver- 
glichen. Die Abszissen errechnen 
sich zu? 


. a:My=140:Re=85 107 A=3.0 
+ Bs thy =154: Rea= 35+ 105 A=57 
oar Mag 115: Rep 45° 10°, A=§0 


yur x 
ce aml (3,14) 
i] 34 


Bei positivem a, ergeben sich 
Verteilungen, die qualitativ den 
-Versuchsbeobachtungen ent- 
sprechen, wie aus Abb.2 zu er- 
sehen ist. Bei dieser Auftragung 


' ee a | eZ 3 i 5 6 7 
wurden fiir vorgegebene Werte 2 
' Pal : Abb. 2. Schubspannungsverteilungen bei Druckanstieg. 
Eh, = 6,/0, und Re, = U, 6,/v mit Kurven a,, by, ¢,: Rechnungsannahme (Q) FO 2r oie 
| ; 4 Kurven a,, b,, c,: Rechnungsannahme (2); «4 = 4,53. 
den Beziehungen der erwahnten Versuchspunkte nach G.. B: Schubauer und P. S. Klebanoff. 


Arbeit? zunachst die Werte A 
-ermittelt und dann unter Vorwegnahme der Resultate des Abschnitts 6 mit a—2,57 im ersten 
Fall und «4,53 im zweiten Fall die 1/t)-Verteilungen berechnet und iiber y/0, aufgetragen. 


| 1 F. Schultz-Grunow, Luftf.-Forschg. 17 (1940), S. 239. 
2 Vel. Gl. (6,15) der in FuBnote 1 von S. 195 zitierten Arbeit. 
3 FuBnote 1 von S. 195, Gl. (6,3), (6,9), (6,16) und (6,19). 
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In diese Darstellung wurden die entsprechenden Versuchspunkte nach G. B. Schubauer und 
P.S. Klebanoff} eingetragen, die vermittels einer Hitzdrahtsonde gemessen wurden. 

Aus den Versuchsergebnissen mu8 man mit einiger Sicherheit schlieBen, da® sich die Schub- 
spannungen iiber eine gréRere Schichtdicke erstrecken als die Abweichungen zwischen der Grund- 
stroémungsgeschwindigkeit U und der Potentialstrémungsgeschwindigkeit U,. Diese Tatsache 
lieBe sich durch entsprechende Anderung der Annahme (3,9) bzw. (3,10) mithelos berichtigen. 
Wir wollen die Betrachtungen iiber die Schubspannungsverteilungen hauptsichlich dazu be- 
nutzen, den Zusammenhang des schon erwahnten Dissipationsintegrals mit dem Profilpara- 
meter A zu ergriinden; fiir diesen Zweck wird die vorgeschlagene Naherung als ausreichend 
erachtet. Die Annahmen fiir den ersten Fall passen offensichtlich etwas besser. 


4, Energiedissipation. Aus der Energiebilanz (2,1) ergibt sich durch Integration tiber die 
Grenzschichtdicke 


,_4au i df 

BLN pore ue 4,1 
[rao 9 [Sdy+ez, [UE dy. (4,1) 
‘ a i 


Das Integral von S iiber y, fiir das der Ausdruck D, 
D= {Sdy (4,2) 
0 


eingefithrt wurde, ist die gesuchte Funktion, die wir jetzt unter Benutzung der im Abschnitt 3 
gemachten Abschatzungen fiir die Schubspannungsverteilung berechnen wollen. Setzt man in 
(4,1) die Schubspannungsverteilung nach (3,5) und die Geschwindigkeitsverteilung nach (3,3) 
fiir den auBeren Teil der Grenzschicht 7 > 6,,/6 bzw. das Wandgesetz? 

U * 

tee (4.3) 


v v 


fiir den wandnahen Teil 0 < 7 <1 ein, so 1aBt sich die Quadratur fiir 


oo ul 


"dU te CRU! yy 
| Ga % | att (4,4) 
0 


ausfiithren. Dabei machen wir wieder von der schon in der friitheren Arbeit vorgenommenen Auf- 
teilung der Grenzschicht in einen wandnahen Teil (0 < 7’ <7) und einen auBeren Teil (7’ > 7) 
Gebrauch, wobei 0,,/6 <4 <1 sei. In dem wandnahen Teil kann man die Schubspannung in 
erster Naherung gleich der Wandschubspannung Tt, setzen; auf diese Weise erhalt man 

1 


n 1 
BONUS oo 7 “dU  F GAO 
—= = aA com SAREE Rs, 
af a ae dy % ie dy’ +t aii dy 
0 0 
1 1 a a a 
SRA hd aera ered eee es) (4,5) 
, Afyp 1% 1 1 4s ! oa MAIR CS 
I “ | 2 I 3 i} 4 5 | hé 


In dem Gebiet, in welchem die Integrationsgrenze 7 liegen soll, gilt das universelle Wandgesetz3 
turbulenter Strémungen 


er vty _ off Wa ee 1 
iw (Zim ! c) ot (Fm : bilnn+ ¢), (4,6) 


v 


wobei C eine von der Wandbeschaffenheit abhangige Konstante ist. Wir wollen der Kiirze 
wegen noch 


a a. a a 
Tie 0 Hegiahiog ie gietee aor et Se ea (4,7) 


* FuBnote 3 von S.198. Diese Messungen sind die einzigen ihrer Art, die bekannt geworden sind. Die 
daraus extrapolierten Wandschubspannungen zeigen erhebliche Abweichungen von den Feststellungen 
anderer Autoren, so daB Zweifel an der quantitativen Richtigkeit aufkommen. Der Vergleich (Abb. 2), 
kann daher nur qualitativ bewertet werden. 

2 Gl. (4,7) der in FuBnote 1 von S.195 zitierten Arbeit. 

3 Vel. J. Rotta, Ing.-Archiv 18 (1950), S. 277. 
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einfiihren. Nach Einsetzen von (4,6) folgt dann aus (4,5) mit (3,14) und Re, = U,4,/y 


a 1 
Lone 1 oot er 
[ea erers(ein Re, ——n—2_ 424 c). (4,8) 
0 
Mit (3,13) ist fiir die beiden untersuchten Annahmen 
5 a 1 a 
1 AW ba es ae 
( ) t 6 I 3 | A ie | ae | 
135). ce 2° Ae a7) 
2) =— =+4 4 
Fir das Dissipationsintegral D ergibt sich dann aus (4,1) 
il 
1+—A 7 
sa 1 1 2 I d 
Diezakar ce haps era t 4,10 
: (in Re, epee ee es c) =| WE dys (4,10) 
0 
Frither wurde fiir das Dissipationsintegral die Form? 
— »x3(1 
Daa, & In Re, + 6) (4,11) 


angegeben. Die dimensionslose Differenz G—C ist unter der eingefiithrten Voraussetzung der 
Einparametrigkeit nur eine Funktion des Profilparameters 


ao B CS U, 2 yor 
C= [| aa iu) (4,12) 
0 


und wurde empirisch aus Versuchen bestimmt. Durch Vergleich von (4,11) mit (4,10) bekommen 
wir jetzt fiir G die Beziehung 


1 oo 
ees 
i! 2 [, il cel 
0 


Wenn man hierbei die GréBe der Turbulenzenergiestromanderung at UE dy als Funktion von A 
0 


ausdriicken kann, dann hat man den gesuchten Zusammenhang zwischen G und A qualitativ 
gefunden. Als Unbekannte tritt nur noch die GréBe a von (3,6) auf. Zwischen den Formpara- 
metern J, und A besteht der Zusammenhang* 


245 4-H 
= i : (4,14) 
rg ( +. 74) 


5. Turbulenzenergiestrom, Zwischen der Anderung des Turbulenzenergiestromes 


und dem Profilparameter A besteht kein eindeutiger Zusammenhang. Hine Abschatzung ist 
fiir die sog. ahnlichen Lésungen méglich. Diese Loésungen wurden vom Verfasser behandelt und 
sind dadurch ausgezeichnet, daB das Geschwindigkeitsprofil langs der Wand nur affin verzerrt 
wird; sie existieren fiir vorgegebene Geschwindigkeitsverteilungen von der Form 


U,=ax™, (5,1) 
wobei a und m konstant sind. Die Dicke der Grenzschicht wachst dabei linear mit der Ent- 
1 Gl. (6,26) der in FuBnote 1 von S. 195 zitierten Arbeit. 


2 Gl. (6,16) der in FuBnote 1 von S. 195 zitierten Arbeit. 
15 


; Baie in kedee ; ; eaicacee 
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(ui) 


2— *Om+1+ mH)’ 


fernung x. Fir die Impulsverlustdicke gilt? 


é (5,2) 


dagegen ist v*/U, von x unabhangig. Fir gegebene Reynoldszahlen Re, stehen m einerseits 

und v*/U,, Hy. und damit A andererseits in eindeutigem Zusammenhang miteinander, der aller- 

dings nicht explizit hingeschrieben werden 

iy raha dy kann. 
0 


10 Zur Abschatzung des Turbulenzenergie- 
stromes wurde die Formel? 
| | ) J UE dy = 0,65 5,U2v* (5,3) 
a 

8 i. angegeben. Eine Verbesserung dieser 
if Beziehung ware mdglich, wenn ahnliche 
/ | Uberlegungen, wie sie zur Ermittlung der 
f is N la | Schubspannungsverteilung angestellt wur- 
e } den, auf die Verteilung der Schwankungs- 
| S 8 / energie angewendet wiirden. Da die Wir- 
i oy Z kung des Turbulenzenergiestromes jedoch 
C ix von nur geringer Bedeutung ist, wollen 
5 7 iv “e wir von einer solchen Verfeinerung ab- 
® sehen und uns weiterhin der Formel (5,3) 
bedienen. Fithrt man darin (5,1) und (5,2) 

| : ein, so ergibt sich 

it} 0) 
VY aoe 
Ha | eed, a (1+3m) Ay 
. Ag | wigg | UE A= 005 gat anti 
: (5,4) 
nn ? 
4 ‘ oe [ In Abb.3 ist dieser Wert fiir verschiedene 
A Re, iiber A aufgetragen. Da dieses Glied, 
Z | wie schon erwahnt, nur geringen Anteil 
(ZB am Energiegleichgewicht und deshalb in 
r (4,13) nur die Bedeutung eines Korrektur- 
liedes hat, wird es fiir den qualitativen 
0 2 4 6 8 70 12 Tee ? 

; Pei ert ‘ Zusammenhang zwischen G und A ohne 
ais Mier raced ai tarry e lia weiteres statthaft sein, die Ergebnisse der 


ahnlichen Liésungen auf den allgemeinen 
Fall beliebigen Druckverlaufes zu iibertragen. Wir wollen es dariiber hinaus sogar als zulassig 
erachten, die Turbulenzenergiestromanderung in (4,13) durch die in Abb.3 gestrichelt ein- 
gezeichnete Kurve 


eld, | UE dy = 0,9-+ 0,04 42 (5,5) 


v*3 dx 
0 


anzunahern. 


6. Auswertung von Versuchsergebnissen, Die Ermittlung des Dissipationsintegrals aus ge- 
messenen Geschwindigkeitsverteilungen in Grenzschichten ist nur iiber eine Differentiation 
moglich. Deshalb diirfte die in der fritheren Arbeit durchgefiihrte Auswertung der Versuchs- 
ergebnisse hinsichtlich der Zuverlassigkeit nicht voll befriedigen. Nachdem wir in den vorher- 
gegangenen Abschnitten den Zusammenhang zwischen dem Dissipationsintegral D und dem 
Geschwindigkeitsprofil bis auf eine empirische Konstante a festgelegt haben, ist es unter Be- 
nutzung der Beziehung (4,13) mit (4,9) und (5,5) jetzt méglich, die Auswertung auf Quadraturen 
zuriickzufiihren, wobei sich die versuchsbedingten Streuungen der Versuchsergebnisse zweifellos 


1 Gl. (7,4) der in FuBnote 1 von S. 195 zitierten Arbeit. 
2 Gl. (6,22) der in FuBnote 1 von S. 195 zitierten Arbeit, 
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nicht so stérend auswirken wie bei dem erstgenannten Verfahren. Man erhialt dabei fiir jede 
MeBreihe einen ,,gemittelten‘‘ a-Wert. 

Bevor wir uns dieser Auswertung zuwenden, wollen wir kurz das Dissipationsintegral bei der 
Platte ohne Druckgefalle mit den Versuchsergebnissen von F’. Schultz-Grunow! vergleichen, die 
noch nach der alten Differentiationsmethode ausgewertet wurden. Mit der Beziehung 


U. 1 Q 
eee Re,+ B, (6,1) pee = 
wobei? 
B= C+ 0,824 | 23 
er tG..) 
eee Ee | a2 | 
In | 
w 
ist, laBt sich D nach (4,11) auf 2 
die Form 
p= v3(%14 fF) 6,3) 
ras v* i 20 
bringen. Hierin ist mit (5,5) 
719 
eae = ate A ; 
«| (6,4) 4 
1 hice aml ce ait 
— 0,9 — 0,04 A?. | Bog 2 2 Pu 25 26 27 28 29 
Bei der Platte ohne Druckgefalle Abb. 4. Dissipationsintegral bei konstantem Druck in Abhangigkeit von U,/v*. 
‘ Kurve J]: Rechnungsannahme (1): F = — 3,736) x = 0,4; 
ist der Wert von F wegen a,—0 Kurve 2: Rechnungsannahme (3; [3 — ANN Jy Al SOO 


von der Konstanten « unabhangig 


und ergibt sich mit 4,—0,9 und x=0,4 bei den unterschiedlichen Annahmen fiir die Schub- 
spannungsverteilung zu 


(1) ee 3,736, | 

(2) F=—4,586.{ 
Die entsprechenden Geraden sind in Abb.4 eingetragen und zeigen befriedigende Ubereinstim- 
mung mit den Versuchsergebnissen. 


Die neuere Auswertung geht auch von dem Energie-Integralsatz fiir Grenzschichten keil- 
férmiger Strémungen aus?: 


(6,5) 


o0 fur ad 
Dans isc RSET gear I dx x + Xp 
0 
Hierbei ist 
SO r U\2 
reece (6 
1 1 


die Energieverlustdicke. Eliminiert man die Entfernung x % der MeBstelle von dem gemein- 
samen Schnittpunkt der Stromlinien mit Hilfe des von A. Kehl* fiir Keilstrémung gegebenen 


Impulssatzes 5 
dos fg fly Mi ey ep? (6,8) 

Wee eae U, ds] U;, 

so erhalt man nach Einfiihrung von ~ 
__ 9s (6,9) 
Fs aes 
ep os 6,10 
ee dy 8 ey, (6,10) 
ih U6, | UE dy U, 2 
0 


FuBnote 1 von S. 199. oe ; 
Gl. (6,9) und (6,19) der in FuBnote 1 von S. 195 zitierten Arbeit. 
Gl. (6,28) der in FuBnote 1 von §. 195 zitierten Arbeit. 

A. Kehl, Ing.-Arch, 13 (1943), S. 293. 
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und Multiplikation mit 2/U} 
d 6, dU 2D ee 
by (Hy — T)= (Hy — 7) (Bis 1) + (Ha — DF) (610 
Die Auswertung ist also bei Keilstrémung von der Gréfe des Keilwinkels unabhangig; daher 
haben auch Sekundarstrémungen in der Grenzschicht in erster Naherung keinen Einflu8 auf das 
Auswerteergebnis. 
Setzt man in (6,11) D nach (6,3) unter Beriicksichtigung von (6,4), (4,9) und (3,6) ein, und 
lést nach der gesuchten GréBe « auf, so erhalt man mit x —0,4 und A, =—0,9, wenn sich die 
Auswertung der Grenzschicht auf das Interval a <% <b erstreckt, im ersten Fall 


eee aC sx EP 6,12 
(De — 0,750 J, + 0,646 Jy + 0,2083 Jy (0,12) 


und im zweiten Fall 


(Hse — T)a 3,608 1,050 0104 Jf, -— 
‘in =I | ee | este ae ae 


= 0/563); 2 0,5 Lapa, oo. 


In 


(2) "a == (6,13) 


Hierbei bedeuten 


i= | 2 Galore ar (6,14) 


ii Ne 2» dx 
lire [las see ta 


Die Auswertung von 15 verschiedenen 
MeBreihen! nach dieser Methode ergab 
im Mittel fiir den ersten Fall 


(1) «= 2,57 mit einer mittleren Streu- 
ung von +0,45 


und fiir den zweiten Fall 


(2) «= 4,53 mit einer mittleren Streu- 


Abb. 5. GréBe G des Dissipationsintegrals in Abhangigkeit des 
Formparameters J,. ung von -+-0,65. 

Kurve 1: Rechnungsannahme (1): « = 2,57 + 0,45 | (6; 
# 


= 5,2; 
ape Diese Werte stimmen gréSenordnungs- 


Iwurve 2: Rechnungsannahme (2): «0 = 4,53 - 0,65 {A 

Kurve 3: Friiher vorgeschlagener Verlauf. i : mabig mit der Abschatzung nach dem 
Versuchspunkte nach friiherer Auswertung. Misch t bei d : ass Al 
Schraffierte Gebiete stellen mittlere Streubereiche dar. ischungswegansatzZ, bel der wir & — 


fanden, iiberein. 
Den mittleren Verlauf der Funktion G nebst ihrem mittleren Streubereich tiber dem Form- 
parameter I, zeigt Abb.5 fiir C=5,2. Zum Vergleich sind die Versuchspunkte der alten Aus- 


wertung und die damals vorgeschlagene Kurve mit eingetragen. Die neuen Kurven liegen etwas 


' Ks handelt sich um z.T. unveréffentlichte MeBreihen aus dem MPI f. Strémungsforschung und um 
MeBreihen der Schrifttumsstellen, FuBnote 3, S. 198; FuBnote 4, S.203 und A, E.v. Doenhoff und N. Tetervin, 


GLa ial of General Relations for the Behavior of Turbulent Boundary Layers, NACA Report No. 772 
1943). 
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tiefer, bestatigen jedoch die fritheren Auswertungsergebnisse im groBen und ganzen. Die Streu- 
ungen sind aber auch bei dieser Methode erheblich. Nun sind die Sekundarstrémungscinfliisse 
bei der Auswertung zwar in erster Naherung eliminiert worden, doch kennt man nicht die GréBe 
des verbleibenden Restes, so da man nicht sagen kann, inwieweit die Streuung auf solche 
Ursachen zuriickzufiihren ist. Es ist aber auch denkbar, daB der Streubereich einfach durch die 
Annahme der Einparametrigkeit bedingt ist. 

Die Streuung ist bei kleinen [,-Werten gering und wachst mit I, an. Das bedeutet, daf die 
Rechnungen, die auf Grund dieser Unterlagen ausgefiihrt werden, fiir geringen Druckanstieg 
oder Druckabfall ziemlich zuverlassig sind, bei groBen Druckanstiegen nimmt die Genauigkeit ab. 
Dies ist rechnerisch dadurch bedingt, da8 sich Fehler in a bei steigendem I, in wachsendem Mae 
auswirken, und physikalisch dadurch, da®B mit Annaherung an die turbulente Abldsestelle die 
Voraussetzung fiir die Giiltigkeit der einparametrigen Zusammenhinge weniger gut erfiillt ist, 
wie schon in Abschnitt 2 hervorgehoben wurde. 


7. Beriicksichtigung der Wandbindung. Bei den vorhergehenden Untersuchungen hatten wir 
die auf Grund der Navier-Stokesschen Bewegungsgleichungen geltenden Wandbindungen (3,1), 
fallen gelassen und auch Griinde dafiir angegeben. Wenn wir nun wenigstens die erstgenannte 

: Ot Op 

wd erage ey) 
doch beriicksichtigen wollen, gleichzeitig aber auch an unseren iibrigen Annahmen (Einpara- 
metrigkeit der Geschwindigkeitsprofile, eindeutige Abhangigkeit der Verteilung der statistischen 
Schwankungen von den Geschwindigkeitsprofilen) festhalten wollen, so bekommt man einen 
eindeutigen Zusammenhang zwischen Geschwindigkeitsprofil und dem in geeigneter Weise dimen- 
sionslos gemachten Druckgradienten dp/ax. Obwohl sich eine solche gegenseitige Zuordnung 
fiir die Berechnung turbulenter Grenzschichten nicht als sehr erfolgreich gezeigt hat, ist eine 
diesbeziigliche Untersuchung im Hinblick auf eine Nachpriifung der hier zugrunde liegenden 
Uberlegungen recht aufschluBreich. 

Nach Einfithrung der Ansatze (3,5) und (3,6) in (7,1) ergibt sich 


1dp,. aa 2 
= qo = 0 (A— A). (7.2) 
Mit Hilfe der Beziehung? 
1 
== A 

St oD 
“Ou* =e n (7,3) 

wird hieraus die GréBe 6 eliminiert, und man erhalt 

i 
1+ A 

es 2 7,4 
diggin: uae adc A,) 5 G (7,4) 


Bei der GréBe A handelt es sich um einen Parameter, der in ahnlicher Weise bereits von A. Buri? 
verwandt wurde. Unter Beriicksichtigung der Wandbindung (7,1) als zusatzliche Bedingung 
148t sich der ,,Buri-Parameter“ A also in Abhangigkeit von A und somit der Zusammenhang mit 
dem Formparameter I, nach (4,14) angeben. In Abb.6 ist dieser Zusammenhang fiir Werte 
a—2,57 und a=4,53 dargestellt. Gleichzeitig sind die aus MeBergebnissen mehrerer Versuchs- 
reihen von H. Ludwieg und W. Tillmann? errechneten Werte J, und A eingezeichnet. 

Wenn die vorgegebene Druckverteilung derart ist, da% der Buri-Parameter / itber eine langere 
Strecke in Strémungsrichtung konstant bleibt, so sollte man erwarten, daB sich das Grenz- 
schichtprofil etwa so einstellt, wie es die Theorie unter Beriicksichtigung der Wandbindung (7,1) 
voraussagt. Tritt dann eine ploétzliche Anderung von A ein, so wird das Geschwindigkeitsprofil 
- diese Anderung nicht sofort mitmachen; die Anderung des Geschwindigkeitsprofiles hinkt nach. 
Hat man nun Druckverteilungen, bei denen J in Strémungsrichtung stetig ansteigt, wie es bei 
den beiden, in Abb.6 durch offene Kreise und gerade Kreuze gekennzeichneten MeBreihen der 
Fall ist, so sollte man demnach fiir I, Werte erwarten, die unterhalb der theoretischen I,, A-Kurve 


1 Gl. (6,15) der in FuBnote 1 von 5.195 zitierten Arbeit. 
2 FuBnote 2, 5. 197. 
3 H. Ludwieg u. W. Tillmann, Ing.-Arch. 17 (1949), S. 288 sowie unveréffentlichte Versuchsreihen. 
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liegen. Die beiden erwahnten MeBreihen liegen in Abb.6 zwischen den beiden rechnerischen 
Kurven, folgen aber qualitativ ihrem Verlauf. Hiernach erscheint die Kurve 1, deren Annahmen 
schon im Abschnitt 3 als die Wirklichkeit besser treffend erkannt wurden, die richtigere von 
beiden zu sein. 

Diese Gedanken werden noch durch die durch schrage Kreuze und Dreiecke dargestellten 
Versuchsreihen erhartet, bei denen J in Strémungsrichtung zundchst ansteigt und dann wieder 
abfallt. Solange J anwiichst, liegen die Auswertungspunkte wieder unterhalb der Kurve J. Wenn 
dann A aber wieder abnimmt, steigt I, zunachst weiter und iiberschneidet die Kurve 1. Eine 
Ausnahme bildet die durch ausgefiillte Kreise gekennzeichnete Reihe mit nur schwachem Druck- 
anstieg. Hierbei bleibt I, auch nach Uberschreiten des Groéftwertes von A etwas unterhalb der 
Kurve I. 


tl; 
IO ] 


05 


-50 0 50 700 750 200 250 300 350 400 


Abb. 6:. Zusammenhang zwischen Formparameter I, und Buri-Parameter A 
Kurve 1:3 & = 2,07 \ # = 0,4; 
Kurve 2: 4 = 4,53 f A, = 0,9. 
Versuchspunkte nach H. Ludwieg und W. Tillmann. 


Die Deutung der Versuchsbeobachtungen ist mit den Uberlegungen der vorhergehenden 
Abschnitte vertraglich, so daB diese Untersuchung als ein Beitrag zur Erhartung der angegebenen 
Zusammenhange gewertet werden kann. Dariiber hinaus gibt aber Abb.6 noch die Méglichkeit, 
die Brauchbarkeit eines vereinfachten Berechnungsverfahrens zu beurteilen, bei dem der Form- 
parameter als eindeutige Funktion von A aufzufassen ware. Ein solches Verfahren entspricht 
dem von A, Buri vorgeschlagenen und bietet den Vorteil, da fiir die Berechnung der Grenz- 
schicht nur eine gewohnliche Differentialgleichung statt deren zwei zu integrieren ist. Nach 
Abb. 6 scheint ein solches Verfahren durchaus brauchbar zu sein, wenn es sich um mabige Druck- 
anstiege handelt. Bei gréBReren A-Werten mit wechselndem Verlauf versagt allerdings ein solches 
Verfahren. 


8. Zusammenfassung, Zweck der Arbeit ist es, den Zusammenhang zwischen dem bei der 
Berechnung turbulenter Grenzschichten mit Hilfe des Impuls- und Energiesatzes benétigten 
Dissipationsintegral und dem Formparameter des Grenzschichtprofiles zu ermitteln. Als Grund- 
lage wird die Hypothese benutzt, daB® sich das Geschwindigkeitsprofil (mit Ausnahme der von 
der Zahigkeit und Wandrauhigkeit beeinfluBten, diimnen Schicht nahe der Wand) als einpara- 
metrige Kurvenschar darstellen la8t, und da® zu jedem Geschwindigkeitsprofil auch ganz be- 
stimmte Verteilungen der turbulenten Schwankungsgeschwindigkeiten gehéren. Die Anwend- 
barkeit dieser Annahme setzt voraus, da die Reynoldszahl hinreichend gro8 und der Grenz- 
schichtzuwachs dd/dx von zweiter Ordnung klein gegen 1 ist. 
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Fitr die Schubspannungsverteilung wird dann ein Naherungsansatz eingefiihrt, mit dessen 
Hilfe unter Verwendung einer schon friiher benutzten Formel fiir das Geschwindigkeitsprofil die 
Berechnung des Dissipationsintegrals bis auf eine empirische Konstante méglich ist. Die Gre 
dieser Konstanten wurde aus einer Anzahl gemessener turbulenter Grenzschichten berechnet. 
Der Fortschritt gegeniiber einer fritheren Auswertung ist dabei, da® jetzt Quadraturen an die 
Stelle der frither notwendigen Differentiationen treten. 

Eine Untersuchung iiber den Einflu8 der Wandbindung é1/dy = dp/dx, deren Beriicksich- 
tigung theoretisch in ahnlicher Weise wie bei dem von A. Buri vorgeschlagenen Berechnungs- 
verfahren eine Beziehung zwischen dem Druckgradienten dp/dx gebildeten Parameter liefert, 
erganzt die Untersuchungen. 


(Eingegangen am 20. Dezember 1951.) 
Aus dem Max-Planck-Institut fiir Str6mungsforschung, Gottingen. 


Anschrift des Verfassers: J. Rotta, Géttingen, Béttingerstr. 6/8. 


Zur Berechnung des EinfluBfeldes der quadratischen Platte. 
Von E. R. Berger. 


In seiner Arbeit ,, Uber die Singularitatenmethode an elastischen Platten‘‘! berechnet Pucher 
die KinfluBfelder der allseitig eingespannten quadratischen Platte nach der Singularitatenmethode, 
mit Erfiillung der Randbedingungen an einzelnen Punkten. Unter anderem rechnet er (S. 87) das 
Einflu8feld des Drillungsmomentes in der Mitte durch den Ansatz 


w =~ (sin p + 3) 4,0?" sin 2ng + 3b, 0?"*? sin 2n¢—y). (1) 


Der Ansatz liefert mit sechs Freiwerten zufriedenstellende Ergebnisse. 
Als Nachstes berechnet Pucher das Stiitzmomenten-EinfluBfeld in der Seitenmitte (S. 88) und 
macht dafiir den Ansatz 


i — + (= cos? pm + S Cn, 0" +! cos @ cos ng) (2) 
n=1 


Er schreibt dazu: .,Die guten Erfahrungen, die mit der beim Drillungsmomenten-EinfluBbfeld 
angewendeten Methode gemacht wurden, lieBen erwarten, daB man auch hier auf ahnliche Weise 
zum Ziele kommen wiirde ... Das Ergebnis war aber unbrauchbar, obwohl die Summe iiber elf 
Funktionen, also verhaltnismaBig viel Glieder, erstreckt wurde. Die Nachpriifung zeigte, daB.. . 
zwischen den Punkten, in denen die Randbedingungen streng erfiillt sind, sowohl] die Funktions- 
werte als auch die Neigung der Randtangenten viel gréBere Abweichungen vom Sollbetrag zeigten, 
als zugelassen werden konnte. Die Ursache ist in der mit steigendem n gréBer werdenden Anzahl 
der Nullinien der einzelnen Funktionen zu suchen, so da’ auch die Hinzunahme weiterer Glieder 
das Ergebnis nicht verbessern konnte.“‘ 

Die hier beschriebene Erscheinung ist ein typisches Kennzeichen dafiir, dafs die im Ansatz ver- 
wendeten Funktionen kein ,,vollstandiges‘‘ Funktionensystem bilden. Sie wiirde in genau der- 
selben Form auftreten, wenn man z. B. eine periodische Funktion in einer ungeraden Anzahl (n) 
von aquidistanten Punkten durch ein trigonometrisches Polynom exakt darstellen will, dabei aber 
im Ansatz nur die Harmonischen gerader Ordnung verwendet. Anstelle der im Ansatz fehlenden 
Funktion cosx z. B. tritt jetzt mit derselben Amplitude die Funktion cos (n — 1) x auf, da sie in 
den vorgegebenen Punkten mit cos x itbereinstimmt. Zwischen je zwei solchen Punkten durch- 
lauft sie aber fast eine ganze Periode, daher kann von ciner Annaherung in den Zwischenpunkten 
keine Rede mehr sein. 

Etwas Ahnliches ist nun offenbar auch beim Ansatz (2) geschehen. Zunachst ist aus dem Text 
nicht ersichtlich, ob unter nin (2) alle natiirlichen Zahlen oder nur die ungeraden Zahlen verstan- 
den sein sollen. Fiir gerade n erfiillen aber die einzelnen Glieder zwar die Randbedingung 


w= 0 fir o=--an/2, (3) 


1 [ng.-Archiv, 12 (1941), S. 76. 
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sie erfiillen aber nicht die Randbedingung 


M0 fir p= ta/2. (4) 
Wir miissen also die geraden n ausschlieBen, und das hat Pucher wahrscheinlich auch getan, da er 
vorher schreibt: ,,Es liegt nahe, auch das allgemeine Integral so anzusetzen, daB es, unabhangig 
von den Freiwerten, am Rande gleichzeitig mit der Ableitung verschwindet.‘* Nun sind aber auch 
die geraden Zahlen Eigenwerte des Problems, und da die zugehérigen Eigenfunktionen im Ansatz 
fehlen, ist das System eben unvollstandig. 


Die Eigenlésungen fiir die geraden Zahlen lauten 
all 
= Anal ay == 
w, = 0 [cos p eos NP sin p sin ng) (n = 2,4,6,...). (5) 
Wir haben daher fiir die Lésung anzusetzen 


il 
w = Z| — oot + S  Cn410"t* cosy cos np + | 
wt nals 8555 cc 


(6) 
SS Opes ont (cos g cos n@ — sing sin no) . | 
4, 


Te oe a Orerere 
Die Rechnung wurde versuchsweise mit diesem Ansatz durchgefiihrt und die Randbedingungen 
in den Punkten A, B und C (vgl. Abb. 1) streng erfiillt. A und C geben je zwei Randbedingungen, 
B drei; wir benétigen also 7 Freiwerte. Das numerische Ergebnis wird weiter 
unten mitgeteilt. 

Im Grunde genommen hat Pucher selbst schon die Funktionen (5) ver- 
wendet, nur in einer anderen Gestalt. Er schreibt namlich weiter: ,,Es 
muften andere Funktionen herangezogen werden, die langs x = 0 dem ein- 
gespannten Rande entsprechen. Auf der Suche nach solchen Funktionen 
wurden Polynome gefunden ...*‘ Diese biharmonischen Polynome (bzw. 
jene darunter, die in y gerade sind) sind aber genau identisch mit den 
Funktionen des Ansatzes (6). Das 1aBt sich im Einzelnen durch algebraische 
Umformung verifizieren, geht aber auch aus folgender Uberlegung hervor: 


Die biharmonischen Polynome sind homogen in x und y, kénnen also auch in der Form ge- 
schrieben werden 


n+l 
Pe 02 DN a, Cost pein Ot) ae 


Das la8t sich umformen in 


n+1 
Py = 07th, cosk@ 


(baw. fiir die in y ungeraden in Yb, sin kg). Da die biharmonischen Polynome Lésungen sind 
von AA P = 0, scheiden alle Glieder aus au®er den beiden 


Pen OCG (Gr cos (n + 1) + €, cos (n — 1) p)= 0"! (d, cosp cos ny + d, sin y sin nq). 


Die Randbedingungen (3) und (4) bilden zwei homogene Bestimmungsgleichungen fiir d und d, die 


aber vertraglich sind (weil eben eine Eigenlésung vorliegt), sie bestimmen das Verhaltnis d:d, und 
damit ist P bis auf einen konstanten Faktor festgelegt. Da nun die Funktionen des Ansatzes (6) 
durch genau dieselben Bedingungen definiert sind, kénnen sie sich von den P nur durch einen kon- 
stanten Faktor unterscheiden. Setzt man etwa x — 1, y = Obzw. 0 = 1,y = Oin die beiden Aus- 
driicke ein, erhalt man jedesmal 1, womit die Identitat bewiesen ist. 


Zur Bezeichnung der biharmonischen Polynome wire noch zu bemerken: es wire vielleicht 
zweckmabig, die Zeiger nicht nach dem Grad der einzelnen Variablen, sondern nach dem Grad des 
ganzen Ausdruckes zu wahlen, und die in y geraden bzw. ungeraden Funktionen durch verschie- 
dene Buchstaben zu unterscheiden, also etwa: 
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Pearls 
Epa PO COSC), 
P= "==. p* Cos" @), 


0; = ay =o cosy sing; 
es : 
Pea = 08 (cos p cos 29 + sin sin 29), (7) 
Q= y= 5 cosy sin29, 
Pi= x4 —3x2¥% = 0 cosgp eos 3Q, 


Opa ye ot: : (cos p sin 3 — sin | cos 39) : 
Allgemein ist 
Py, =o" cosy cos (2n —1)q, 


2n— I] , 1 : 
O35 = OF aay [coe sin (2n— 1) ony Sin Y Cos (Qn —1)¢] . 
ae [eos p cos 2n@ +5 sing sin 20] ; Cy 


1 ; 
Ore = 0777 5 Cosy sin 2n@p 


(ausgenommen fiir P,) P, und Q,). 


Fir die Durchrechnung der quadratischen Platte wahlt Pucher die Lange der halben Quadrat- 
_seite mit 0,4. Er schreibt dazu (S. 91): ,,Infolgedessen bleiben diejenigen Bereiche, in denen bei 
den Polynomen héheren Grades die Tangente 0 P;/0x bzw. 0 P,/dy sehr steil ansteigt, auBerhalb des 
Bereiches. Das ist zweckmafig, weil auch die aus der Singularitat folgenden Randneigungen 
innerhalb kleiner Grenzen bleiben und somit die Angleichung der Randwerte des Regularteiles an 
die des Singularteiles erleichtert wird.“* Diese Wahl des Ma®istabes kann nun vielleicht die prak- 
tische Durchfiihrung der Rechnung etwas erleichtern, hat aber auf das Ergebnis prinzipiell keiner- 
lei Einflu8. (Das wird grundlegend anders, wenn man Integrale verschiedener Dimension addiert, 
wie es Pucher im nachstfolgenden Abschnitt tut [S. 93, Formel (27)]. Andert man dort den MaB- 
stab, so Andert man das Gewicht der einzelnen Integrale und damit auch das Ergebnis). Da die 
einzelnen Summanden homogen sind, bewirkt die MaBstabainderung lediglich, daB die Koeffizien- 
ten c, sich durch (0,4)" dividieren. Man kann also am SchluB der Rechnung ohne Miihe wieder 
auf das Quadrat mit der halben Seitenlange 1 zuriickgehen, indem man die erhaltenen c, mit 
(0,4)" multipliziert. Mit den so reduzierten Koeffizienten von Pucher wurde das Ergebnis der 
oben angefiihrten Rechnung verglichen: 


Wie man sieht, stimmen die ersten ae Pane 
Perdensueniaenten, uusciahr ubeL- 

ein, im weiteren Verlauf wird der C5 + 1,159806 + 1,1658293 
Unterschied immer gréfer. Op — 0,489015 — 0,5147134 
= : im : Cy — 0,277883 — 0,1346960 
Fir die gréBten Abweichungen a +. 0.228991 +. 0.1278343 
in den Zwischenpunkten ergibt sich C6 + 0,019572 — 0,0154916 
(vgl. S. 92): Cy — 0,006522 — 0,0082812 
Cy + 0,001844 + 0,0022277 

|Aw| << 0,036/z bzw. 0,036/a F Cy aS 0,0003407 oa) 

é ; c — 0,0002643 — 

fiir die Abweichung der Steigung on +. 0,00002994 io 

(reduziert auf Seitenlange 1): oe + 0,000002798 2s 


4 on <0,29/m baw. 0,44/7 


schlieBlich fiir das Stiitzmoment bei Gleichlast: 
—0,0519 pi? bzw. — 0,0484 pl. 


Die Genauigkeit des 7-gliedrigen Ansatzes wiire also noch unzureichend. 
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Zusammenfassend kann man also feststellen: Der von Pucher verworfene Ansatz ist nicht 
prinzipiell ungeeignet, sondern nur unvollstandig. Die von Pucher verwendeten biharmonischen 
Polynome sind mit dem vervollstandigten ersten Ansatz identisch. Die Wahl des Mafistabes ist 
fiir das Ergebnis ohne EinfluB. 


(Eingegangen am 25. Januar 1952.) 


Anschrift des Verfassers: Dr. E. R. Berger, Wien, Technische Hochschule. 


Zuschrift 
zuder Arbeit von E. R. Berger, ,,Zur Berechnung des HinfluBfeldes 


6G 


der quadratischen Platte”. 
Von A. Pucher. 


Die Aufklarung, die Herr Berger nunmehr iiber die Ursache des Versagens des von mir seiner- 
zeit verwendeten Ansatzes 


= aS cos?g + Cn Ot! cosy cos ng) (1) 


W=1, 35 
gebracht hat, ist sehr zu begriiBen, da sie unsere Kenntnisse iiber das vollstandige Integral von 
AAw = 0 vertieft hat. 

Die Identitaten der biharmonischen Polynome mit den trigonometrischen Reihen 

n+1 n+1 
P,,. 0") Dia, vos* pein 2) 0" tee Deosiye 
wurden schon von Herrn K. Zweiling in den Jahren 1942—43 im Rahmen einer umfassenden 
Arbeit iitber die biharmonischen Polynome gefunden und mir mitgeteilt. Leider ist keine Ver- 
éffentlichung dieser wertvollen Ergebnisse erfolgt. 

Ich habe mich mit dem Versagen des Ansatzes (1) nicht weiter beschaftigt, da ich zunachst 
mit Hilfe der biharmonischen Polynome, die in kartesischen Koordinaten ausgedritckt, fiir die 
Berechnung rechteckiger Platten naturgemaB viel bequemer sind, als in Polarkoordinaten dar- 
gestellte Funktionen, die gesuchten EinfluBfelder mit ausreichender Genauigkeit berechnen 
konnte. Im Verlaufe meiner weiteren Untersuchungen fand ich jedoch die ,,modifizierten Singu- 
laritaten“, die erstmalig 19431, sowie 1950? und 19513 mitgeteilt wurden, und die die Berechnung der 
KinfluBfelder noch viel einfacher gestalteten. Mit Hilfe dieser modifizierten Singularitaten ist es 
mdglich, die verschiedenen EinfluBfelder der Rechteckplatten mit so gut konvergierenden Reihen 
zu berechnen, daf in den meisten Fallen nur ein Reihenglied zu beriicksichtigen ist. 

Demnach sind die Ausfiihrungen von Herrn Berger vorwiegend fiir die theoretische Erkenntnis 
wertvoll, fiir die numerische Ermittlung der EinfluBfelder der Rechteckplatten jedoch weniger 
von Bedeutung. 


(Eingegangen am 3. April 1952.) 
Anschrift des Verfassers: Prof. Dr.-Ing. A. Pucher, Graz, Burgring 16. 


1 A. Pucher, Ing.-Arch. (1943), S. 246. 

* A, Pucher, Die EinfluBfelder des Plattenstreifens mit zwei eingespannten Randern; Beitrige zur 
angewandten Mechanik, Federhofer-Girkmann-Festschrift. Wien 1950. 

8 A, Pucher, KinfluBfelder elastischer Platten, Wien 1951. 
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Magnetische Werkstoffe. Von Dr, techn, Ing. Franz Pawlek, o. Professor an 


der Technischen Universitat Berlin-Charlottenburg. (Reine und angewandte Metallkunde 
in Einzeldarstellungen. Herausgegeben von W. K éster, Band 11.) Mit 270 Abbildungen. 
VII, 303 Seiten. 1952. SS Ganzleinen DM 42.— 


Inhaltsibersicht: Einleitung. Dauermagnete. Allgemeines: Charakterisierung der Dauermagnete. 
Theorie der Koerzitivkraft. Theorie der Remanenz, —Dauermagnet-Werkstoffe: Geschichtlicher Uber- 
blick. Durch Gefigeumwandlung gehartete Werkstoffe. Durch Ausscheidung gehiartete Werkstoffe, Legierungen, 
deren Dauermagneteigenschaften mit der Ausbildung einer Uberstruktur verknipftsind. Pulvermagnete. Oxydmagnete, 
Magnetisch weiche Werkstoffe. Allgemeines; Charakterisierung der magnetisch weichen Werkstoffe. Theorie 
der Remanenz. Theorie der Koerzitivkraft und Anfangspermeabilitét. — Magnetisch weiche Werk- 
stoffe: Geschichtlicher Uberblick, Reines Hisen., Technische Eisensorten, Lisenlegierungeg, Eisen-Nickel- 
Legierungen. — Magnetische Werkstoffe fiir Pupinspulen und Hochfrequenzkerne. All gemeines: Charakteri- 
sierung der Werkstoffe fir Pupinspulen und Hochfrequenzkerne. — Massekerne: Eigenschaften der -Masse- 
kerne. Herstellung und Eigenschaften der Massekerne. — Iso perme: Ausscheidungsisoperm, Textur-Isoperm, 
Plattierungs-Isoperm. — Nichtmetallische magnetische Werkstoffe: Theorie des Ferrimagnetis- 
mus. Higenschaften der Ferrite. Verwendung, — Werkstoffe mit groBer Magnetostriktion. — Werkstoffe mit be- 
Sonderen physikalischen Eigenschatten auf Grund verborgener magnetischer Vorgange. Allgemeines: Theorie 
der Werkstoffe mit geringem Temperaturkoeffizienten der Ausdehnung. Theorie der Werkstoffe mit geringem 
Temperaturkoeffizienten des E-Moduls. —-WerkstoffemitkleinemAusdehnungskoeffizienten: 
Invarlegierungen. Glas-Einschmelz-Legierungen. — Werkstoffe mit kleinen Temperaturkoeffi- 
zientendes E-ModuIs. — Magnetische Werkstoffe mit stark temperaturabhangiger Sattigungsmagnetisierung. 
_ Unmagnetische Werkstoffe: Nichteisenmetalle. Unmagnetisches GuBeisen, Unmagnetischer Stahl, — Anhang. — 
Namen- und Sachverzeichnis. 2 5 rae 


Die Edelstiihle.. Von Professor Dr.-Ing. Franz Rapatz, Stahlwerk Gebr, Béhler 
& Co. A.-G. Kapfenberg (Steiermark). Vierte, verbesserte und erweiterte Auflage. Unter 
Mitwirkung von Dr.-Ing. Helmut Krainer, und Dipl.-Ing. Josef Frehser, Stahlwerk Gebr. 
Béhler & Co. A.-G., Kapfenberg (Steiermark). Mit 338 Abbildungen und 121 Zahlentafeln. 
V, 730 Seiten. 1951. ; Ganzleinen DM 49.50 


Ausden Besprechungen:... Das Buch ist einfach und klar geschrieben, und alle Spezialausdriicke werden 
auf das genaueste erklart. Wenn der Leser einigermaBen gute Kenntnisse in Physik und Chemie besitzt, so vermag 
er ohne Schwierigkeit der Darstellung zu folgen. Wenn auch ein unerhért groBes Wissen auf relativ engem Raum 
zusammengedringt ist, so ist das Buch doch leicht lesbar. Es ist zweifellos ein besonders’ wertvolles Nachsehlagewerk, 
in dem man praktisch genommen alles theoretisch und praktisch Wissenswerte auf diesem Gebiet wenigstens er- 
wahnt findet. Ein reichhaltiges Literaturverzeichnis macht es dem Leser méglich, ohne besondere Miihe die Angaben 
des Buches zu ergainzen. Es reicht bis ungefahr zum Jahresweehsel 1949/50 und umfaBt nicht nur Deutschland, 
sondern auch USA und in gewisser Hinsieht GroBbritannien... »leknisk Tidskrift**. 


Die Zerspanbarkeit der metallischen und nichtmetallischen 


Werkstoffe. Von Dr.-Ing. habil. .Karl Krekeler, Professor an der Technischen 
Hochschule Aachen. Mit 148 Abbildungen. XII, 358 Seiten. 1951. Ganzleinen DM 34.50 


hungen:... Krekelers neue, wirklich umfassende Arbeit stellt — ohne jedoch praktische 
Foteoinanen a abemeken — yerstindlicherweise mehrheitlich auf wissenschaftliche _Untersuchungs- 
ergebuisse ab und arbcitet diese, die Erfahrungen anderer Zerspanungsspezialisten und -forscher mitverwertend, in 
einer fir die Praxis direkt brauchbaren Form aus. Weit mehr Hand- als Lehrbuch ist dasselbe dem Zerspannngse 
fachmann ohne jegliche Beschrankung zu empfehlen... — » Technische Rundschau“, 


Hydraulische Schmiedepressen und Kraftwasseranlagen. Kon- 


struktion und Berechnung. Von Ernst Miller, Duisburg. Zweite, verbesserte 
Auflage. Mit 167 Abbildungen. VII, 208 Seiten. 1952. . Ganzleinen DM 28.50 


i i n Schmiedepressen stehen neben Walzwerken und Himmern unter den Maschinen zur spanlosen 
Tete, to Meas an eaten Stelle. Sie haben im Laufe de® Entwicklung, hauptsichlich nervorgerufen durch 
den Bedarf an nahtlosen, geschmiedeten Hochdruckkesseln und immer gréBer werdenden’ Schmiedestiicken fir den 
Schiffs- und Maschinenbau sowie fiir die chemische Industrie, Abmessungen erhalten, die sie in die Reihe der groBten 
Bauwerke der -Maschinenindustrie stellen. Die erste Auflage dieses Buches erschien kurz vor Ausbruch des zweiten 
Weltkrieges. In der Zwischenzeit ist die Entwicklung der Konstruktionen weitergegangen, die sich im wesentlichen 
auf die Vervollkommnung der Antriebe und Hilfseinrichtungen fiir die Pressen beschrinkte. Dampf und Luft als 
Betriebsmittel fir hydraulische Schmiedepressen wurden immer mehr zugunsten des rein- und elektro-hydraulischen 
Betriebes aufgegeben. Wahrend in England fiir schwere Pressen vornehmlich der direkte Pumpenantrieb weiter- 
entwickelt wurde, beyorzugte man in den USA und Deutschland den reinhydraulischen Antrieb in Verbindung mit 
Druckiuftakkumulatoren. Auf den letzten Seiten des Buches werden aus allen Lindern der Schwerindustrie die 
gréBten SchmiedepreBwerke aufgefihrt, wobei auBer dem Druckvermégen der Pressen auch die Betriebsart und das 
Baujahr angegeben werden. GroBe Fortschritte sind auf dem Gebiete der Mechanisierung des Blocktransportes durch 
weitgehende Anwendung von Schmiedemanipulatoren zu verzeichnen, Sie hat zu groBer Produktionssteigerung und 
Einsparung von Blockwirmen und Arbeitskraften gefiihrt. In der vorliegenden neuen Auflage sind die einzelnen 
Kapitel dem heutigen Stand der Technik angepaBt und entsprechend erweitert worden. 
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Statische Berechnung von Kesselboden. von p:.-ing. Maria EBlinger, 
Saarbriicken. Mit 21 Abbildungen. VIII, 100 Seiten. 1952. DM 10.50 


Aus den Besprechun gen: ... Auf Grund der aufgestellten Gleichungen und Rechen- 
schemata sind eine Reihe von Beswelen durchgerechnet, von denen drei solchen Béden ent- 
sprechen, die Siebel experimentell untersucht hat. Die Ubereinstimmung zwischen Rechnung 
und Versuch ist — unter Beriicksichtigung der méglichen und unvermeidlichen Mefifehler bei 
solchen Spannungsmessungen an Hohlkérpern — so iiberraschend gut, daB man den Rech- 
nungsergebnissen und den Schluffolgerungen, die EBlinger daraus zieht, volles Vertrauen ent- 
gegenbringen kann. Das Rechenverfahren liefert daher ein wertvolles Hilfsmittel, um in 
schwierigen Fallen eine sichere Entscheidung tiber die erforderliche Wanddicke eines Bodens 
treffen zu kénnen, und eine wertvolle Bestatigung fiir die bewahrten, auf Grund der Versuche 
von Siebel u, a. aufgestellten Gleichungen der Werkstoff- und Bauvorschriften fir Land- und 
Schiffsdampfkessel. >» Hnergie’. 


SPRINGER-VERLAG / BERLIN- GOTTINGEN - HEIDELBERG 


Berechnung mechanischer Schwingungen. von Dipi-Ing., Dr. techn. 


Fritz Séchting, iit. a. o. Professor an der Technischen Hochschule Wien. Mit 140 Text- 


abbildungen. X, 325 Seiten. 1951. -Ganzleinen DM 32.70 — 


Die standig wachsenden- Anspriiche an die Maschinenleistung und das Bestreben, durch immer 
knappere Abmessungen Gewieht, Raum und Werkstoff einzusparen, haben die Berechnung 
der an Maschinenteilen auftretenden Schwingungen und Resonanzfrequenzen zu einer wichtigen 
Aufgabe fiir den Techniker gemacht. Das vorliegende Buch wendet sich ebenso an den Stu- 
dierenden, dem es den Weg zur selbstandigen Weiterarbeit weisen will, wie an. edi in der 
_ Praxis stehenden Ingenieur. 


- 


Mechanische Technologie. von vipi-tng. A. Kopecky, Steyr, und Dipl.-Ing. 
Dr. techn. R. Schamschula, Wien. Mit 424 Textabbildungen. X, 363 Seiten. 1951. 
DM 18.—; Ganzleinen DM 19.80 


Aus den Besprechungen: ... Auf verhaltnismaBig kleinem Raum kann anf einem weit- 
gesteckten Gebiet wie dem der mechanischen Technologie nur ein Wissen vermittelt werden, 
wenn in der Tiefe der Darstellung Beschrankung geiibt wird. Dies ist den Verfassern ge- 
lungen. Je mehr sich die industrielle Fertigung bemitht, die Endform in voraufgehende Ar- 
beitsgange vorzuverlegen, wird sich der Betriebsmann und Konstrukteux mit den Eigenarten 
und Méglichkeiten dieser Arbeitsverfahren vertraut machen miissen, Er wird sich dabei nicht 
nur auf die ,klassischen‘ Werkstoffe des Maschinenbaues beschranken kénnen. Dies haben 
die Verfasser beriicksichtigt. > lndustrie-Anzeiger**, | 
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